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The Band Spectrum of GeF 


By E. B. ANDREWS anp R. F. BARROW 
Physical Chemistry Laboratory, University of Oxford 


IMS. received 10th August 1949 


ABSTRACT. ‘Two band-systems attributed to GeF have been observed in emission 
from a discharge through a stream of GeF, vapour. Vibrational analyses of both systems 
have been made and the following constants derived : 


State Ve We XeWe (cm~*) 

ey 35007-4 800°4 4-1; 

2 ey, 23316°8 413-5 1-0; 
Aleve 935 665-8 2°85 
Alla 0) 664°5 DT 


A brief discussion of the results, and a comparison between the force-constants of the mono- 
and tetra-halides of the group IV—B elements, is given. 


SUN ROD LOIN 


may be expected that the study of the variation in bond properties in going 

from fluorine to iodine and from carbon to lead—both among themselves 
and in comparison with the tetrahalides—may be of some interest. Gel is one 
such molecule for which there appears to be no experimental information. It 
might be expected to give rise to an ultra-violet band-system, ?& —*I1, analogous 
to those of GeCl and GeBr (Jevons, Bashford and Briscoe 1937), and also to a 
system, perhaps ?A—#II, similar to that recently found for GeCl in the near 
ultra-violet (Barrow and Lagerqvist 1949). 

The present work shows that GeF does indeed emit two band-systems. One,. 
in the ultra-violet region, probably corresponds to the 7& —?IH systems of GeCl 
and GeBr; the second, in the visible region, has however no analogue in GeCl, 
but appears to be similar to the least refrangible system of SnF (Yuasa 1939). 
The lower state of both transitions is almost certainly the ground-state, *I]. For 
the *II5.—Iy. separation we find 935cm~1, in excellent agreement with the 
value of 940 cm“! predicted by Howell (1945). 
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T= diatomic halides of the group IV—B elements form a series in which it 


186 E. B. Andrews and R. F. Barrow 


§2. EXPERIMENTAL 

A simple quartz discharge-tube fitted with molybdenum wire electrodes was 
used. GeF, was produced by heating solid BaGeF, (Dennis and Laubengayer 
1927) contained in a side-arm. A discharge at about 0:25 amp., 1,000 v., was 
sufficient to produce an intense violet-coloured source in which the bands 
attributed to GeF were observed. 

The systems were photographed in about 20 minutes on a Hilger E1 spectro- 
graph with glass and quartz optics. One system, degraded to the violet, 
lies in the region AA 2600-3100a., the other has its strong bands in the region 
dA 3900-5200 ., and is degraded to the red. 

Interference due to other spectra was slight, and limited to SiF and, to a lesser 
extent, OH. GeO bands were observed on the short wavelengths side of the 
ultra-violet system. Ge, H, Mo and Si lines which appeared on the plates were 
used as standards in compucing wavelengths: their values were taken from the 


M.I.T. Wavelength Tables (1939). 


§3. APPEARANCE OF THE BANDS 


Repreductions of spectrograms of the two systems are given in the Plate. 


(i) The Ultra-violet System 


The system splits clearly into two sub-systems, so that each band consists of 
two sub-bands. Each sub-band possesses two heads, one P-like and the other 
Q-like. Near the origins the separation of these two heads is about 10cm}. 
The rotational structure is only partially resolved near the origins. Away from 
the origins, the heads become progressively more complex in appearance and more 
than two heads are observed for each sub-band. ‘The separations of these heads 
are of the order of magnitude expected for a vibrational effect arising from the 
Ge isotopes. 

(ii) The Visible System 
This is similar to the ultra-violet system, except that the bands degrade to 


the red. Near the origin, the R-like heads are separated from the Q-like heads ~ 
by about 5cm “1. | 


§4. VIBRATIONAL ANALYSIS 


No particular difficulties were encountered in the analyses, which were facili- 
tated by the strengths of the 0,0 bands and by the vibrational isotope effect. 
The data are set out in Deslandres schemes in Tables 1 and 2. In deriving the 
vibrational constants given in Table 3, only the figures for the Q-like heads were 
used: the v, values were got from these heads of the 0,0 bands by correcting for 
zero-point energy. 

A Birge-Sponer extrapolation gives 4:9 ev. for-the dissociation energy of the 
ground state. ‘This fits reasonably with the value of 4:3 ev. for GeCl. Similar 
extrapolations for the upper states lead to considerably higher dissociation limits, 
but they are much less certain and it is hardly possible to say more than that the 
higher excited state probably gives one or more excited atoms on dissociation. 
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§5. DISCUSSION 
There is little doubt that the bands arise from transitions between states of 
GeF. ‘The nature of the source, the appearance of the bands and the fact that the 
derived constants are of the order of magnitude expected all lead to this conclusion. 
Whilst the lower state is almost certainly the 2I] ground state, all that can 
definitely be concluded about the upper states is that both appear to show very small, 


Table 1. Band-head Data for Ultra-violet System of GeF 


(Figures in parentheses are visual relative intensities) 
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if not zero, doublet separations: to identify them as 22 would be consistent with 
Some support to the view that these excited 
states of GeF are 2X comes from a comparison with SnF. GeCl is remarkably 
similar in spectroscopic properties to SnCl, and we might expect the fluorides 
Of the five stable excited states of SnF known (Jenkins and 
Rochester 1937, Yuasa 1939), it appears that the two lowest excited states are 


the apparent structure of the bands. 


also to be similar. 
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probably analogous to the two states of GeF, thus: 


Both these states of SnF have been supposed to be 7%, so it appears that the 
similarity of the tin and germanium compounds extends at least partly to the 


e>o 


GeF 


Poy 
4:34 
2-89 


SnF 
RR’ Ve, EV. Rk” 
1:45 4-23 1-34 
0-39 2:50 0:51 


k is the force-constant. 
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Prism spectrograms of (a) ultra-violet system, (0, c) visible system, of GeF. 
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fluorides. It is therefore probable that transitions involving more highly-excited 
states of GeF at wavelengths less than 2500 A. await discovery. 

Finally it is of interest to compare the force-constants of the diatomic halides 
with those for the corresponding tetrahalides. On general grounds it might be 
expected that the valency bonds should increase in strength in passing from 
MX to MX,, because the hybridized orbitals in MX, are likely to be stronger 
than those in MX. ‘This is supported by the experimental values of the inter- 
nuclear distances in SiF (7) =1-60., Eyster 1937) and in Sif, (r= 1-54.a., Maxwell 
1940). However, as Heath and Linnett (1948 a,b) have emphasized, repulsion 
between non-bonded atoms is important, particularly of course when the central 
atom is small, and the derivation of the bond-stretching constant k, from the totally 
symmetric vibration vy, becomes dependent upon the law of interaction between 
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the X atoms. In the table below are given for k,(MX,), (1) the values calculated 
with the simple valency force field, ‘S.V.F.F.’ (Herzberg 1945), (ii) the values 
calculated by Heath and Linnett (1948b) using the orbital valency force field, 
‘O.V.F.F.’, and including the effect of repulsion between non-bonded atoms. 

It is feed possible to conclude from these figures alone which set of force- 
constants for the tetrahalides is the more likely to be correct. But either set shows 
what seems to be a significant decrease in the ratio k(MX,)/A(MX) in the direction 
F to Br for a given metal; moreover the ratios determined with the (O.V.F.F nr 
repulsion) force-constants decrease in addition in the direction Sn to Si for a given 
halogen. No doubt there are considerable changes in the bond-type in the two 
series of molecules, but it seems significant that both these trends are in the direction 
expected if there is repulsion in the tetrahalides, varying with the sizes of M and of 
X—repulsion which weakens the bond in the tetrahalides relative to that in the 
diatomic halides. The figures for the carbon-chlorine bond follow the trend 
observed with the (O.V.F.F. +repulsion) figures in the series Sn to Si, except 


Table 4 
. R(MX,) 
Values of k Ratios ai 
Ge Sn Si Ge Sn 
==] 4-89 3-92 3°31 
—F, (a) 7°16 Deo 4-83* (i) 1:46 —_ — 
—F, (11) Sei a — = Gny teil — — 
=(Cil 2:63 2°33 1:96 
—Cl, (i) 3°75 3-27 2°80 (i) 1°43 1-40 1:43, 
— Cl, (ii) 2°58 2:44 2:30 (ii) 0-98 1:05 1-17 
= Br; 2:21 1:97 1-73 
— Br, (i) 2:92 2:58 2:28 (i) 1-32 1-31 1:32 
— Br, (ii) 2:03 1:95 1°85 Gi) 0-92 1:01 1-07 


* Estimated, with k(i)(MX,4)/k(MX)=1-46; this gives vy} ~715 cm for GeF,, 655 for 
SnF,. 
+ Heath and Linnett (1948 a). 


that here k, for CCl, is much smaller—as calculated in either way—than would 
be expected by comparison with CCl (Horie 1939). The ratios k,(CC1,)/A(CCl), 
are (i) 1-11 according to 8.V.F-.F., (ii) 0-47 according to (O.V.F.F. + repulsion). 
The bond in CCl, appears detiore to be abnormally weak—a conclusion which 
is also reached by the force-constant analysis of Heath and Linnett (1948 a). 
It is unfortunate that similar comparison cannot, for lack of data, at present be 
made for the other carbon halides. 

A further attack on this problem might be made by comparing the energy 
changes for the processes MX,, > M,+4X,, (£,), with those for MXN oe 
(Z£,). Unfortunately the es at present available are not Siaicendy feliabie te 
detect any general trends such as are revealed by the force-constants. Never- 
theless, the average value for the ratio E,/E, for Si-F, Si-Cl, Si-Br, Ge-Cl, 
Sn-Cl and Sn-Br is 4, Although the interpretation of E,/E, is not simple 
(Skinner 1949), it ees again that the bonds in these tetrahalides are somewhat, 
but not much, stronger than those in the monohalides. This result may therefore 


be regarded as some support for the Heath—Linnette force-constants given in 
Table 4. 
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On the Interactions of Dislocations and Solute Atoms 
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ABSTRACT. An extension of the reciprocal theorem of Colonetti is given and with it a 
general expression derived for the elastic energy of a multiply-connected solid containing 
sources of internal stress and subjected to external forces. ‘The terms in the energy which 
depend on the relative position of a selected source of stress with respect to the field of all 
the others can be separated. The expression is used to estimate approximately the inter- 
action energy between parallel edge dislocations with arbitrary slip vectors, between parallel 
screw dislocations end between an edge dislocation and a solute atom causing a spherical 
distortion of the lattice. 


Sdn IN TRO DUCEION 


ISLOCATIONS and solute atoms are sources of internal stress in a crystal 
and interact because of their associated stress fields. The interaction 
energy of a given source in a stress field is the change in energy of the 

system, including both work done by external forces and changes in the elastic 
energy of the crystal, which accompanies the introduction of the source. When 
the interaction energy depends on the position of the source, forces of strain 
interaction arise. Estimates of the interaction energy may be obtained by 
treating the crystal as a continuous elastic solid and using the linear elasticity 
theory. This theory does not permit the estimation of the energy near to the 
dislocation lines, since the displacements become singular there. ‘These lines 
are therefore excluded by small cylinders drawn about them, and we assume with 
Koehler (1941) that the energy of the material inside these cylinders is unaffected 
by the presence of other sources of stress, so that the interaction energy arising 
from changes in the elastic energy of the solid is determined by elastic energy 
changes outside these regions. We also assume that the nature of each source 
of stress is unaffected by the presence of the others; that is, each source will 
produce in the surrounding medium certain displacements and stresses, and these 
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can be combined everywhere by the Principle of Superposition. The interaction 
calculated in this way will be a sufficient approximation if the sources are not too 
close together. 

Interaction energies have been deduced by Cottrell by imagining that one 
source is ‘expanded’ in the field of another. Thus, for a solute atom causing a 
spherical distortion of the lattice in a stress field o,, (the boundary conditions 
being unspecified) the interaction energy W is given (Cottrell 1948) as 

W = = $rep*(o41 + Ogg Seq), ee (1) 
where p is the radius of the solvent atom and p(1 +e) that of the solute. This is 
obtained by assuming that W is equal to the work done by the local forces due to 
the stress field acting over a spherical surface of radius p surrounding the solute 
atom when the radius of this surface becomes p(1 +e). ‘The stress field is assumed 
constant over the spherical surface, and also constant during the ‘expansion’. 
When g;; is due to a positive edge dislocation of strength Ay situated at the origin 
in an infinite medium, and (¢,) are the polar coordinates of the solute atom, 


(1) becomes 
1+/)\ sin 
W =$p€p*Ay (; +) _ 5. oe Le, ieee (2) 
where p is the shear modulus and f is Poisson’s ratio. For a positive edge 
dislocation at the origin and a parallel negative one at (t, 4), both of strength Ao, 
Cottrell (1949) finds for the interaction energy 


A 2 
W= FTA (in —Intseos?s). | sa) ae (3) 

where L is the radius of a large cylinder containing the dislocations. "The second 
dislocation is imagined to be made by making a cut perpendicular to the slip 
vectors (that is, in the x,3 plane) and introducing an extra half-plane of atoms. 
W is assumed to be the work done by the forces of the first dislocation field acting 
on the surfaces of the cut when the two sides suffer the required relative displace 
ment. ‘The expression (3) leads to interaction forces: 


pee Oe Soe! 
1 Ob a re ee ee r 
oe VOW A) SiG SIZ ” 
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The F, force is identical with that obtained by Koehler (1941), but the function 
of (z,) in the F’, force is multiplied by a different constant. The result (4) has 
been confirmed by Eshelby (private communication). In an unpublished note 
Cottrell has shown that a result differing from (3) only by an additive constant 
is obtained by taking a section parallel to the x, x, plane and sliding the two sides 
of the cut over each other in the slip direction. 

In these estimates stresses for an edge dislocation in an infinite medium 
(L— ) with tractions applied to the inner boundary were used (see Koehler 1941). 
It was considered that other contributions to the interaction energy were of a 
higher order. We now derive a general expression for the elastic energy of a 
multiply-connected solid containing scurces of internal stress under definite 
beundary conditions on the stress fields. | Under certain conditions integrals of 


the type considered by Cottrell arise as the leading interaction terms in this 
expression, 
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§2. GENERAL EXPRESSION FOR THE ELASTIC ENERGY 


Suppose that an elastic solid B occupies a cyclic space I’, N times connected, 
which is bounded by a surface S. Let A”, r=1,2,3,. .(N—1), be a set of 
HN — 1) surfaces lying in I’ such that, if barriers are arava along them, the space 
I is reduced to a singly-connected one. If a cut is made in such a solid along 
one of the surfaces A, the two surfaces of the cut given a small relative displacement 
-as rigid bodies, and the two sides then rejoined, elastic material being added or 
removed if necessary in this operation, the solid contains a ‘ dislocation’ (Volterra 
‘ distortion’) across A” and is ina state of strain even if no external forces are acting. 
It has been shown by Volterra (1907) that if the body is in a state of strain 


Ou; Ou 
ceil a) (neal be ats 


-due to given dislocations across the surfaces A, and to the action of external 
body forces f;(«,, x2, x3) per unit volume and of external surface forces, the elastic 
-energy E may be written 


(d (N—1) 
2F = | | | fal + | I. tOudS+ > i £ fu —u\dA, 
F A(t) 


r=1 


Here u(x1, %2, x3) and t/(x,, x9, #3) =0,;v; are the components along the three 
-coordinate axes (¢=1, 2,3) of the displacement, and the traction across a surface 
whose normal is v respectively, and o;,(x1, x, 3) is the stress tensor. In order 
‘to avoid confusion with tensor suffixes, letters distinguishing a particular quantity 
of a similar set are placed as superscripts in brackets; letters in brackets after the 
‘tensor suffixes denote the usual functional dependence. The letters « and 
denote the two sides of each of the surfaces A” and the positive normal v is 
-assumed to be directed outwards from the solid at each point of S, and hence on 
the surfaces A” from the side « to the side P. 

Extensions to such a multiply-connected solid of the Reciprocal Theorem of 
Betti have been given by Volterra (1907) and by Colonetti (1915). It is not 
difficult to prove a general reciprocal formula which includes these extensions as 

-special cases. This result is required in the deduction of the expression for the 
-energy and does not seem to have been stated explicitly before. 

Suppose that two states of strain e,“) and ¢;,,°) are separately produced in the 
body B by given systems of dislocations across the surfaces A®, and by the action 
of given body forces and surface tractions, f,?, f,° and t,°, 4°, respectively. 
Let u{ and u;” be the corresponding dies Then the expression 

=1) 


Me fu; qT + UL te Dy, ods +. = at ee es ae dA” 
Aw 


ae 


may be transformed by the divergence theorem into 


(2) 
{If oe aie ine ay eal, oe aaa 


where W(1) is the strain energy function expressed in terms of e,(, and the suffix 
contractions 1=11, 2=22, 3=33, 4=23, 5=31, 6=12 are made. This last 
expression is, by a general property of ledbats functions, symmetrical in the 
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two sets of strain components e,”, e,@. Thus, denoting the expression (6),. 
regarded as a function of its numerical superscripts, by F(1,2), we have 


FU, =FO 1), | ee (7) 


which is the formula required. If f%=f=0 in I and tV=7,%%=0 on 8, 
this reduces to Volterra’s equation (2), (Volterra 1907, p. 432). If f= 0 in 
l, ¢,°°2=0 on S, u, is single-valued, and u;®) is discontinuous across only one 
of the surfaces A”, the theorem of Colonetti is obtained. 

We now suppose that the multiple-connectivity of the solid B arises from the 
removal of small cylinders of material about certain lines (which may extend across 
the whole solid or be closed in themselves), and also that there are within the solid 
singly-connected spaces not occupied by elastic material, bounded by closed 
surfaces G. The cylindrical surfaces C and the surfaces G together with the 
external surface H form the outer boundary S of the solid. In such a solid 
elastic displacements can be considered which are multi-valued about the lines — 
enclosed by the cylinders, and singular on these lines, and also singular at points 
within the surfaces G. A given source of stress (p) is then specified by defining 
the tractions or displacements on one of the surfaces G, or the discontinuity in 
displacement across one of the surfaces A” (or the tractions on H if it is an external. 
source of stress). The displacement uw; defining the state of the solid when this. 
source acts depends, of course, on the boundary conditions on all the other 
surfaces, and on any body forces f;,”) which are present. The energy E“) when 
the source (p) acts is given by equation (5) with u, replacing u,, f,” replacing 
jf, and t,°-”) replacing ¢,® respectively. If the displacement functions wu, are 
known for a number of sources separately, the energy & when they all act 
simultaneously is given by a similar expression in which the body forces, dis- 
placements and surface tractions are everywhere the sum of the corresponding 
quantities for each source separately. If we now select one of the sources defined 
by uw, and denote the sum of the displacements of all the others by u,;“, and the 
corresponding body forces and surface tractions similarly, the energy & may be 
written down by means of (5) and is easily transformed using the reciprocal: 
relation (7) into the two equivalent expressions: 


6 =EOL EM 4 | { { _fiOusad + | [ 4 eutas 
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Here E and E denote the energy when the displacements are 1 and u,(8” 
respectively. The energy when all the sources act together is thus expressed 
in terms of the elastic fields associated with the separate action of the sources of 
stress (a) and (A). In using the formulae to calculate the energy in any problem, 
the arbitrary boundary conditions for all the sources must be suitably chosen 
so that, on superposition, the final boundary conditions are those required. In 
the examples usually considered the body forces are relatively unimportant and. 
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it is assumed in the subsequent discussion that they are zero. The energy is then 
expressed as the sum of a number of surface integrals, and it is possible to examine 
the relative importance of these terms, and to separate those whose value depends 
on the relative position of the sources (a) and (R). The interaction energy can 
thus be found without evaluating a volume integral of the strain energy directly. 

When the boundary conditions are chosen in a special way a simple inter- 
pretation can be given to the various terms. If uw; is the displacement produced 
by the source (a) alone in an infinite medium, it does not depend on the position 
of (a) with respect to the surfaces G, C, although the maintenance of such a 
displacement when these surfaces are present requires the application of tractions 
to them. Thus £™ is a true ‘self-energy’ term. Under similar conditions on 
u;, E® is also a self-energy term, while the remaining terms are interaction 
terms. In general, however, all the terms give rise to interactions because u,, 
for example, depends on the relative position of the surfaces G, C and H with 
respect to the source (a). Thus in a finite solid sources will'in any case interact 
indirectly as they change their position with respect to the outer boundary. 

It may be noted that the formulae (8) and (9) may be obtained by imagining 
that cuts are made initially along all the surfaces A® so that the solid B becomes 
singly-connected. Then to produce the displacements uw” additional equal and 
opposite tractions have to be applied to the two faces of each cut to keep them in 
juxtaposition (or ina dislocated position if u;” is due to a dislocation across one 
of them). <A direct appeal to the reciprocal theorem for singly-connected solids 
can then be made, and the expressions for the energy reduce to those given, 
since the extra integrals vanish over the surfaces to which the additional tractions 
are applied. 


§3. THE INTERACTION ENERGY OF TWO PARALLEL EDGE 
DISVOCALTIONS Wilh ARBITRARY SiEIP VECTORS 
Dislocation (1) is taken along the axis of a cylindrical crystal of radius L, and 
dislocation (2) at a position x, =X, =tcos %, x, =X, =tsin x, where (x1, x)= (R, ) 
are coordinates with the centre of the cylinder as origin (Figure 1). Outside 
small cylinders C,, C, of radius 7, drawn about the dislocation lines, the crystal 
is treated as a continuous elastic solid which becomes singly-connected if cuts 
are made along A® and A®, and the interaction energy is calculated in the limit 
when Loo. The total elastic energy per unit length along the dislocation lines 
when the two dislocations are present is, by (8), 


€=ED+EO+ | [ £0 Vu, Pd + | ie. t, V[u,@® — ued A®@, 


(H+C,+ C2) 


where the superscripts refer to the fields when the two dislocations are present 
separately, and H is the outer cylindrical surface. It is assumed that tractions 
are applied over the plane ends of the cylinder in order to maintain plane strain 
(ux =ux®=0) and the integrals over these surfaces are therefore zero. The 
integral over the surface A® is zero because u,) is not discontinuous there. 
The boundary conditions on the fields of stress o,;‘) and o,; are chosen as 
follows. We assume that neither dislocation produces a discontinuity in the 
field of the other, and hence that tractions are required on C, to maintain w; and 
on C, to maintain u,. This assumption is made because, although Hooke’s 
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law is no longer valid inside C, and C,, there are in reality no holes in the lattice 
there. We assume that ¢;”=0 on C, and H, and that 7,” =0 on C, and H. 
Taking the material to be isotropic, the state of stress (1) is then defined by a 
stress function (Timpe 1905): 


R Gere af ; 
—(A,l, sin®), 
2( LE 763) © 2 Ear Arh ) 


O=7{ Rin R+ 


where the slip vector of dislocation (1) has components (A,4,,0,0), and where 
y=m(A+p)/(A+2u)z, A and pw being Lamé’s elastic constants. We denote the 
three terms in yx) by P, H and C, the last two being introduced to make ¢;”? =0 
on H and C, respectively. The slip vector of dislocation (2) is assumed to have 
components (Ag, Amz, 0). 


Figure 1. 


When Loco, E™ and E®) become self-energy terms, and with these assump- - 
‘tions the interaction energy arises from the A® and C, integrals. Taking the 
surface A® at an angle Q to Ox, and neglecting terms of the form 7,/t, the 
‘contribution to the integral from the P term in y is, when L—>oo, 


The H term gives 
W.=A,Agl y[m{ —$ —sin? O}+me{$sin2Q}],  ...... (13) 


while the C term is at most of the form 7,2/t?. 

To evaluate the C, integral we require the values of u,® on C,. We assume 
that when Loo, the value of the C, integral calculated from these values does 
not differ significantly from that calculated using values of u,;2) corresponding to 
a stress function of the same form as x, but with (A,J,sin®) replaced by 
A,(m, singd—m,cos¢), and R replaced by 7, where (7,¢) are polar coordinates 
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with (Xj, X,) as origin. Thus w;®) on C, is of the form {In1, +f(¢)} when Loo. 
On C, the tractions 7," are of the form 1/t, so that the C, integral is at most of 
the form 7/t. 


Thus, neglecting terms of order ry/t, the interaction energy is 
W=,Ahy[m,{ const. —Int—sin®?s}+m,{Fsin 2p], 2.2... (14) 


_ which is independent of Q. When m,=0, /,=+1, m,=—1, A,=A, =A, this 
expression leads to the same interaction forces as (3), which is obtained from (12) 
by making the additional substitution Q=7/2. The forces are 


jo AyAghmy 
t 
F, = “A247 tm, sin 24h —m, c08 2] hy: (15) 


Thus when m,=0, and the slip vectors are perpendicular, the leading term in 
the radial force is at most of the form (7)/t)(1/t). 

If we use only the P term in x“, the C, and H integrals in (10) no longer 
vanish but they contribute to the interaction energy terms of higher order than 
the A@ integral (12), which thus gives, apart from additive finite constants, the 
interaction energy. ‘This now depends on Q, but gives the correct forces for any 
constant {2 except Q=,7. ‘To obtain the energy in this case, it is necessary to 
include the H term in x which gives (13).* The energy of the system depends 
on where the surface A®) is taken when only the P term is used, because equilibrium 
can then be maintained only when external forces are acting. When such forces 
are present, the integrals of the tractions across any two surfaces A® and A” 
which are equivalent (that is, which can be transformed into each other by con- 
tinuous variation without passing out of the solid) are not necessarily equal, since 
these tractions are no longer the only forces keeping the elastic material between 
the surfaces in equilibrium. 

In this example the leading term in the interaction energy is the A® integral 
in (10). If we put 


x, =qg,cos O—g,sinQ; Xy = sin Q +4, cos Q; 
Q, =t cos (f—Q); Q, =tsin (4 —Q), 


so that (q,, g2) are coordinates parallel and perpendicular to A®, then the leading 
term in the force in the Q, direction (Q, constant) when the dislocation (2) is at 
(O,, Qs) can be written down by inserting the values g,=Q,, q2=Q, in the A® 
integrand. This gives for the force in the direction of the surface A® 


FENONO\uer aa ee 9 a ae (16) 


By choosing two surfaces at right angles the force in two perpendicular directions, 
and hence the force in any direction, can be found. For this rule to be valid no 
contributions of comparable order must be made to the energy by integrals other 
than a ‘cut’ integral, that is, an integral over a plane surface such that a cut along 
it would serve to introduce the dislocation. Further, if the force is to be indepen- 
dent of the surface, the energy must not depend on where this surface is taken. 


* The author is indebted to Mr. F. R. N. Nabarro for this remark. 
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§4. THE INTERACTION ENERGY OF TWO PARALLEL SCREW 
DISLOCATIONS 


We assume as before that neither dislocation produces a discontinuity in the 
field of the other. Then, taking the dislocation lines as in Figure 1, we assume 
the following stresses and displacements (Burgers 1939): 

us) =(A,/277)® ; us = (Ap/27)(p — Q); 

on) =(udy/2n)(cos®/R); 99,6! = —(pAy/2n)(sin ®/R). 

vigte (17) 
The other stresses and displacements are zero. Tangential tractions are required 
on the plane ends of the cylinders, H, and Hg, to maintain wu” but 7? =0 on H 
and C,. Tractions are required on C,, H, H, and H, to maintain wu. We again 
assume that the interaction energy calculated using this u,? does not differ 
significantly when L—>oo from that found when a displacement making t{"? =0 
on H and C, is used. The energy & stored outside the cylinders C, and C, is 
now given by (10) where the tractions and displacements are those obtained from 
(17) and the S integral extends over H, +H, +C,. On Hy, t? =0,3%, and thus 
this integral and the corresponding one over H3, are both zero. On C, the 
stresses o;;“) are of the form 1/t, so that the C, integral is at most of the form 79/t. 
Neglecting terms of the form 7,/t, therefore, the interaction energy reduces to a 
‘cut’ integral, and the leading term of the interaction force in the ‘cut’ direction 
may be deduced from (16). Taking A® at an angle Q to Ox, the force in the 
Q, direction is 
HA,dy cos (if — Q) 


7 ee 5: (18) 


giving a force in the x, direction proportional to cos(%#) when Q=0 and an x, 
‘force proportional to sin(s) when Q=7/2, and hence a total force 


ad du 2nt ios | 2 (19) 


in the % direction. Evaluating the A® integral and taking the limit as L— oo, 
we find, neglecting terms of the form 7)/t, that its value is 


(wA,Az/27)(In L —1n t). 
‘This is independent of Q and leads to the force (19). 


$5. INTERACTION OF AN EDGE DISLOCATION AND A SOLUTE ATOM 
It is assumed that the dislocation lies along the axis of a cylindrical crystal 
(Figure 2); the plane ends H, and H, are x,= +M,, andx,=—WN,. The centre 
-of the solute atom is taken at 
x, =Rcos ® =tcos¢, 
%,=Rsin® =tsing, 
x3 = 0! 


‘Outside a small cylinder C of radius 7) about the axis the continuum approximation 
is made and the crystal treated as an elastic solid which becomes singly-connected 
if a cut is made along a surface A. The stress function for the dislocation alone 
is taken to be that of (11) with A,j/;=A). Thus no discontinuity is assumed to 
arise in the dislocation field due to the presence of the solute atom, that is, the 
elastic constants of the latter are assumed to be the same as those of the surrounding 
iattice. ‘Tractions on H,, H, and G, the surface of a spherical hole of radius p 
about (¢, #), are required to maintain this state of stress. The model of the solute 


On the Interactions of Dislocations and Solute Atoms 199 


* atom, assumed spherical, is of a larger sphere of (similar) elastic material forced 
into the hole of radius p about (¢,%). We assume for the field of the solute atom 
alone ; 


‘ Ji 
OES (% ) cig FE ae aes (20) 


where y=rsinOcos¢?, y.=rsinOsingd, y,=rcos@  ...... (21) 


are coordinates with (t,%) as origin. No discontinuity in the solute atom field 
thus occurs on C. ‘Tractions on G, C, H, H,, Hy, are required to maintain this 
displacement, and we again assume that when L, M,, N,—> o the interaction 
energy calculated using uw; does not differ significantly from that obtained using 
a displacement making 7;”?)=0 on H, H, and H;. We show later that the energy 
integrals arising from the presence of tractions on these outer surfaces tend to 
zero as L, M,, Ni >. ‘The constant 6 =«p3, where « is the parameter defining 
the final radius of the solute atom in an infinite medium (Mott and Nabarro 1940). 


Figure 2. 


For the energy of the material outside G and C when both the dislocation and 
solute atom are present, (8) and (9) give 


6 =EO+EO + i | 4 Pu2dG + | i Peds (22) 
G 


(H, + H2) 
=EOLEO+ [| 40? uMdS + [J ae taer—uemda, (23) 
S) A 


where in (23) S=C+G+H+H,+H,. It can be seen immediately from (22) 
that the interaction term is given by the integral considered by Cottrell, provided 
that the second integral makes a contribution of higher order. On H, the stresses 
o,; are of the form 1/R and the displacements u, of the form 1/M,?. Thus the 
H, integral (and similarly the H, integral) tends to Zer0 as L, M,, N,> ©. On G 
we expand the stresses o;,() by the binomial theorem in powers of p/t, the expansions 
converging uniformly in @ and 4 when p<t(/2— 1). Integrating and neglecting 
terms of the form (p/t)(1/t), the expression in equation (1) results, the stresses in 
this formula being o;; (¢, #). . 
Using the reciprocal expression (23) to calculate the energy, we obtain the 
same result. The important contributions come from the integrals over G and 
A, the others being either zero or of higher order. After partial integration of 
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the tan“! term, the G integrand may be expanded inpowers of p/t. When L— o,. 
and terms of the form (p/t)(1/¢) are neglected, this integral is 
16 du?A, sin’ 
A Brae oats 
Neglecting similar terms, and taking the surface A at an angle Q to the x ,%s. 
plane, the A integral is 


45 uAg = eye (25) 


which is independent of Q. In the C ale we perform the x, integration, 
take the limit, and expand the ® integrand; thus the integral is at most of the form 
(r,/t)(1/t). The integrals over H, H, and H,, which arise because we use a 
nee which does not make ¢;” 2) =( on these surfaces, all tend to zero as. 
L, M,, N,— ©; for the stresses o4@) are of the form 1/r? and the displacement 
u;)) at most of the form In L. 

Thus the energy &, neglecting terms of the form (p/t)(1/2) and (7,/t)(1/t), is. 
the sum of (24) and (25) which is equal to the expression in equation (2). The 
formula (2) thus gives that part of the elastic energy of the material outside the 
surfaces G and C which depends on the position of the solute atom. The 
displacement of the surface of the original spherical hole when the sphere of 
similar material is inserted is u,® when the dislocation field is negligible, and in 
this condition a certain elastic energy is stored within G. When the dislocation 
field is no longer negligible, the displacement of this surface is u{?+u,?, and 
u;,”) is a function of the position of the atom in the dislocation field. Since the 
final tractions are t,”- +7,°-®) the change in the elastic energy of the material 
inside G due to the presence of the dislocation field is 


{ { 1,0» u,OdG, 
G 


where on account of the theorem of virtual work we have neglected the work 
done by the tractions ¢,”-). (This work is easily shown to be of higher order in 
this example.) This expression is exactly the G integral in (23) with the sign 
reversed, since the normal now points outwards from the centre of the sphere, 
and is thus equal to minus the expression (24). Therefore, when changes in 
the energy of the material inside G are included, (25) gives the interaction energy. 
This is greater by a factor [3(1 —f)/(1 +f)] =3/2 when f= 1/3, than the expression. 
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ABSTRACT. The magnetization of a paramagnetic salt is a relaxation process with a 
relaxation time depending on the external field and temperature. Measurements have been 
made with the object of determining relaxation times at various fields and temperatures 
in the liquid helium range. A description of the apparatus is given. 

The results show that paramagnetic relaxation is a complex phenomenon, and suggest 
that in any one specimen it must be described in terms of a distribution of relaxation times 
rather than a single relaxation time. In general we find that relaxation times are of the 
order of 10~*sec., with a temperature dependence according to a power of 1/T between 
2 and 3, while increasing with the external field H according to a power of H between 0 and 1. 


§1. INTRODUCTION 


N the statistical theory of paramagnetism the state of magnetization of a 

paramagnetic substance under given conditions of temperature and magnetic 

field is calculated by considering the statistical equilibrium of an assembly 
of magnetic dipoles. Since the classical work of Langevin (1905), the theory 
has been elaborated to take into account the quantization of the magnetic energy 
levels and the interactions of the dipoles with one another and with their surround- 
anges (Van Vleck 1932, 1937, Hebb-and Purcell 1937.) It is essentially an — 
equilibrium theory; to each set of external conditions there corresponds an 
equilibrium distribution of the dipoles, and an equilibrium magnetic moment of 
the assembly. If now the magnetic field is suddenly changed, a new equilibrium 
must be set up. ‘This requires changes in the energies of the individual dipoles, 
which can occur only if some mechanism for energy transfer exists. In general 
the new equilibrium will not be attained instantaneously; the rate of approach 
to it can be described in terms of one or more relaxation times, introduced by 
equations of the form 

dM M,—M 


dt T 


Here M is the magnetic moment of the specimen along a given direction, My is 
the equilibrium value of M, and 7 is a relaxation time. 

Interest in these time effects first arose from their analogy with dielectric 
relaxation. Waller (1932) was the first to make a study of specifically para- 
magnetic relaxation. Subsequently Kronig (1939) and Van Vleck (1940) 
introduced further ideas which have given a general qualitative understanding 
of the subject. The pioneer of experimental research in this field was Gorter 
(1936), whose first experiments were carried out at Haarlem in 1935 and 1936. 
The introduction of the adiabatic demagnetization method of attaining very low 
temperatures (Giauque and MacDougall 1933) stimulated interest in relaxation 
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effects at very low temperatures, since demagnetization experiments would not 
be fruitful if relaxation times were very long, which was the prediction of Heitler 
and Teller (1936) from an elaboration of Waller’s theory. 


§2. THE MECHANISM OF PARAMAGNETIC RELAXATION 


In paramagnetic crystals the magnetic dipoles to be considered are the electron 
spins of the magnetic ions, the interaction of the electric field of the crystalline 
lattice with the orbital magnetic moments being so strong that no changes of 
orbital moment can occur at ordinary temperatures, so that the orbital magnetism 
is ‘quenched’. Paramagnetic relaxation in these crystals includes two separate 
phenomena: spin relaxation, which is concerned with the attainment of equili- 
brium within the spin system, and spin-lattice relaxation, which is concerned 
with the interchange of energy between the spin system and the crystal lattice 
of the salt. Each spin is subject to a fluctuating magnetic field due to its 
neighbours; interchange of energy among the spins may be regarded as arising 
from absorption or induced emission in this field. The resulting spin relaxation 
time is of the order of 10-® sec. independent of temperature. The spin-lattice 
relaxation time, on the other hand, is highly temperature dependent; it is of the 
order of 10-7 sec. at room temperature and 10-* sec. at 1°. for most paramagnetic 
salts. In the first suggested mechanism for spin-lattice relaxation (Waller 1932), 
the changes in the electric field of the salt’s crystalline lattice, due to its distortion 
by thermal vibration, were considered to act directly on the interaction field of the 
spins. ‘This gave relaxation times of the order of 10"/H? sec. at liquid hydrogen 
temperatures—far longer than are found experimentally. Kronig (1939) 
however showed that the indirect action of the crystalline electric field on the 
spins through the spin-orbit coupling is much stronger than its direct action, 
and derived theoretically relaxation times of the correct order of magnitude. 

Transfer of energy between the spins and the lattice can take place by two 
processes—a ‘first order process’ in which a whole quantum of the lattice waves 
is absorbed or emitted, and a ‘ Raman process’, analogous to the Raman scattering 
of light waves, in which inelastic scattering of lattice waves takes place with 
exchange of energy between the spin system and the lattice. According to 
Fierz (1938) and Kronig (1939) the number of Raman processes is proportional 
to T* for TS 6 and to T’ for T<6, where @ is the Debye characteristic temperature 
of the lattice. The number of first order processes is proportional to T. Kronig 
(1939) found that for caesium titanium alum (for T’<@) the first order processes 
gave rise to a transition probability B= 10-!°H*T sec, and the Raman processes 
a probability D=10-°H?T"sec. The relaxation times so predicted are extremely 
short under the experimental conditions usually employed; the experimental 
values are however considerably shorter. Caesium titanium alum is a less typical 
case than potassium chrome alum, on which calculations have been made by 
Van Vleck (1940). Here the calculated relaxation times are of the correct order 
of magnitude both at liquid air temperatures, where the Raman processes are the 
more important, and at liquid helium temperatures, where the direct effects 
should predominate; but neither the temperature dependence nor the field 
dependence is correctly predicted. The relaxation time at liquid helium temper- 
atures varies more rapidly with temperature than the law r= A/T where A is 
independent of temperature, and it increases with increasing magnetic field, 
instead of decreasing, as the theory predicts. 
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§3. THE SCOPE OF THE PRESENT EXPERIMENTS 


The experiments described in this paper are concerned with spin-lattice 
relaxation only, being carried out at frequencies in the audio-frequency range 
at liquid helium temperatures. At these frequencies the spin system may be 
regarded as being always in internal equilibrium. It was hoped to obtain more 
definite results on the temperature and field dependence of the spin-lattice 
relaxation times of representative paramagnetic salts, and also to determine the 
magnetic specific heats of these salts. As the experiments progressed it was 
found that the phenomena were more complex than has been supposed till now. 
Not only do different samples prepared from the same material give different 
relaxation times, but the behaviour of a single specimen cannot be described 
in terms of a single relaxation time. On the other hand, experiments on different 
specimens of the same material, having quite different relaxation behaviour, 
yield the same value for the specific heat. 


§4. THE MEASUREMENT OF RELAXATION TIMES 


The theory of spin-lattice relaxation measurements is based on the fact that 
provided the spin system is in internal equilibrium (as is certainly the case in 
audio-frequency measurements) it can be treated as a thermodynamic entity 
separate from the crystal lattice, whose temperature may differ from that of the 
lattice (Casimir and du Pré 1938). Measurements are made by applying a large 
steady magnetic field and a small parallel oscillatory field to the salt, and measuring 
the oscillatory component of the susceptibility. (The presence of the large 
field is essential, since otherwise the energy changes associated with magnetization 
are small, and can be absorbed in the spin system without spin-lattice relaxation 
effects coming into play.) The temperature of the spin system oscillates with the 
field, but if the period is sufficiently long compared with the relaxation time, 
the spin temperature will never differ materially from the lattice temperature 
(supposed constant) and the isothermal susceptibility yp will be measured. If on 
the other hand the period is too short to allow of the exchange of energy of 
magnetization between the spin system and the lattice, another value will be 
obtained, which we may call the adiabatic value ys. At intermediate periods 
relaxation effects will be observed. 

The amplitude of the temperature oscillations of the spin system and thus 
the value of yg will clearly depend on the heat capacity of the spin system. 
‘There are two important thermodynamic relations which are independent of 
any further assumption about the relaxation time: 


ee Cee are Sth 1 
ant (1) 
ol \2 
; Ca Cf a Pear? Pe ee ee 2 
and H I (5), | (2) 


Here / is the magnetization, and Cy and C, are the specific heats of the spin system 
at constant field and at constant magnetization respectively. It is clear that if 
Xr, Xs and (1/87), are determined, the specific heats Cy and C; can be found. 
This is the basis of the specific heat measurements to be described in a later paper. 
The relaxation measurements and the specific heat measurements were carried 


out in the same set of experiments. 
14-2 
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(i) The Casimir—du Pre Equations 


Casimir and du Pré (1938) have considered the case when the period of the 
oscillatory field is comparable with the relaxation time. Their treatment assumes 
that the relaxation may be ascribed to a single relaxation time 7. If the magnetic 
field is represented by H +he'®', the susceptibility may be expressed as x’ —7x". 
Then 


x =x0(4 +(1-A) 55) vale (3) 


; WT 
and Shi (a -4) Tome a) ee (4) 


where A =C;/Cy = X5/Xr- | | 
As Casimir, Bijl and du Pré point out (1941), x’, x”, xg and x7 satisfy the relation 


(Xen Xs) FX = aa Xs) 7 eee (5) 
which represents a circle in the x’, y” plane with its centre at (4(x~+xs), 0} and 
radius 3(yp—xg). From equations (3) and (4) another relation may be obtained, 
viz. 

Vp —X SWTX. Mi, ee, ee eee (6) 


Equations (5) and (6) may be used to test the validity of the Casimir—du Pré: 
assumptions. A graph of x’ against y” for different frequencies at the same field 
and temperature should give a semi-circle, while a graph of yp— x’ against wy” 
should give a straight line. If the assumption of a single relaxation time is valid, 
its value can at once be found. 


(11) The Measurement of x’ and x" 


The measurements to be made, then, are determinations of x’ and y” at 
different frequencies, which cover as much of the relaxation range as possible. 
Liquid helium temperatures are particularly suitable for the investigation of 
spin-lattice relaxation. The values of susceptibility are very much bigger than 
at liquid air or hydrogen temperatures, a consideration which is quite important 
for such magnetically dilute salts as the cupric Tutton salts. Also, the values 
of the relaxation time at helium temperatures are normally such that most of the 
relaxation range can be covered at low audio frequencies, at which it is a relatively 
simple matter to make simultaneous measurements of y’ and yy”. 

The method used was essentially that of Casimir, Bijl and du Pré (1941). 
The specimen was surrounded by a mutual inductance, the primary of which 
carried the oscillatory magnetizing current. Changes in the magnetization of 
the specimen produced changes in the induced voltage in the secondary. The 
voltage due to the in-phase component of magnetization was compensated by a 
second mutual inductance with its primary and secondary respectively in series 
with those of the coil round the specimen. The difference AM between the 
value of the compensating inductance and its zero value in the absence of the 
specimen gave x’ in arbitrary units. The voltage-due to the out-of-phase com- 
ponent of magnetization was balanced by a potentiometer device in which a 
voltage was fed into the secondary circuit from a tapping across a resistance in 
the primary circuit (Hartshorn’s method). When the difference of this resistance 
from its zero was AR, and the frequency w/27, AR/w gave y” in the same units. 
as x’. 
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(iii) The Values of xy and xg 


In general ys could be determined without difficulty from relaxation measure- 
ments in which x has nearly the adiabatic value, yy. Usually a small extrapolation 
has tobe made. This is done on the (x’, x”) diagram, some specimens of which are 
shown in Figure 3; these give typical examples of extrapolations to yg. 

As the relaxation times were too long, it was not possible to obtain yy in other 
than zero magnetic field, when, as remarked above, spin-lattice relaxation does 
not come into play. ‘The difference between y, in zero field and its value in field 
H was always small, and was calculated from the Brillouin function: 


Xt = XT 7=0 (1 ee Ee i) 


~ pes Grea! 
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Here yu is the Bohr magneton, g the Landé factor, Rk Boltzmann’s constant and j 
the effective angular momentum quantum number, which, owing to ‘quenching’ 
of the orbital moment, is usually different from the value for the free ion. 

It was desirable to obtain x» thus rather than by an effectively static measure- 
ment (e.g. by a ballistic galvanometer method) as it was inconvenient to change 
over from A.C. to D.c. measurements and, further, it was difficult to hold the field 
current for the field H sufficiently steady to take accurate ballistic galvanometer 
readings. 


where 


7 


$55 DHE APPARATUS 
The form of apparatus used is shown in Figure 1; Figure 2 shows the electrical 
circuit diagrammatically. Small crystals of the salt investigated were normally 
contained in a spheroidal container A of Perspex, 40mm. long and 10mm. 
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Figure 1. The apparatus. 


diameter internally, which stood in the glass cryostat B, suspended by two tubes 
C and D which were joined by platinum-—glass seals to the expansion chamber of 
a small Simon expansion helium liquefier. (In the measurements on the copper 
Tutton salts, in which demagnetizing effects were unimportant on account of 
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the small susceptibility of the salts, the Perspex container was omitted, and a 
larger volume of salt was used.) The cryostat and liquefier were enclosed in a 
vacuum case, of which the part surrounding the cryostat was of glass, attached to 
the metal part above by a copper-glass seal. The entire vacuum case was 
surrounded by liquid hydrogen in a Dewar vessel. Susceptibility measurements 
were made by means of a compensated mutual inductance (immersed in the liquid 
hydrogen) of the type described by Casimir, Bijl and du Pré (1941). This system 
employs two primary coils with the numbers of turns in the inverse ratio of their 
cross-sectional areas, and three secondaries of which the two end coils are wound 
in opposition to the centre coil and compensate its mutual inductance with the 
primaries. ‘The arrangement has the advantages that in the absence of a specimen 
of salt the mutual inductance between the primary and secondary systems is 
almost zero, that the external field is very small, so that eddy currents in the 
surrounding solenoid are minimized, and that it produces a highly uniform 
field at the specimen. 
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Figure 2. The electrical circuit. 


The mutual inductance measurement was made by Hartshorn’s method 
(Figure 2), the mutual inductance being compensated by an equal and opposite 
inductance provided by a system consisting of a 4-decade mutual inductance set 
(graduated in arbitrary units of approximately 0-153 H.), and a potentiometer 
device. The latter consisted of 24 steps of 0-1 ohm and a 0-1 ohm slide wire, 
shunted to give dividing factors of 10 and 100. The shunts served the double 
purpose of reducing to a negligible quantity the inductive feed from primary to 
secondary and increasing the accuracy of reading. The primary circuit could be 
fed with frequencies of 12-5, 20, 33-3, 75, 175, 500 and 950cjs. fromea 
“Birmingham Sound Reproducers” beat-frequency oscillator, the frequencies 
being checked against a standard 50c/s. The secondary circuit was connected 
to a detector consisting of a “‘ General Radio”’ cathode-ray null detector preceded 
by a two-stage battery-operated pre-amplifier. Both the null detector and 
pre-amplifier were tuned by twin-Tee feedback networks. The limit of sensitivity 
of the detector was l»v. The primary current used was generally 0-1 amp. 
which produced an oscillatory field of about 4 oersted at the specimen. Parallel 
steady fields up to 1,000 oersted were produced by an external iron-free solenoid, 
the limit being set by the heating of the coil. The solenoid was excited by a 


battery of accumulators, and the current was measured by a potentiometer and 
standard resistance. 
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$6. PROCEDURE 

Each experiment began with the measurement of the susceptibility of the 
paramagnetic salt at different temperatures. The mutual inductance of the 
measuring coil was first determined with the specimen at liquid oxygen temperature, 
then at liquid hydrogen temperature. Helium was then liquefied in the expansion 
chamber of the Simon liquefier, and sufficient helium condensed into the cryostat 
to cover the specimen. ‘l'he mutual inductance was then measured at different 
temperatures in the helium range, the helium vapour pressure being measured by 
a manometer connected to tube D, which dips into the liquid, while the temperature 
of the helium bath was reduced by evaporating the helium by a pump connected 
to tube C. By connecting successively a mercury manometer, an oil manometer, 
and a McLeod gauge to tube D, the measurement of the susceptibility of the 
specimen could be made at temperatures down to about 0:95°K. This part of 
the experiment serves to establish the exact law of variation of the susceptibility 
of the specimen with temperature, after which it could be used as a secondary 
thermometer which could at any time be checked against the vapour pressure 
thermometer. In most of the experiments a complete calibration was made only 
at a single frequency, check measurements at the other frequencies being made at 
the hydrogen boiling point, the helium A-point, and at the temperatures at which 
relaxation measurements were made. 

With the cryostat maintained at a constant temperature by pumping helium 
through tube C at a suitable rate, mutual inductance measurements were made 
at the different frequencies in zero external field and in two or three different 
values of steady field. For the specific heat measurements, which depend on the 
ratio (x~—Xs)/xXs, fields were selected whenever possible to make this ratio 
approximately unity, so as to give the greatest accuracy in the measurement of 
the ratio. Measurements were not always made at all the available frequencies, 
but at enough to give the value of y,, and to vive information on the spin-lattice 
relaxation time. At the higher temperatures of measurement it was sometimes 
impossible to secure adiabatic conditions with the frequencies we were able to 
employ; in this case measurements at a series of frequencies usually allowed 
a value of x, to be found by extrapolation with an accuracy better than 1%. 


87s RESU LAS 
In the course of the experiments it soon became apparent that the observed 
phenomena could not be explained in terms of a single relaxation time. Thus 
precise determination of values of relaxation times and their dependence on 
temperature and external field proved impossible, and conclusions are of a 
qualitative rather than a quantitative nature. 


(i) Test of the Casimir—du Pré Relations 


Tests of the Casimir—du Pré relations (3) and (4) are made in the ways already 
described. Plots of y” against x’ lie on a semicircle for only one salt at one 
temperature of measurement—manganese ammonium sulphate at DTK. 
The curves obtained for the other salts may be roughly described as flattened 
semicircles which tail off to yp. Some examples are given in Figure 3. Even 
in the case of manganese ammonium sulphate at 2:17, however, the fact that a 
semicircle can be drawn accurately seems to be largely accidental, since it does 
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not pass through yp, and the positions of the experimental points (x’, x”) along the 
semi-circle do not agree with formulae (3) and (4); this comes out clearly in the 
plot of (yp — x’) against wx” (Figure 4) when a straight line is not obtained. None 
of the salts investigated gave a straight line on this method of presentation. In all 
the cases the curves obtained had a steadily decreasing slope. ' 

Casimir, Bijl and du Pré (1941) obtained for the plot of yx” against x’ for 
potassium chrome alum a semi-circle which, however, did not pass through x+y. 
‘They did not consider this significant and suggested its cause to be a small amount 
of ferromagnetic impurity or a slight mechanical vibration of the coils. Neither 
of these suggestions is easily acceptable. 
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Figure 3. x’, x” relaxation diagrams. 


If relaxation phenomena were due to a single relaxation time 7, we should 
find from equations (3) and (4) r=(x_—x’)/wx”. Let us call this ratio 7,—the 
apparent relaxation time—in an actual case. In any actual case it is found that 
T, 1s not constant but depends on w. For any particular specimen, then, H and 
T are not sufficient to fix +,—another parameter is required; we have chosen 
x" (Xr — x’), which specifies the position on the x’, x” relaxation diagram. Let us 
call this ratio p for convenience. If there is only a single relaxation time 7, then 
p=l1/wr. It may be expected that a reasonably true comparison of relaxation 
times will be obtained if comparisons are made of 7, at similar positions on the 
‘circle diagram’. Table 1 yives values of 7, with corresponding values of p for 
iron ammonium alum at 2:16°K. This set of results gives a typical example 
of the behaviour of 7,. 

A typical set of results (for a single crystal of iron ammonium alum) is given 
in Table 1, and a summary of all the results obtained is given in Table 2. 7, is 
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Figure 4. (xn—y’), wx” relaxation diagrams. 


Table 1. Relaxation Data for Iron Ammonium Alum (Sample 3) 


(63) 

Hai GuiKe) a H (oersted) Xerox wy” a Seca al Om) Dp 
2:16 950 400 41-4 9-9 4-2 0-040 
570 60-2 12-8 4-7 0-036 
748 TEE 11-5 6:3 0-027 
500 400 40-5 9-6 4-2 0:076 
570 59-2 10-6 5:6 0-057 
748 Willer, 10-4 70 0-045 

O25 400 33-9 6°25 5 +4 0-32 

570 52°6 7°74 6:8 0-25 

748 67-0 7°79 8-6 0-20 

S898; 400 23°7 2°84 8:4 Oc 7 

570 38°8 4-04 9-6 0-50 

748 S27 4-59 ililes 0-42 

DS 400 14:3 0-96 15 0°85 

570 24-9 1:56 16 0-80 
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given for values of p as close to 0-1 as possible to make the values of 7, comparable. 
All the samples were of Analar material. In the case of copper sulphate penta- 
hydrate the relaxation time was too short for calculation in the same way as for 
the others—the ratio (y;—x’)/wx” was rather indeterminate, both numerator 


Table 2. Relaxation Times—Summary of Results 


Salt* ii (GR) H (oersted)t Ta (Secu Om) ili Pp 
Iron ammonium alum : Dil} 285-968 0:2— 0-4 0:5 
Sample 1 0-93 285-968 1:5-— 4 0-1 
Iron ammonium alum : 1-30 400-795 0:3— 0-6 0-2 
Sample 2 0:91 343-841 1-4— 3-5 0:2 
Iron ammonium alum : 2:16 400-748 4 —755 0-1 
Sample 3 0:94 400-748 44 -—80 0-1 
Iron ammonium alum : 2:16 400-748 1:5— 2°5 0-1 
Sample 4 0:96 400-748 10. -—20 0:02. 
Manganese ammonium 2°16 514-852 5 -7 0-1 
sulphate 
Chromium ammonium 2:16 684-908 4 -5 0-1 
alum 0:98 684-908 17 -—20 0-1 
Chromium caesium alum ANS 396-678 0-8 -1-0 0-1 
0:96 396-678 13 -16 0-01 
Nickel fluosilicate 0:96 400-748 0-9— 1:3 0-1 
Gadolinium sulphate 2:16 570-851 7 -8 0-1 
octohydrate 0:99 570-851 50 0-1 
Copper ammonium 2:18 286-428 10 0-1 
sulphate 0:97 286-428 40 -50 0-05: 
Copper potassium 2:16 226-339 3 -4 0-1 
sulphate 1:0 226-339 18 -22 0-1 
Copper rubidium 1-56 172-284 5 0-1 
sulphate 0-98 172-284 15 0-1 
Copper sulphate 1:0 700 ~10-2 = 


(pentahydrate) 


* Sample 1 of iron ammonium alum was a fine powder. 


Sample 2 of iron ammonium alum was of small crystals, of average size about 10 mm‘*. 
Sample 3 of iron ammonium was a single crystal. 


Sample 4 of iron ammonium was a spheroid of compressed powder. 


t+ Where, for example, H is given as 400-748 oersted and Tas 1°5-2:5 x 10-3 sec. it may be taken: 


roughly that 7, varies from 1:5 x 10 sec. at 400 oersted to 2:5 210-3 sec. at 748 oersted. 


and denominator being nearly zero. It was therefore obtained from the Casimir— 
du Pré formula for x”, equation (4), using the value of C; (0-23 at 1°.) obtained by 


Ashmead (Duyckaerts 1942). + is too short for a proper determination and it 


can only be given as being of the order of 10-5 sec. or less. 
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(1) Effect of a ‘ Mixture’ of Relaxation Times 
Let us consider the effect on 7, of various different relaxation times throughout 
the sample. Suppose in the first place that the sample consists of 1 equal 
portions whose relaxation times are 7,,7,....7 


Then from formulae (3) and (4), 


nt 


n wr? 
ig © 
y Xr—-Xx 1 p=, 14+ 27,2 
Ls ie ? 
wy” ay, wT, 


n 1 
DL) = > = S=-; 
a 
UG; ai 
1 
(2) when aa > 1 this reduces to 
pee 
ess iz” 


where 7 etc. indicate average values. 

It will be apparent from these expressions that in case (1) a value of 7, shorter 
than the average, 7, will be obtained, while in case (2) it will be larger than the 
average, t.. Whatever distribution of relaxation times is assumed this will 
still hold. 

The departures from ideal behaviour indicated on the various methods of 
presentation are all of this nature. 7, is long for the low frequency measurements, 
decreases steadily and becomes much shorter at high frequencies of measurement. 
This is most clearly seen on the (yp — x’, wx”) presentation or from the behaviour 
of r,in Table 1. The finding of widely different relaxation times (by a factor of 
20-30) for different samples of iron ammonium alum makes it appear very 
reasonable that any one sample should not have all its constituent parts with the 
same relaxation time. It is not clear how an analysis of the relaxation results 
into a distribution of relaxation times could be made when the (x’, x”) relaxation 
diagrams do not fit any simple formula. While relaxation times are no doubt 
dependent on the chemical purity of a material, the indications are that differences 
between specimens cannot be attributed solely to this. Gorter (1947) has pointed 
out that the relaxation time has been found to alter on re-crystallization, and 
Casimir, Bijl and du Pré (1941) suggest that it depends on the ‘ physical purity’ 
of the material. 

The range of measurements and the complexity of the phenomenon do not 
permit any detailed conclusions to be drawn as to the dependence of 7 on H and T. 
The actual values of H and T at which measurements were made are shown in 
Table 2. However, using r, at a value of p of 0-1 approximately, it appears that 
zt, increases with H according to a low power of H, between 0 and a value less 
than 1. The dependence on 7, also, can only be described very roughly. The 
variation of +, would fit roughly a dependence on 1/T according to a power 
between 2 and 3, being nearer the former for the cupric Tutton salts. The 
results for the four samples of iron ammonium alum do not show any clear 
dependence on the state of aggregation of the salt. 
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§8. CONCLUSIONS 

It is not possible to describe spin-lattice relaxation effects in the paramagnetic 
salts studied at liquid helium temperatures by a single relaxation time. It is 
tentatively suggested that it might be possible to describe them by a ‘ mixture’ 
of relaxation times. Dependence on H may be said roughly to be according to 
a power between 0 and 1, varying considerably from one salt to another, and on 
T according to T-”" where n is between 2 and 3. As proportionality to 1/T 
would be expected if the relaxation were the result of ‘first order’ processes only, 
this indicates that ‘Raman processes’ must still play a considerable part in 
relaxation between 1° and 2°K. 


ACKNOWLEDGMENTS 


We wish to thank Professor Lord Cherwell for extending to us the facilities 
of the Clarendon Laboratory, Professor F. E. Simon for his interest and encourage- 
ment, and Mr. L. E. Drain for many helpful suggestions. Grants were made by 
the Royal Society for the purchase of the beat-frequency oscillator and the null 
detector used in the experiments. One of us (R.J.B.) is indebted to the Ministry 
of Education for a maintenance grant. 


REFERENCES 


Casimir, H. B. G., BijL, D., and pu Pr, F. K., 1941, Physica, 8, 449. 

Casimir, H. B. G., and pu Pr, F. K., 1938, Physica, 5, 507. 

DuyckaeErTs, G., 1942, Mém. Soc. Roy. Sct., Liége, 4° série, tome VI, 195. 

FiERz, M., 1938, Physica, 5, 433. 

GrauquE, W. F., and MacDoucaLt, D. P., 1933, Phys. Rev., 43, 768. 

Gorter, C. J., 1936, Physica, 3, 503; 1947, Paramagnetic Relaxation (Amsterdam : Elsevier 
Publishing Co.), Chapter III. 

Hess, M. H., and Purce Lt, E. M., 1937, 7. Chem. Phys., 5, 338. 

HEITLER, W., and TELLER, E., 1936, Proc. Roy. Soc. A, 155, 629. 

*Kronic, R. DE L., 1939, Physica, 6, 33. 

LaNnGEVIN, P., 1905, 7. Phys. Radium, 4, 678, Ann. Chim. (Phys.), 5, 70. 

Van VLECK, J. H., 1932, Electric and Magnetic Susceptibilities (Oxford : University Press), 
Chapter XI ; 1937, 7. Chem. Phys., 5, 320; 1940, Phys. Rev., 57, 426. 

“WALLER, I., 1932, Z. Phys., 79, 370. 


21% 


Specific Heats of some Paramagnetic Salts at 
Temperatures near 1° x. 


By R, J. BENZIE anp A. H. COOKE 
Clarendon Laboratory, Oxford 


MS. received 29th August 1949 


ABSTRACT. The magnetic specific heats of gadolinium sulphate and of several salts of 
the iron group of transition elements have been measured at temperatures near 1° K. by a 
method depending on paramagnetic relaxation. The results obtained are analysed into 
contributions due to magnetic dipole coupling and interaction of the magnetic ions with the 
electric field of the crystalline lattice. Cases in which these contributions are insufficient 
to account for the observed specific heats are discussed. 


§1. INTRODUCTION 


INCE the first successful application of the magnetic method of cooling 
by Giauque and MacDougall (1933) there have been many experiments 
employing the magnetic method to reach temperatures below 1°K..,. 

using a variety of paramagnetic salts as cooling agents. As is well known, the 
temperatures attained are measured in the first approximation on the ‘ magnetic 
temperature scale’, obtained by the extrapolation of Curie’s law of paramagnetism 
known to be valid at high temperatures, to the temperatures reached by adiabatic 
demagnetization. The two most important properties of the salts below 1°k. 
are therefore the entropy, whose temperature variation determines the temperature 
reached on demagnetization, and the magnetic susceptibility, which determines 
the relation of the magnetic temperature to the absolute temperature of the salt. 
In this paper we describe specific heat measurements on a number of paramagnetic 
salts in the temperature range 1° kK. to 2° K., which provide the data for calculating 
the propeities of the salts below 1°. by the theoretical formulae yiven by Van 
Vleck (1937) and by Hebb and Purcell (1937). This information is obtained by 
the analysis of the specific heat into the contribution due to magnetic dipole 
interaction of the magnetic ions and that due to che splitting of their energy 
levels by the Stark effect of the electric field of the crystalline lattice. At the 
temperatures at which we have made measurements both effects contribute 
terms to the specific heat which vary as the inverse square of the temperature. 
The magnetic dipole interaction is calculable from the susceptibility, which can 
be measured at high temperatures. Our experiments then give the Stark splitting 
of the eneryy levels. Values of the splitting are derived for iron ammonium alum, 
manganous ammonium sulphate, gadolinium sulphate, nickel fluosilicate, two 
chromic alums and three cupric Tutton salts. It will be seen that besides the 
salts commonly used in magnetic cooling work, this list includes a number of 
others which have been recently investigated by radio-frequency spectroscopy 
(Bagguley et al. 1948). The interest here lies in the comparison of the spectro- 
scopic determinations of the energy levels (made at relatively high temperatures) 
with the values derived from the specific heat measurements. In some of the: 
cases in which this comparison is possible (chromium ammonium alum and. 
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chromium caesium alum) the values are in fair agreement, but in the copper 
Tutton salts, for which there is no Stark splitting of the spin-degenerate energy 
levels, the specific heat is much too large to be accounted for by magnetic inter- 
_action alone. The excess specific heat is attributed partly to exchange coupling 
of the magnetic ions and partly to hyperfine splitting of the nuclear energy levels. 


§2. METHOD OF MEASUREMENT 


‘The magnetic method of measuring specific heats (first described by de Haas 
-and du Pré (1939)) which we have employed, consists essentially in the 
determination of the susceptibility of the magnetic ions of the salt under isothermal 
conditions and under adiabatic conditions. From these measurements the 
specific heat of the system of magnetic ions is obtained by a simple thermodynamic 
relation. The isothermal measurement is made in a straightforward way. The 
salt, normally in the form of small crystals, is immersed in a bath of liquid helium, 
and observations of the susceptibility of the salt are made from its influence on 
the mutual inductance of two coils that surround it. Measurement of the 
adiabatic susceptibility by isolating the salt from the helium bath would be 
difficult experimentally, and would, moreover, have the disadvantage for our 
purp-se that the specific heat so determined would include the lattice specific 
heat of the salt, which in the helium temperature range is of the same order of 
magnitude as the specific heat of the magnetic ions. We can, however, measure 
the susceptibility under conditions in which no heat transfer is possible between 
the magnetic ions and the crystalline lattice of the salt, and so obtain the specific 
heat of the ions alone. To do this, use is made of the finite time of relaxation 
for thermal transfer from the magnetic ions to the surrounding heat bath formed 
by the crystal lattice and liquid helium. When a large steady field and a small 
parallel oscillatory field are supplied to the salt, the ‘temperature’ of the magnetic 
ions will oscillate with the field, rising as it increases and falling as it decreases. 
These oscillations of temperature will be reduced by heat transfer between the 
ions and the heat bath, and if the period of oscillation is sufficiently long the 
variations of magnetization will take place under practically isothermal conditions. 
If however the period of oscillation is very short compared with the time constant 
for heat transfer, the variations will be adiabatic. By measuring the influence 
of the salt on the mutual inductance of two coils surrounding it at low and at 
high frequencies it is thus in principle possible to determine the isothermal and 
adiabatic values of the magnetic susceptibility, y,=(0l/0H)» and ys=(0//0H)z, 
under identical conditions of field and temperature. In practice there are 
difficulties in doing this directly, but they are avoided as described below. The 
specific heat is related to these quantities by the thermodynamical equations 


Cy/Cy = Xr/Xs 


and 
2 
Cues (5) i 
oT H Xp 
from which we derive 
ol \2 
= eon Xs 
= 
NO ee) 


where J is the magnetic moment per cm’, 
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If the salt behaves as an ideal paramagnetic, i.e. if 7=f(H/T), then 


é ay SE lint dae 
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FA? xrxs 
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(The specific heat at constant magnetization, C\, is a more convenient quantity 
than that at constant field, Cy, since for an ideal paramagnetic C, is independent 
of the magnetic field at a given temperature. It is equal to the value of Cy at 
H=05 

In practice the assumption that the salt behaves as an ideal paramagnetic 
is not realized. Interaction and demagnetizing effects cause the magnetization 
to obey the equation J=f{H/(T—f)} rather than T=f(H/T). Moreover, it is 
normally difficult to determine the value of y» in a magnetic field due to the 
relaxation time being too long, and so y, is usually measured in zero field, where 
relaxation effects do not affect the measurement, since yp=yg at H=0. It is 
then necessary to correct for the effect of saturation in order to arrive at the value 
of yxy; in a field H. This correction is always small and has on occasion been 
checked by a ballistic galvanometer comparison of yy in field H with its value in 
zero field. If the relaxation time is too short, it may not be possible to obtain 
the value of yg from a measurement at a single frequency, but it may be derived 
without difficulty from measurements at several frequencies which nearly 
give the adiabatic condition by making an extrapolation to yg as described in 
the preceding paper on paramagnetic relaxation (Benzie and Cooke 1950). 

For a substance obeying the equation y;=C/(T—f), 

(es gee ca 
( t B)* Xr Xs 
Including the saturation correction 
Rice ate) 


and ignoring terms in A to a power higher than the first, we have 
Ci; Zs Xs 


giving C,= 


Ce Tapp gal PART) 
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2 
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where C is the Curie constant for the salt, » the Bohr magneton, g the Landé 
splitting factor, k Boltzmann’s constant, and 7 the effective angular momentum 
of the ion, which owing to ‘ quenching’ of the orbital component has a ‘ spin-only’ 
value for the salts we have investigated. R is the gas constant, 

It will be seen that since we know the value of the Curie constant, which is 
well established for all the salts we have used, the experiment is reduced to the 
determination of the ratio of yg to xp,9. It is not necessary to know f and A very 
accurately as they occur only in correction terms. 
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§3. THE APPARATUS AND PROCEDURE 
Details of the apparatus and procedure during the experiment are given in 
the preceding paper (Benzie and Cooke 1950). 


§4. RESULTS 

The main source of the specific heats which we have measured is usually the 
Stark splitting of the spin-degenerate energy levels of the magnetic ions in the 
crystalline electric field of the lattice (Hebb and Purcell 1937). ‘To derive the 
values of such splittings we must first calculate the component of specific heat 
due to magnetic dipole interaction from the formula given by Van Vleck (1937), viz. 

Cr T?/R 7 Or; 
where 7/3 is the Curie constant per unit volume and Q is a factor depending on the 
lattice arrangement of the magnetic ions. ‘Taking this to be a face-centred 
cubic arrangement for the alums, Q= 14-4 and C,,T?/R=2-47?. Assuming that 
the Tutton salts have the same crystal structure as magnesium ammonium 
sulphate (Hofmann 1931), with crystallographic axes in the ratio 3:4:2, Q is 
lo-Sgiving Chl R=2377; 

In the absence of any other contribution to the specific heat, the remainder 
after subtracting the magnetic interaction contribution can be attributed to the 
Stark splitting. It is shown below however that the values of Stark splitting 
so derived are consistently greater than those determined spectroscopically. 
The discrepancy may be explained by the presence of exchange interaction, 
which contributes to the specific heat. Evidence of the presence of exchange 
interaction is provided by the magnetic behaviour of the salts. Several of them 
are found to obey a Curie-Weiss law with a Weiss constant which is not wholly 
attributable to magnetic dipole interaction. 


Iron Ammonium Alum. 


As an example of the method of working, the results of an experiment on a 
specimen of iron ammonium alum are given in some detail. ‘This specimen was 
a spheroid of fine compressed powder, for which the value of Weiss constant f, 
calculated from the Lorentz internal field and demagnetizing effects, was 0-05°. 
This agreed with experiment without the addition of any Weiss constant due to 
exchange effects. 

Experiments on two other specimens, one a loosely packed fine powder, and 
the other consisting of small crystals of about 10 mm average size, gave values 
of C,T?/R at 0-93°K. of 0:0144 and 0-0143 respectively. In both cases the 
relaxation time was too short for values to be obtained at higher temperatures. 

The Fet+* ion is in a §S state. If we assume that the water molecules 
surrounding each ion produce an electric field of cubic symmetry, the six-fold 
degenerate level will split into one of two-fold and one of four-fold degeneracy 
(Van Vleck and Penney 1934). The partition function per ion is then 


Z=2{1+2exp(—d/RT)} 
where 6 is the energy separation of the levels. At temperatures high compared 
with 6/R this leads to a contribution to the specific heat given by 
CT?/R = (2/9)(8/k)?. 
With this arrangement of energy levels the deviation of C,T?/R from this square 
law value is +3% at 2°K. and +6% at 1°K. if 8 is taken positive and —3-5% and 
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—7% if} is taken negative. Our measurements, however, show a variation of less 
than 1°, between 2-2° k. and 0-96°k. This is more consistent with the assumption 
of the splitting of the six-fold spin degenerate level into three equally spaced 
doublets, giving a partition function 
Z=2{1 +exp(—t3/RT) +exp(—S/RT)}, 
which would result ina calculated falling off of C,T2/R by 0:5°%, from 2°. to 1°K. 
Further evidence of the near-constancy of C,T?/R down to 1°% is found in the 
agreement of the present measurement with the most recent determination—by 
Broer (Gorter 1947)—at liquid air temperature; Broer’s value of C,T2/C is 0-27. 
With the assumption of three equally spaced doublets the high temperature value 
of CT*/R due to the splitting is 3(3/k)?, where 8 is the overall splitting. This 
leads to a value of 0-23 for 8/k. 
Table 1 


IPCs) tal Xr. 9 Xs Xs r — — 


(oe.) 


2:16 400 42) 66:5 0:639  ~=0-001 2-50 ) 
570 61:5 47-5 1-295 O-0025 2°56 = 81-074 143 e270 
748 TESS Sa 2-161 0:004 2-53 J 


) 
| 
| 
isis 400 73°8 111-25 0-663" 0-003: 2:42) | | 
570. 106-6 78-4 12360: +0-006 2°40 > 171123 1-43 2-70 > 0:20 
748 128-7 56-3 2-284 0-011 2-40} | 


0:96 400 104-3 149-0 0-700 0-006 2-30 } 
570 = 148-7 104-6 1-421 0-013, 2-312 1-1:72 1422-69 
7458) 17827 74-05 2-395 G-02i 2208) J 


The values given for the y’s are actually mutual inductance measurements, but these are 
proportional to the susceptibility, and are therefore taken as susceptibilities in arbitrary 
ea es 5 me WEPIIEE 
units ; 1=(Xq 9 ~Xg)/Xg3 7,=r—AH?/T?. 


A crystalline field of symmetry less than cubic would lead to a greater removal 
of degeneracy than that found by Van Vleck and Penney, and would make the 
splitting into three doublets possible. That the crystalline field in the alums 
has a component with symmetry less than cubic is clear from the existence of 
a splitting of the spin degenerate levels of the chrome alums—there would be 
no such splitting in a purely cubic field (Hebb and Purcell 1937). 

It would be unwise to draw definite conclusions as to the arrangement of the 
split levels from our results, but it does seem certain that the split levels are three 
doublets, rather than one doublet and one quadruplet, and that the spacing is 
not far from being equal. 

Preliminary spectroscopic measurements by Bagguley e¢ al. (1948) suggest 
however that the overall splitting is less than 0-1 cm™!. If this is the case 
we must conclude that the theory we have assumed in deriving 3/R is inadequate 
but we are unable to offer any explanaticn of the discrepancy. It is hardly 
possible that exchange interaction could account for the difference. 

Table 2 summarizes the results of specific heat measurements on different 
substances. Besides the value of C,T?/R, we give the value of C,;T?/C, where 
C is the Curie constant for the substance. Gorter in his book Paramagnetic 
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Relaxation (1947) gives most specific heat results in terms of this quantity—it is 
his ““b/C”’. The value of the crystalline splitting of the energy levels, deduced 
by the same procedure as for iron alum, is also given. The results for the 
different substances are discussed in the notes below. 


Table 2 
Cote CRI on GC: 
Substance a (GRE) a a k (°) Spin = 
x 108 
Iron ammonium alum 2:16 0-0143 0-270 ) 5/2 — 
Sample C - 1-31 0:0143 0-270 > 0-20 
0-96 0-0142 0-269 
Sample A O°98 0-0144 
Sample B 0593 0-0143 
Manganous ammonium Do, 0-0331 0-63 : fe ee 
sulphate 0:95  0-0330 0-63 \ ne ue 
Ammonium chrome 2:16 0-0370 1-64 0-38 3/2 — 
alum 1-54 0-0365 1-61 
S27) 0:0361 1-60 
0:97 0:0354 eSsz 
Caesium chrome alum 2:16 0:0109 0:481 0:20 3/2 — 
1-82 0:0108 0-480 
1-60 0-0108 0-480 
1-02 0-0106 0-473 
Nickel fluosilicate 0-95 0-0139 0-89 0-23 1 — 
Gadolinium sulphate 2:165 0-305 3-24 ibo3}5) 7/2 0-065 
1:50 0-274 2-92 0-122 
1-20 0-242 2°58 0-168 
1-055 02228 2-42 0-204 
0-994 -— 0-217 2°30 0-220 
7 bes (Opal Cyl Weiss : (Ci—Ch) a 
Sukstance EXC TRS R —G constant Spin ae 
rollOme >< le SO 
Copper ammonium Als 8-701 1:62 0) 1/2 7:4 
sulphate 1-81 8:70 1-62 7-4 
1-59 8-69 1-G2ae 7-4 
1-30 8:55 1-59 7:2 
0-97 8-43 (S57 Jil 
Copper rubidium 1:56 3-36 0-61 0) 1/2 Dall 
sulphate 135 6283 0-61 2) 
1:104 3-30 0-60 2-0 
0:98 3°30 0-60 2:0 
Copper potassium 2:16 ES fleilil 0-035 12) 
sulphate 1-49 6:01 ied +0:005 s & 
1:28 5-93 au 
1:016 5-98 lea J 
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Manganous Ammonium Sulphate. 


The specimen was a loosely packed powder contained in an ellipsoidal shell. 
The value of 0-63 x 10° for b/C is in satisfactory agreement with the value of 
0-64 x 10° quoted by Gorter (1947) presumably obtained at liquid air temperature. 
The value of 5/k is obtained by assuming that the ®S ground state is split by the 
crystalline electric field into a two-fold and a four-fold level. The value obtained 
agrees with that deduced from measurements below 1°xk. (Cooke 1949), but as 
Bleaney and Ingram (to be published) have detected a nuclear hyperfine structure 
in the paramagnetic resonance spectrum of much diluted manganous salts, 
the value of 6/k may be too high, since our experiments would not distinguish 
between the electronic and nuclear splittings. In addition, there will probably 
be an appreciable contribution from exchange interaction, though it is likely 
to be small compared with the contributions from the splittings. The resolution 
of the nuclear and electronic specific heat anomalies would require experiments 
at temperatures of 0-1°k. or lower. 


Chromic Alums. 


The quenching of the orbital magnetic moment of the Cr+++ ion by the 
crystalline field converts the 4F ground state effectively into 4S. A crystalline 
field of symmetry lowerx than cubic causes splitting of this state into two doublets 
of energy separation 6 (Hebb and Purcell 1937). As this is the maximum extent 
of splitting which an electric field can produce (Kramers 1930) there can only 
be one splitting associated with each ion. ‘The specific heat associated with the 
splitting will be given by 

ee Oe Aexp (ORL) 
es =1(5) Tremp RT 
At high temperatures (6<RT) this reduces to }(6/k)?. 


Ammonium Chrome Alum. 


The specimen was of loosely packed powder. 

After allowing for the Lorentz and demagnetizing effects, and the effect of 
the Stark splitting on the susceptibility, a Weiss constant of —0-03 was required 
to express the temperature dependence of the susceptibility of this salt. 

The value of C,T?/R fell with temperature, which is to be expected in the 
temperature range of our measurements, since Cy7?/R is constant only for 
temperatures high compared with 3/k. The fall from 2° to 1° is 35%, in close 
agreement with 3° calculated theoretically. 

Paramagnetic resonance experiments at 77°K. and 20°xK. by Bleaney and 
Penrose (1948) and Bleaney (1949) showed two splittings—0-45° and 0-34". 
Since there can be only one splitting associated with each ion, we have to conclude 
that different ions in the unit cell have different splittings. ‘To reconcile our 
results with those of Bleaney and Penrose would require three ions having the 
larger splitting to one ion having the smaller. 


Caesitum Chrome Alum. 

Our experiments were made with a single crystal ground to a rough sphere. 
After the appropriate allowance had been made for the Lorentz and demagnetizing 
fields a Weiss constant of — 0-02 was found. 

Our value of 0:20° for 5/k compares with 0-19 from paramagnetic resonance 
data (Bleaney 1949). 
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| Nickel Fluosilicate. 


The specimen was in the form of a loose powder. A Weiss constant of 0-1 
was found after allowance had been made for the Lorentz and demagnetizing 
effects. 

This substance is of special interest for paramagnetic resonance work, as it 
has only one ion per unit cell. The ground state of the Ni*+** ion is ?F,, but as a 
result of quenching of the orbital levels it is effectively in an S state but with a 
Landé g-factor of 2:27, due to the influence of the higher orbital levels. Assuming 
that the electric crystalline field splits the three spin levels into a singlet 7=0 
and a doublet 7 = +1, we find the separation §=0-23°. Penrose (Bagguley et al. 
1948) found a single splitting of 0-17° at liquid hydrogen temperatures from para- 
magnetic resonance experiments. The Weiss constant of 0-1° for the salt implies 
the existence of exchange coupling, which will contribute to the specific heat 
(Opechowski 1937, 1948), accounting partly at least for the difference between the 
spectroscopic and our value of 5/R. 


Gadolinium Sulphate Octohydrate. 


The specimen, in the form of a Icose powder in an ellipsoidal container, was 
prepared from gadolinium oxide kindly supplied by Dr. J. R. Marsh, of the 
Dyson Perrins Laboratory, Oxford. 

For a substance of the high susceptibility and large Stark splitting which 
gadolinium sulphate possesses, the usual Curie-Weiss expression for the suscepti- 
bility isinadequate. It was found possible, however, to express it as C/{T — B(T)}, 
where A(T) is a semi-empirical function of 7, corresponding to a behaviour 
intermediate between what is given by the Onsager (1936) and Van Vleck (1937) 
improvements to the Lorentz internal field theory, together with a correction 
for the effect of the Stark splitting (Hebb and Purcell 1937). This necessitates 
a further correction factor (1—dB/dT)? to the expression for C,T?/R. 

Our values for the specific heat are in good agreement with those of other 
workers (Clark and Keesom 1935, de Haas and du Pré 1939, and van Dijk and 
Auer 1942) at liquid helium temperatures. 

The ground state of the Gdt+t++ ion is 8S, which is split by an electric 
field of cubic symmetry into two doubly degenerate levels and a four-fold 
degenerate level between them, the ratio of the two separations being 
3:5 (Penney and Schlapp 1932). The partition function for one ion 
Z=2{1+2exp(—35/8RT)+exp(—d/RT)} where 6 is the overall splitting, 
assuming that the four-fold level is nearer the lower of the other two, as this gives 
the better agreement with experiment (van Dijk and Auer 1942). The con- 
tribution to C,T?/R arising from the splitting 

Dy 3 d? 

Grad ea Soe A TP 

cannot be replaced by its high temperature approximation between 1° and 2° x. 

as 6 is too large, and the full calculation must therefore be made. The magnetic 

interaction specific heat has been calculated from an expression given by Hebb 

and Purcell (1937) as the approximate Van Vleck formula for C,,T?/R given 
above also becomes inadequate. 

Our value cf 1-35° for 6/R is chosen to fit the experimental value of specific 
heat at 2-16°, as the calculated magnetic specific heat should be quite reliable here. 


(T1nZ) 
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The observed specific heat, however, falls off more rapidly than that calculated 
using this splitting, C,/R being 0-22 instead of the calculated 0-25 at 1°. This 
may be either because the calculated value of the magnetic interaction specific 
heat is in error or because the assumed disposition of the split levels is wrong. 
This trend is further confirmed by the value of b/C found by de Vrijer, Volger 
and Gorter (1946)—3-9 x 108 at room temperature and liquid air temperature— 
which is equivalent to a splitting of 1-4°. 


The Copper Tutton Salts. 


The ground state of the Cu*+ ion is 2D, but the orbital degeneracy is 
completely removed by the crystalline electric field. The distance of the higher 
orbital levels is such that they have a second-order effect on the susceptibility 
in that g=2-19+0-01. Otherwise the susceptibility has the spin-only value. 

The copper salts are of particular interest because, the Cut+ ion having only 
one electron spin, there can be no splitting of the spin levels by the crystalline 
electric field (Kramers 1930). There is therefore no Stark specific heat. It 
was thought that the specific heat could therefore be analysed into contributions 
due to exchange interaction and to magnetic dipole interaction, and that, as the 
magnetic interaction contribution is calculable with fair accuracy, information 
on the exchange interaction could be obtained. In fact it transpired that the 
situation is less simple. The discovery by Penrose (1949) of a nuclear hyperfine 
structure with an overall splitting of the order of 0-02 cm™ in very dilute cupric 
salts made it clear that this splitting will contribute appreciably to the specific heat. 


The Ammonium, Potassium and Rubidium Copper Tutton Salts. 


All the specimens were in the form of loosely packed powders. After allowing 
for the Lorentz and demagnetizing effects the ammonium and rubidium salts 
obeyed Curie’s law exactly down to 1°x. A Weiss constant of 0-01 would have 
been easily detectable. In the case of the potassium salt a Weiss constant 
of 0-035 + 0-005 was found. 

The value of C,T?/R found for the potassium salt (6:0 x 10~*) disagrees with 
the value of 6-8 x 10-* found by de Klerk (1946) below 1°x., but it is in good 
agreement with 6-1 x 10-4 found by Garrett (1949), also below 1° x. 

In all three cases the specific heat considerably exceeds the value attributable 
to magnetic interaction (C,,7?/R=1-3 x 10-4 approx.), varying from about seven 
times in the case of the ammonium salt to 2} times in the case of the rubidium 
salt. It seems reasonable to suppose that the contribution from the nuclear 
splitting is approximately the same in all three salts, and therefore, in the ammonium 
and potassium salts at least, the greatest part of the specific heat must be attributed 
to exchange interaction (Opechowski 1937, 1948). 

A further series of experiments, using diluted copper Tutton salts, has been 
undertaken to elucidate the nuclear contribution to the specific heat. 
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ABSTRACT. The spectra of the photoelectrons ejected from various targets by the 
y-rays of Th(B+C-+-C”) have been measured and the emission probabilities of the y-rays 
of energy 0-238, 0-277, 0-299, 0-510, 0-582, 0:726, and 0-859 mev. have been estimated. 
Using the known intensities of the internal conversion lines, coefficients of internal 
conversion haye been found. The four low-energy rays seem to be magnetic dipole in 
character with some electric quadripole admixture, but the strong 0:582 Mev. ray of 2°8Pb is 
predominantly electric quadripole. The possible spins of the nuclear levels are discussed. 


§1. INTRODUCTION 


HE aim of the present work was to find the intensities of certain y-rays of 
| thorium (B+C+C’") in order to elucidate their multipole character and 
to check some contradictions in the literature. For example, Itoh and 
Watase (1941) found the 0:510 and 0-582 Mev. y-rays of 7°°Pb to have nearly equal 
intensity, whereas Thibaud (1926) and Johansson (1947) found their ‘external’ 
photoelectron lines to be in a ratio of about 1:3. The two latter authors have 
assigned to the G (0:277Mev.) y-ray of 2°%Pb a photoelectron line of height 
about one-third that of the line due to the strong F (0-238 Mev.) y-ray of ThB— C. 
So large an intensity would make the G y-ray electric dipole in character. 

The intensities have been calculated by using a modified form of the relation 
of Ellis and Aston (1930), which applies in the case of a thick target. ‘The relation 
can be tested at quantum energies hy of 0-238 and 2-62 Mev. because independent 
estimates of intensity for these two strong y-rays are available from f-ray data. 
The theory seems to overestimate the intensity ratio of the high to the low energy 
y-ray by about 50°%, perhaps because a difference in the angular distributions of 
the photoelectrons ejected by hard and soft y-rays may affect the probability of 
detection with our arrangement of source and target, shown in Figure 5(6). In 
favour of this view is the fact that an even larger discrepancy has been found in 
some measurements using the alternative method of very thin targets. 

The source consisted of about 60 mc. of radiothorium in thorium oxide carrier 
held in a copper tube 0-09 cm. thick and 3/16in. outside diameter with soldered 
ends. A flat square target was waxed to a flattened part of the tube. The surface 
of the target, about 1-6cm. wide, was inclined so that it made an angle of 22° 
with the central ray of the focused beam. 

The magnetic spectrometer described by Martin, Richardson and Hsii (1948) 
was used, the magnetic field strength being obtained from the voltage generated 
in a coil rotating at constant speed. The f-ray counts were recorded with a 


scale of 32. 
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The spectra of the ejected electrons are shown in Figures 1 and 2, in which the 
notation resembles that used by Ellis (1932) for the internal conversion lines. 
For example, the F y-ray of 0-238 mev. ejects the F, I and J photoelectron groups 
from the K, L and M shells of the target atoms. The background due to y-rays 
penetrating the lead shielding was about 1,800 counts/minute. To reduce 
statistical errors about 40,000 electrons were counted for each plotted point. 

40 = , oe 


hundreds) 


Counts /Minute 


Curves 6,7 


| ease Ob 7 DANEEL |e | aes es ali le a 17 
l" I f 
fo a a 090 — 309 —— shy ag ag — 
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Figure 1. The photoelectron and Compton electron spectrum of RdTh and its products from 
targets of (1) 200 mg/cm?, (2) 3-7 mg/cm?, (3) 0-72 mg/cm’, (4) 7-4 mg/cm? gold, 
(5) thorium, (6) polythene and (7) copper. The y-ray background for each curve is indicated 
by the letter b. The space between ordinate lines represents 400 counts/min. 


The height of the line due to the 0-51 Mev. (L) y-ray is about 4 that of the 
0-582 Mev. (M) ray while that of the 0-277 Mev. (G) ray is very small, about 1/10 
that of the 0-238 mev. (F)iray. © The I and Hilines unresolved with gold, are 
Just distinguishable with thorium, which also shows the J Jines “Cutvech 2a 
and 4 show sharper lines from thin layers of gold mounted on an aluminium base. 
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As the thickness of gold is increased the added photoelectrons build up the low 
energy side of the line because they come from deeper layers in the target. The 
Compton recoil electron groups are well shown in curve 6 for polythene which, 
because of its low mean atomic number, emits few photoelectrons. With copper, 
curve 7, the F photoelectron line is already appreciable. 

Curve 


12,4 67 Curve 
33 5 


25 


Coun ts/Minute (hundreds) 


eet ee ee ee Se 
8400 9000 9600 10200 10,800 11,400 
Approx. Ho 


Figure 2. The spectrum near the 2:62 Mev. y-ray lines. 


‘Pre-absorption’ of the 0-238 Mev. y-ray in the thick target of 200 mg/cm? of 
gold is shown by the fact that the F line is lower with that target than it is with one 
of 3-7 mg/cm? of the same material. 


$3; ANALYSIS OF THE LINES FROM A THICK TARGET 


The heights # of the line in counts/min. obtained after subtracting the 
underlying background are given in Table 1. ‘The X line is corrected for the 
incompletely resolved Xa line from the L shell and the M line is corrected for the 
nearly coincident line ejected by the 0-51 Mev. y-ray from the L shell, assuming 
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7, the atomic absorption ccefficient for the K shell, tobe 4 x7z._ Above 0-34 Mev., 
7, has been interpolated from the curves of Hulme, McDougall, Buckingham and 
Fowler (1935). For the 0-238 ey. y-ray the formula of L. H. Gray (1930) for 
Z = 82 has been used: 


log (1:25 x 7x) =log A + 0-48( log A)? — 24-8684, 


where A is the wavelength in x.u. The values for gold and thorium were found 
by multiplying by (Z/82) and, for the G and H rays, by interpolating between these 
values and the curves of Hulme et al. 


Table. 1. Heights h of Photoelectron Lines from Thick Targets and 
Calculated y-Ray Intensities BJ’ in Quanta per Disintegration 


Nucleus ‘hCG. 2Pb Thc. 2ePh Ph Stin@ ae Epeea eo 
y-ray e) G H L M O 1 xX 
hv (Mev.) 0-238 0-277 0:299 0-510 0-582 0-726 0:859 2-62 
Tx x 10”? Au 1225 8-8 ED 1:90. (1:39: ~~ 0:86 0:6 0 208 


(cm?/atom) \e Rhee Zit 14:8 12-2 Oil 2°38 1-49 1:06 0-184 


. . fAu=1630 160-144” 304 9948) 176 encoun 
h (counts/min)" ~* rTh 1200). 4112) 76 ee 204 mens ee ee a 


% lossinmica fAu 3-6 Zo NGS 0 — — -— — 
window Te 25.0 3-0 Dee, 0 — — — = 
-, ~fAu 0:50. :0:57 (0-61 0-78 0:82 © 0:83" "0.855 10.995 
y-tay transmisston\'Th 0:35. 0:43 «0-49. 0-73 9 -0-76n 60-811 =O 6cmmm atm 
Au 31-4 26-4 246 17:0 16:2 15:0 14:5 14-6 
Q (Th 35-6 28-7 261 17-2 16:0 ° 14:8 14:3 14:3 
I’ (relative no. Au *0°85 - 0-075. ‘1 (30-495 (0-665 O-11um O-140mmet oo 
of quanta) Th 0-94 0:068 0:10 0-17 0:64 0-11 0-11 1:00 
Mean BI’ 0-301 0-072 0-033 0-18 0-65 0-057 0-125 1-00 


B corrects for the branching of the radioactive series at thorium C. 


seve Rm*c‘B4 ae 
o=8 (nape bP): 


The intensities /’ in quanta per second relative to the 2:62 Mev. y-ray have 
been calculated using a modification of the formula of Ellis and Aston (1930). 
Their relation for the height hy of the line from the K shell may be put in the 
following form, derived elsewhere (Richardson 1950): 


= constant x J’7, HpB? 


fn aan ee 
n is the number of atoms per cm? in the target of atomic number Z and J is 13-5. 
electron volts. Hpe/c is the momentum of a photoelectron of velocity v=fe. 
R=A(Hp)o/Hp where A(Hp), is the distance in the spectrum between the peak of 
a line and its high energy limit. 

As R is in the logarithm, its value is not critical. For our apparatus it was about 
1/40 for all the lines, as can be seen from Figure 3, where the lines are plotted with 
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the relative displacement in momentum as abscissa, and the fraction of peak 
height as ordinate. Identity in shape of the lines was shown by Ellis and Aston 
to imply that no large error will be produced by neglect of electron scattering in 
equation (1). The shapes can only be similar if scattering has reduced the heights 
of all the lines in the same ratio. 


NV (Hp) & (Hp) 


0 0-025 0:05 0:075 0:10 


oak (Hoy (1 Hp), 


Figure 3. Comparison of the shapes of the strong photoelectron lines. (Hp) is the 
high energy limit of the line. 


The table shows the correction for loss of B-rays in the mica window of the 
counter and also the fraction of the y-rays transmitted to the target surface. 
This ‘ pre-absorption’ is large for the soft y-rays and is discussed in an appendix. 
It is subject to error owing to uncertainty in the density of the thorium oxide 
carrier. ‘Tests with inactive samples of powdered oxide gave a most probable 
density of 2-7 gm/cm*. An error here of 25°% would give one of 10% in J” for: 


the F y-ray. 
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§4. LINES FROM THIN TARGETS 


Deutsch, Elliot and Evans (1944) have pointed out that K photoelectrons 
from a very thin target should give lines in the momentum spectrum of area A 
given by 

I’ =cA/rx, wie) of enetel (2) 


where c is a constant for the apparatus. Now rx as usually defined involves 
averaging the emission probability over all angles and (2) will only be true if any 
-anisotropy in the photoelectron emission is evened out by the apparatus. 


Table 2. Relative Intensities /” deduced by Equation (2) from 
the Lines from Thin Targets 


Target I’ with I’ with 
Line neice ei piece ra M line as X line as 
@ueiem’) transmission area of line Son Coen 
F 0:72 0:61 2050 0:49 — 
F 3-7 0:61 7860 0:58 0-70 
IL, 3-7 0°83 316 0-11 0-14 
M 0:72 0:84 447 0:65 = 
M 3:7 0:84 1400 0:65 0:78 
M 7:4 0:84 2090 0:65 0:59 
xs 3°7 0-93 141 0:83 1:0 
x 7:4 0:93 280 1:13 1:0 


Values of cA/r_x found with three thin gold targets are given in Table 2 with 
‘the constant c adjusted, first, to give the same J” for the strong M ray as the mean 
for thick targets in Table 1, and then, to give /”=1-0 for the X y-ray. 

The values of J” for the F and L rays are markedly less than those found with 
thick targets. Except in the case of the 2-6 Mev. ray the area of the line does not 
increase in proportion to the target thickness, indicating an apparent absorption 
-of photoelectrons coming from the deeper layers. 


§5. COMPARISON WITH y-RAY INTENSITIES DEDUCED FROM THE 
NATURAL £-SPECTRUM 

The values of J” found for the 0-238 Mey. F y-ray are less than that deducible 
‘from the known intensity of the f-transition to the 0-238 mev. level of ThC, 
which has been estimated to be 0-83 + 0-03 electrons per disintegration (Martin 
-and Richardson 1948). Subtracting from this the experimental number of internal 
‘conversion electrons per disintegration, which is 0-39 according to Table 3, we 
are left with 0-44+0-03 for the number BI’ of y-quanta escaping unconverted. 
This is 46% larger than the mean value 0-30 in Table 1 for thick targets and 2-2 
times the mean value for thin targets in Table 2. 

The larger discrepancy with thin targets suggests that a line of low energy 
suffers a greater reduction in area by scattering than does the comparison line of 
higher energy. ‘The intensities from thick targets should be more reliable 
because the effect of scattering seems to be effectively corrected in equation (1). 

Some of the 46% discrepancy with thick targets may be due to the variation 
‘with energy of 6, the most probable angle of ejection of the photoelectron from 
‘the atom. Fora K photoelectron @, measured from the direction of the incident 
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quantum, is 90° for a soft y-ray but, as Latyshev (1947) has shown, it falls to 9° 
for the 2-62 Mev. y-ray in lead, because much forward momentum is transmitted 
to the electron. ? 

The arrangement of source and target in Figure 5 may favour the detection 
of electrons emitted at small angles with the incident radiation. This would 
explain why the F line from the thinnest target of 0-7 mg/cm? is exceptionally 
weak because in such a target multiple scattering has little chance to even out 
the anisotropy. 

In the absence of information on the energy variation of the discrepancy 
Table 3 gives (a) lower and (b) upper limits to the nuclear transition probabilities. 
per disintegration p” where 


pi BL peop a ee 8) be hi: (3) 


The p® are the numbers per disintegration of internal conversion electrons from 
the K, L etc. shells. The upper limit (b) assumes that all except the 2:62 ev. 
y-ray have true intensities BJ” 46° higher than the mean value in Table 1, while: 
the lower limit (a) applies no increase to BI’. 


Table 3. Upper and Lower Limits (a) and (b) of the Experimental Internal Con- 
version Coefficients «, compared with Theoretical Values for (c) Magnetic: 
Dipole and (d) Electric Quadripole Radiation 


y-ray 1) G H L M O 2 xX 
Nucleus ThC 28Pb ThC %8Pb 28Pbh ThC’ 28Ph 2°8Ph 
pf Bis 2-64 117, 184 1-34-0106 0262, 0-178 
B-line strengths pe 5-7 0-42 02 ss 0-36 
he Spo) 388) es 2 — 1-44 2-18-1982 0127' 0-314 -0:014 
TEST ODS Ty 10 4-7 20 67 355) 100 
eet 1) 4 83 14 6-2 29 97 S50 10 100 
x, from (a) 45 26 24 $) 2:0 1:8 2) 0:18 
experiment (b) SY el 19 6:4 1:38 ilss} 1-4 0-18 
(c1) 32:5 24 21 71 
(iy 125) S5y5 7 22 DA 136) 00920015 
Ga fromtheory (9) 44 34 32 10:2 7:5 54 3:1 0-20 
(d2) 10 ee Deed et 5S 1-13 0-77 010 


Values expressed as percentages. 


§6. THE INTERNAL CONVERSION COEFFICIENTS ag 

In Table 3 the experimental and theoretical coefficients are compared, using 
the definition «% =p," /p’. 

The values of p’ for the H, [+Ja, J + Ja, Jb and N lines are from Flammersfeld 
(1939) and for the rest from Martin and Richardson (1948). Feather, Kyles and 
Pringle (1948) have found p’ = 28-8 % for the F line which gives a,(a)=43% and 
BD) = o4-7 Vos 

The theoretical values in rows (cl) and (d1) are from the curves of Hulme, 
Mott, Oppenheimer and Taylor (1936), while (c2) and (d2) are from the compre- 
hensive unpublished tables of Rose, Goertzel, Spinrad, Harr and Strong (1950). 
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We have assumed that the tables give the quantity «,°=p,*/BI’ and that 
Og = Oy BI’ /p’ =a49/(1 +1-24a,°). 

Both the experimental values (a) and (b) of «, for the F y-ray are larger than 
the 32:5°% for pure magnetic dipole (M.D.) radiation given by the earlier theory. 
The discrepancy is of opposite sign according to the new theory in which the 
value of 44°% in row (c2) is 17% higher than the experimental value (b) of 37:7%. 
The difference renders doubtful the hypothesis of Feather and Richardson (1948) 
that the ground level of ThC has zero spin quantum number J, because this makes 
the F y-ray pure magnetic dipole, the only radiation allowed by the spin and parity 
selection rules for a 1-0 transition. However, the difference is probably 
within the limits of error. 

Turning to the y-rays of 2°8Pb, the new theory indicates that the 2-62 Mev. 
(X) y-ray is pure M.D. corresponding to AJ=1—-0, rather than pure electric 
quadripole (E.Q.) as would be required by the alternative possibility that 
AI=2—0*. Thus the first excited state should have J=1. 

This result makes it possible to eliminate two of the three alternative sets of 
values of J previously suggested for the levels of ?°°Pb (Richardson 1948), all 
assumed to have the same parity. The remaining set of values is shown in Figure 
-4 which gives spin changes AJ=0, -1,—2, —2 for the G, L, M and P y-rays, 


Spin 
Energy 
in MeY. JE 
371 4 
3-48 3 
3:20 3 


2-62 


Figure 4. 


respectively. According to the new theory the experimental values of «, for the — 
G and L rays are definitely lower than for pure M.D. radiation. This is consistent 
with the admixture of some E.Q. radiation which is to be expected for transitions 
between levels of high spin. ‘The lower value of «, for the strong 0-582 Mev. 
(M) ray is in satisfactory agreement with the value for E.Q. radiation appropriate 
to the larger spin change of two units, which forbids the emission of M.D. radiation. 
The experimental values of a, for the 0-859Mev. (P) y-ray which also has 
AI=3-— 1 are, however, not in good agreement with the value for E.Q. radiation, 
but suggest a mixed radiation corresponding to AJ= +1 or 0. It is probable 
that the discrepancy is within the experimental error in measuring the weak 
internal and external conversion lines due to this y-ray. Similarly the weak 
0-726 Mev. (O) y-ray of ThC’ has a multipole character which may be either pure 
FQ» of mixed) 


> Weare indebted to Dr. B. B. Kinsey for drawing our attention to this point and to the existence 
-of the new calculation. ; 
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ALP PE NeDalexX 
THE CORRECTION FOR y-RAY TRANSMISSION LOSS 


In the general case of y-ray absorption multiple Compton scattering must be 
considered, but in our problem it can be ignored because y-quanta which have been 
scattered once will usually have lost enough energy to prevent their photoelectrons 
from contributing to the height h of the lines. 

We have for the total linear absorption coefficient 

PAIR, 
where 7 represents the total photoelectric absorption assumed to be 5/4 times the 
absorption in the K shell, « is the coefficient for pair production and o is that for 
loss by scattering from a parallel beam of y-rays, as given by the Klein—Nishina 
formula. « and o have been taken from the data of Sizoo and Coumou (1936) 
and Sizoo and Willemsen (1938). 

We assume that photoelectrons ejected from any two points on the flat 
surface of the target have equal probability of entering the counter so that the 
number counted from an element of area do of the surface is proportional to Ido, 
where J is the intensity of radiation at the element. J is reduced by absorption 
in the thorium oxide, the copper tube and the target. {Jdo is evaluated over the 
surface, firstly for the actual case with absorption and secondly for the ideal case 
with no absorption. ‘The ratio of the two integrals is the transmission term 
given in Tables 1 and 2. 

To evaluate [[do the volume source is replaced by five parallel linear sources 
whose mid-points are shown at A, B, C, D, E in Figure 5(c). The target 
surface is also broken up into seven equidistant lines parallel to the tube axis. 
Consider a source line m and a target line m a distance d,,,, apart (Figure 5(a)). 
Then J at Q due to an element dy, of the source line at P is proportional to 


{dy, exp (— Zp,d,sec 0)}/d2m sec? 6, 
where Xp;d;sec @ is the sum of the absorption exponents for the media between 
Pand Q. The rate N,,,, of emission of photoelectrons from the mth target line 
due to the uth source line is 


N 


constant 
nm — a 


lL, L 
| dy, | d0 exp(—Xp,d;sec6).  ...... (4) 
0 0 


Equation (4) has been integrated numerically for each value of nm for the 0-238 
and 0-582 Mev. y-rays. 

In combining the JN,,,,, the accuracy may be improved by using, where possible, 
a numerical integration formula instead of simple addition. ‘Thus the seven 
Nm for each source line n are combined by the Newton—Cotes formula, and 
Simpson’s rule has been used for combining the three discrete sources at 
B, E and D to find the emission due to the area marked 1 in Figure 5(c). 
The rule approximates to the true z-variation of the emission due to a source on 
the line BD by fitting a parabola through the three known points. 

The x-variation of the emission is almost linear, the value of Ny,, being near 
the mean of N,,, and N,,,. Thus the combined areas 2 give an emission propor- 
tional to that of a line source at B’ whose emission can be read off from the fitted 
parabola. A similar substitution can be made for the areas 3. Weights propor- 
tional to the areas have been used in combining the zones 1, 2 and 3. 


nin 
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The transmission found in this way did not differ by more than 3% from the 
value given for a single line source through E, the origin of coordinates. For 
the remaining y-rays therefore the transmission was interpolated by a more 
approximate method, equivalent to using a single line source through E. 
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Figure 5 (a). Figure 5 (6). 


Figure 5 (c). 


The source is first replaced by asingle radiant point at E and we find the number 
NN, of photoelectrons from the part abcd of the target in contact with the flattened 
part of the tube, where the absorber consists of parallel layers. Let ¢ be the angle 
between the ray from E and the axis (the normal to the target surface). Then 


N,=constant x | {do exp(X—p,d,secd)\/22sec2d  ....., (5) 


eG) ees) 2 ak wnt 1 
anddo=2rz*tandsec? ddd. Denote the integralin (5) by Zand put» = Xpjd,sec d. 
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Integrating over a disc between 0 and 4, 
T=2n|" e-"dn|n 
214d; 


SF NB be Lec faa 3) et Ae ra (6) 
where Ei{ —x] = | e~*dx/x, the exponential integral. 
x 


In the case of no absorption, J =27ftan¢dd =27Insec ¢. 

To find N over the rectangle abcd, J is found for discs of the three radii shown 
dotted in Figure 5(b) and the three contributions to N, are weighted in propor- 
tion to the fraction of the complete disc or annulus lying inside abcd. To extend 
the calculation from a radiant point to a radiant line, N, is next found for a single 
point at one end of the line. To do this J is evaluated for a fourth disc of radius 
equal to the length of the source and the contributions of the four zones are 
reweighted and added. The values of N, for sources at the centre and ends of 
the line are then integrated by Simpson’s rule to represent the line source. 

The number N, of photoelectrons from the parts of the target not in contact 
with the copper tube was found by numerical integration of 


g b eee meee 
oe (: “A cal Gayest ee (7) 


where a= Xp,d; for the target and b= Xp,d; for the source and tube. The result 
for a radiant point at the origin was extended to a line in the way used for Nj. 


dx dy 


N,=constant x [oq 
e+ yp4tez 
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Ellis and Aston’s Theory of the Height of a Photoelectron Line 
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ABSTRACT. The formula of Ellis and Aston relating the intensity of a y-ray to the 


height of its photoelectron line produced from a thick target 1s revised using the quantum- 
mechanical theory of stopping power. The effects of multiple scattering, bremsstrahlung 
and anisotropic emission of the photoelectrons are discussed and a general comparison of 
thick and thin target methods is made. 


Sale JON AT ROI UNE WON: 


quantitative method for estimating the intensities in a y-ray spectrum by 


LLIS AND ASTON in 1930 in a paper to be referred to as I, described a — 
E using the photoelectrons ejected from a thick target. ‘They gave the | 


following relation between the y-ray intensity J’ and the height h, of the photo- 
electron line recorded in a magnetic spectrometer : 


hg=l'c, pe f(s) < constant.) ee (1) 


Here hf, is the number of K photoelectrons recorded per second at the peak of 
the line, and J’ the number of quanta emitted per second. 7x is the photoelectric 
absorption coefficient for absorption in the K shell of the target atom. B=v/e, 
where wv is the speed of the electron of momentum p=Hpe/c. (8) is a slowly 


varying function of 8 taken from the classical theory of stopping-power of Bohr 
(1915). 


/(8) can now be given a quantum-mechanical form and it is of interest to see © 


if the additional stopping-power due to radiation losses may be neglected. ‘The 
error due to electron scattering can also be estimated; it is clearly small if the 
resolving power of the apparatus is sufficiently high, because then the peak of the 
line is formed by photoelectrons coming only from a thin layer at the surface of 
the target. Should the target be too thin to build up the peak to its maximum 
height, equation (1) does not apply and we have the alternative ‘thin target’ 
method. In both methods serious error may be caused by the variation with 
quantum energy of the anisotropy of photoelectric emission. 

Equation (1) may be of use in interpreting experiments on the natural f-ray 
lines from sources of low specific activity. The height of a conversion line is 
given by (1) if P’7, is replaced by P,, the intensity of the line in electrons per 
disintegration. For this application the source must be thicker than Ax, the 
‘contributing depth’, given by equation (5) below. 


§2. THE QUANTUM-MECHANICAL FORM OF EQUATION (1) 


(8) is related to the rate of loss of momentum of a photoelectron by collisions 
with the electrons in the target by 


— d(Hp)/dx=constant/ 67/(8\50 5 9) eee (1a) 
where x is measured along the path. An approximate relativistic expression for 
the collision loss in ergs/cm. is (Heitler 1936, p. 219) 

WE?v2 


d 3 me? 2 
ce! seas 2 AUS In Hmct +1- 8 | ailenelentette (2) 


‘ 
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Here n is the number of target atoms of atomic number Z per cm? and J is 13-5 
electron volts, so that /Z is an average ionization energy for the atom. 


by = (877/3)(e2/me?)? =6-57 x 10-25 cm? 


is the Thomson cross section and E( =emce?) is the photoelectron energy including 
the rest energy mc®. 

W is the largest single collision loss which contributes to dE/dx in (2); for our 
purpose it is chosen so that a larger loss could not affect the observed retardation 
of the peak of an initially mono-energetic group of electrons because it would 
remove the electron from the neighbourhood of the peak. Thus a collision is 
excluded if it causes a fractional loss of momentum exceeding R=A(Hp),/Hp, 
where A(Hp), is the interval in the photoelectron spectrum between the peak of 
a line and its high energy limit. A larger loss does not concern us because it is 
certain to bring the photoelectron to the low energy side of the line. 

From the relations p=)mce=mzv/(1 — f?)! and dp/dE =1/v, we get W=vRp 
and 

d(H, mc? Rm?*c'*64 
— = 8nZ$5 | In aa +1-2", ee 3) 


Writing F(R, 8, Z) for the expression in the square brackets and making use of 
(1a), we find that equation (1) becomes 


I’rx Hp B? x constant 


bg RZ) (4) 


The fraction R is a measure of the resolution of the apparatus and, in an 
instrument of constant p, is nearly the same for all the lines. It equals A(Hp)/Hp 
if scattering is absent, where A(Hp) is the width of the base of the narrow line 
formed by unretarded photoelectrons from the outer surface of the target. 
A(Hp)/Hp is fixed by the finite width of the target and slits. Its constant value 
may in practice exceed R(=A(Hp))/Hp) because the actual broad line may be 
sharpened by scattering which reduces the more retarded contributions coming 
from greater depths in the target. 


ave DE Poh wAvn © Pa Liha © ON TRIB UD ING SUR IOANC RIE Aer R 


If scattering is small the path x in the target will equal the depth of origin 
below the surface measured in the direction of the beam focused by the spectro- 
meter. Only electrons starting in a thin layer at the surface can contribute to 
the height of the peak because the others will be retarded by more than A(Hp)p. 
The maximum contributing depth Aw is 


d(Hp) _ 4 RnB? 
Ax = —A(Hp)y / Se Poa (5) 


where 7 = Hp/1704-2 and 3nZ¢,/4 =2-32 cm for gold. 

23-2Ax/7R is plotted against loge (=log E/mce?) in Figure 1 for a gold target. 
If the target is thicker than Ax the height of the line will not be increased but its 
low energy side will be broadened, thus obscuring the resolution of close lines, 
If possible therefore the target thickness should equal Aw for the hardest y-ray 
studied in the experiment and any underlying mounting of the target should 
have small Z in order to add only Compton electrons to the spectrum. 


16-2 


236 H, O, W. Richardson 


84) THEE PERE Cal © Has C Amer ReDING: 


It was remarked in I that if there were no scattering there should be little 
apparent absorption of photoelectrons coming from the deeper layers in the target. 
The theory of retardation and straggling would then predict a line shape in the 
momentum spectrum having a steep rise from the high energy limit up to a long 
plateau extending to low energies. In contrast, the observed lines are well 
defined and their steep fall on the low energy side must be mainly due to scattering, 
which would seem to operate so as to give them all nearly the same shape. ‘This. 
is clear from Figure 3 of the preceding paper (Martin and Richardson 1950), to 
be referred to as II, where several lines have been plotted with the relative dis- 


placement in momentum as abscissa. If the shapes are identical it follows that — 


the heights must all be reduced in the same ratio by scattering and the relative — 


intensities will be given correctly by equation (4). 
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Figure 1. Energy dependence of 23-2 Ax/7R (continuous curves) and of %8 (broken curves). %8 is 
the mean projected deflection in degrees, in gold (Z=79). The thicknesses Ax of the 


eens layer are such as to give mean retardations R=8(Hp)/Hp of 1/20, 1/40 and 
/100. : . 


It is of interest to see if the similarity in shape can be due to compensation of 


the greater scattering per centimetre of slow f-rays by their smaller contributing 
depths Ax. 


A convenient treatment of multiple scattering is that of Williams (1940) in 
which the angles are measured in a unit § given by 


d=2275(nt)i/B7e radians, 59) See (6) 


where 7) = e"/me”, ¢ is the thickness of the scatterer and Z’ =(Z? +Z)?, 


4 
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: The deflections are given projected on a fixed plane containing the initial 
direction. The mean projected deflection % (in units 8) is 
=O 80in)t P1456, oS (7) 


where M =0-647nt(Z’)**(h/me)?8-?, the mean number of collisions. The prob- 
ability of a projected deflection « is 


G(a) = (1 —17/2$9?)(2/r70m) €XP(—07/t7y)?), nee (8) 


where 4,=5-1(In M)!—4-0 and «,, =(% —7/#,)/(1 —77/2¢,?). 
q The formulae are only accurate for deflections less than 10° for which 
sinad ~a6, but they give a rough guidance at larger angles. 
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Figure 2. Energy dependence of the probability P of scattering through more than 90° (P(90)), 

and more than 72° (P(72)), in a path x giving retardations R of 1/20, 1/40 or 1/100. 


Putting ¢=Ax, from (5), in (6), it will be seen that there is some compensa- 
tion of the energy variation. We have 


SO RZ ee Me ge ee ee (9) 


Thus 6 varies quite slowly, nearly as <?. 

For the ideal case of zero straggling we may trace a narrow outgoing pencil of 
focused rays backward into the target. Then those rays having the same retarda- 
tion 5(Hp) will all have the same path-length x and, owing to scattering, their 
points of origin will lie in a cup-shaped volume having cylindrical symmetry 
about the pencil. Owing to the reversibility of elastic scattering, conditions 
inside the target will be the same as for a pencil of incoming rays each scattered 
in a total path-length x. The cup-shaped volume will extend rapidly with x 
and eventually part of it will protrude forward through the target surface. This 
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part, having no sources of photoelectrons in the actual experiment, will be non- | 
contributing. The fraction of rays excluded in this way will obviously depend | 
on the inclination of the focused pencil to the target surface, but it is clear thatsome | 


line sharpening will always occur and that it will be more pronounced when the 
focused pencil leaves the target nearly at grazing angles. 

There may be other less direct ways in which scattering reduces the low energy 
side of a line, but it seems likely that they will cause little variation with energy 
and that rays scattered out of the beam will often be replaced by rays scattered in. 
In view of the possible presence of such further factors, the probability P(90) of 
a projected deflection greater than 90° may underestimate the loss by scattering 
even for normally emergent pencils and so P(72) has also been plotted in Figure 2 
where 


POD). Gaye 


72°/6 


The values are given for three path-lengths x causing three retardations 6(Hp)/Hp 


of 1/20, 1/40 and 1/100 in gold (Z=79). 

In so far as extrapolation of (8) to such large angles is valid, the curves suggest 
that in an apparatus with the high resolution of 1/100 the scattering loss is. 
negligible, but for the low resolution of 1/20 it is quite large. Even for this resolu- 
tion it varies fairly slowly with energy so the heights of close lines will have the 
correct ratio. At 0-15 Mev. a retardation of 1/40 in gold may cause a scattering 
loss of 25%, whereas at 2-5 Mev. the loss is negligible. The decrease in line 
heights at 0-15 Mev. with a resolution A(Hp))/Hp of 1/40 will, however, be less 
than 25°, because the mean path x of the contributing photoelectrons is about 
Ax/2, giving an average retardation of only 1/80. 

A variation of 20°% in scattering loss should cause a small change with energy 
in the shape of the photoelectron lines near their peaks. Close study of Figure 3 
of II shows that the lines of the 2-6 Mev. y-ray are perhaps 0-7°% of Hp broader 
than the others in the region just below the peak. (The 0-47 Mev. M line in 
thorium also shows a small broadening which can be given another explanation, 
namely the presence of the unresolved weak line of 20kv. more energy ejected 
by the 0-51 Mev. y-ray from the L shell.) 


The error due to neglect of the small difference in the sharpening of the peaks - 


by scattering could be avoided by comparing the slopes of the lines on their 
high energy sides, instead of the heights. The electrons here come from such 
small depths that %3 will be less than about 20°. In the experiments of II no 
significant difference could be found between the ratios of the slopes and the heights. 


$3: LARGE COLLISTON TOSSES 


Collision losses greater than the critical value W, which occur in the contri- 


buting layer Ax, will remove photoelectrons which would otherwise contribute 
to the height h. 


The probability in Ax of a loss greater than W is given by the classical 
Thomson formula (Janossy 1948, p. 100, Bohr 1948). 


2/1 
P(W)Ax = 8nZ¢y aa (F es Fr) Agia i cau teen (10) 
max 


In our case the reciprocal of W,,,,,, the maximum possible loss, can be neglected. 
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Then since Ax = — W(dE/dx), we have from (2) and (3) 
ved (Oy Aaa NN Cal Ih Deceit ee ar ree (11) 
For R=1/40 and Z=79, P(W)Ax varies from 1/8-4 to 1/13-7 in the range 
hy =0-24 Mev. to hv=2-6Mey. Thus the fraction of photoelectrons originating 


in the contributing layer but removed in large collisions is about 8°/ and is too 
independent of energy to cause error. 


§6. BREMSSTRAHLUNG LOSSES 
At high energies, radiation losses make an appreciable contribution to the 
stopping-power of media of large atomic number. These losses must reduce h, 
first, by increasing the number of single losses smaller than W which cause the mean 
retardation of the peak, and second, for losses greater than W, by removing 
electrons which would otherwise contribute to its height. The stopping-power 
of the losses less than W is 


dE\ <¥ W/(E— me?) 
s (=) =n i kb, d{k(E—me2)}, ws. (12) 
X rad 0 

where k= hy and ¢,d{k/(E—mce?)} is the cross section of a target atom for a 
radiation loss between k and k +dk. 

The dependence of 4; on E and k has been given by Heitler (1936, fig. 14, 
p. 170). For k<W<E, the small losses are numerous and are due to collisions 
of large impact parameter for which the screening of the nucleus by the atomic 
electrons is nearly complete. For this case the quantity k¢,,/E¢ is independent 
of both the energies E and k. We have 


on Lee ate 228 ak? 
$= (SS) Jeo: lel Us em 


and 
hk, Eb =(4/3)(41n 1832-18 4 1/3) = 15-8 for Z=79. 


Integrating (12) 
i SV ie 
-( ) =15-8nd6EW|(E—me?)=A(e+1),  ....-- (13) 


dx rad 
where A = 15-8n¢dRmce?. 

For «> 1 the stopping power is effectively Ae and fore=1litis2d. Fore=6 
and R =1/40 the contribution of (13) is about 2% of that due to inelastic electronic 
collisions given by (2) so that the effect can be ignored except at very high energies. 

The second effect, due to large losses, is given by the probability P(W),.q Ax 
of a radiation loss greater than W occurring in the contributing layer Aw. 
We have (Heitler 1936, p. 171) 

Ai 


PUW mq Ax =x | $,,4{k|(E—me*)}. 


W/(E— mc?) 
An approximation for ¢, is k-1E¢ exp { —2-2k/(E —mce?)} which gives 
_ ngAxk i re, 14 
OMG) Siem imager cn ie ary eae zz (14) 
where z=2-2k/(E—mc?). 
For E=6me? and R=1/40 the probability (14) is about 2 x 10~ so that the 
reduction in height is negligible. 
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§7. THE THIN TARGET METHOD AND THE ERRDC Ob 
PHOTOELECTRIC ANISOTROPY 

Deutsch, Elliott and Evans (1944) have described a method in which the target 
is thinner than Ax so that hy does not reach its maximum value and is not given by 
equation (4). 

In this method the intensities of two y-rays may be compared by finding the 
areas of their lines in the momentum spectrum. However there is now no 
compensation to reduce the energy variation of scattering and the line of lower 
energy will lose more of its area than the comparison line unless «6 is kept small. 
In practice (In M)? is fairly constant so that it suffices to choose a target thickness ¢, 
giving 5 less than about 0-25 radians. For example, for the F-line of ThB in gold, 
if t=1-33x10-4cm. (2-6mg/cm?), we have from (6) that 6=0-228, giving 
P(90) =2%, and PBG2)=6%. 

For this thickness the scattering is negligible, however poor the resolution, 
but lines due to hard y-rays will be of very small intensity because of the decrease 
of the absorption 7, with increasing energy. 

In contrast, the thick target method gives a valuable intensification of the 
high energy lines due to the growth of the contributing depth Aw and, if R is kept 
small by using a high resolving power and a small target, the scattering error is 
again negligible. 

In both methods the anisotropy in the photoelectric emission will not be 
smoothed out under conditions of small scattering. 

It therefore seems desirable to use a collimated beam of y-rays and to replace 
Tx by the differential photoelectric cross section integrated over the solid angle 
focused by the spectrometer, as has been done by Latyshev, Companeetz, Borisov 
and Gusak (1940) in the case of Compton electrons. Where collimation is not 
possible the calculation should still use the integral taken over a limited solid 
angle because in practice the y-ray flux is not isotropic. 

When a thin magnetic lens is used the predominant flux is in the direction of 
the focused electron beam so that small angles of emission @ will be collected 
preferentially. ‘This intensifies hard y-rays relatively to soft ones, for which the 
probable value of @ is 90°. In the case of the 180° spectrometers of I and II the 
y-flux has a larger component across the beam, but even here there is little or no 
contributing flux in the backward direction so that hard y-rays giving forward 
momentum to the electron are still favoured. Thus in all practical cases the 
errors due to anisotropy and scattering both over-emphasize the hard y-rays. 
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ABSTRACT. Diffusion equations governing the fluctuations of nucleon cascades and also 
of electron—photon cascades are given. The electron cascade is treated at first in the usual 
approximation, but it is shown that the approximations are not essential for the treatment. 
The diffusion equations appear in the form of integro-differential equations which can be 
solved directly using numerical integration. It is hoped to present numerical solutions in a 
later publication. 


Sie EN DRO Dil Can ON, 


N two previous publications (Heitler and Janossy 1949) we have considered 

collisions of nucleons. A primary nucleon is considered to suffer a number 

of collisions; in each collision a recoil nucleon is supposed to be emitted 
and the primary loses some energy. ‘The number of particles thus increases 
by one in every collision and a cascade develops. In the present paper we 
consider the mathematical problem of the fluctuations of such a cascade. 

It will be seen that the problem of the fluctuation of the nucleon cascade is 
mathematically very similar to that of the fluctuation of an electron—photon 
cascade. We present a formalism for both types of cascades in the present paper 
and hope to give numerical results in a future publication. The theory of such 
stochastic processes has been outlined, for example, by Arley (1949). 

In §§2 and 3 we deal with generating functions of the actual distribution, 
in §4 we discuss the determination of the distributions from generating functions, 
and in §5 we show how certain restrictions in the assumptions can be eliminated. 


§2. NUCLEON CASCADE 


We assume, as in the previous publications, that the cross section for nuclear 
collision is a homogeneous function of primary and secondary energies. ‘Thus 
we assume 


dorsto koe. (1) 


y / E’ dk’ dE” 
w(Ey; E’, E")dE' dE" dx =w ( ) 


is the probability that a primary of energy E, suffers a collision along dx so that 


the resulting particles are in the intervals E’, E’+ dE’ and E", E" + dk”. We 


assume 
A, Re he ere: (2a) 


w(e’,<”)=0, eer Mire a! 0) od TE Speen (25) 
and thus the total collision probability per unit path is 


“1 pl 
| | wle',<’)de'de"=4. en ee (3) 
040 


The cascade shall be described by a function 


h(E, 2, Moy X). 
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This function gives the probability that after a depth « a primary of energy 1a, 
has given rise to 7, particles energy above «Zp, and n, particles energy below «Ey; 
for brevity we describe the m, particles as ‘fast’ and the m, particles as ‘slow’. 

To obtain a diffusion equation for ¢ we consider the first collision the primary 
particle suffers. The probability that the primary travels a distance « —€ without 
collision is exp{—a(x—€)}. The probability that along the next interval dé 
there should be a collision dividing the energy into amounts inside intervals de’ 
and de" is dé w(c’,<") de’ de". In order to end up with a cascade of m, fast particles. 
and n, slow particles, the particles of energies «’ and «” have to develop into 
cascades giving rise to n{ and n‘ fast particles and nj and nj slow particles with 


N / u” 
ni tnj=n, and ny +ny=Nyp. 


Summing over all possible depths of dé and all possible divisions of the cascade 
we obtain 


na -0O 7 oO € . 
bleomy mx) =| exp(-a(w-}.dé = | | b(Ssmimsé 
Y 0 ntni=m/0 “0 € 
NE+NZ=Ne 


iS 1 7] ” N 
<b (Somiont é) we, €") de der. nal (4) 
We introduce a generating function 


Gle,x, U)=E6(NyxUN eae (5) 


in which we have written, for the sake of compactness, U for u,,u,; N for 1, 193. 
UO for us 


From (4) and (5) we find 


Conte Hi exp{—a(w—£)} Af i ‘ G (5. é, vu) G (S ae v) w(e’,<") de! de”, 
Lae (6) 


If we multiply (6) by e%*, differentiate with respect to x, and multiply by e°*,. 
we find 


G(e,%, Uy +aGlex, U) = | 


| G(S.% v) G(S. vu) sole’, <”\\de' de", 

0. 20 € € 

where'G,=0G/0x. 1) =) (7) 
Equation (7) is an integro-differential equation which can be solved by 

numerical integration, provided an initial condition is imposed. The initial 

condition is obtained from the assumption that at x=( there is one particle of 

energy Hp, i.e. with « =1, and no other particle. Thus 


{ 8(1—m,)8(m,) for «<1, 
(€, N,N, 0) — 1 ; : 
| 8(m,)0(1 —my) for «>1, 


f 1 
and, therefore, Ge, uy, Uy, 0) = | ee eS (8) 
(ig fone = 1. 


The solution of the system (7) and (8) can be found by the well-known methods 
of numerical integration. ‘The procedure is quite straightforward, except that 
for fixed values of the w, and u, one must compute for all values of ¢ ranging from 
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0 to co simultaneously. This procedure is greatly simplified by the fact that for 
e>I1 G(<, U,x) can be obtained explicitly. We note that there are no particles 
with «>1 and we can therefore write G(e, x, U)=g(x, U) for e>1. Introducing 
this into (7) we find 


£2(x, U) +ag(x, U) =ag?(x, U) for «>1; 


thus 
1 Le ( 1 )p1 
e(x,U) &ex\2(x, U)) 
and 
1 
ae poe a ee ae oy weet eee 9 
(x, U) = za (9) 
From (8) and (9) we have 1/(A +1) =u, and 
e Xe 
— SS Se ee C2 ae «6 8p 56 6: ce Le es 1 
Gie>1, x, U) =2(x, U) Paatiee) (10) 


Equation (10) gives thus the values of G for e>1 and only the values «<1 have to: 
be determined by numerical integration. Developing the right-hand expression. 
of (10) into powers of uw, we find incidentally 


b(e> es ny ==): Ng =n, x) SE ee LOAVES 


The above expression gives the probability of finding exactly m particles (of any 
energy) at x. 


SS DEE CRON SOAS Cube 


1. We describe the electron cascade by a function $c, N,«) where 7=1 for 
a primary electron, 7 =2 for a primary photon, N=n, mo, M3, M4; 11,2, = number 
of electrons of energy >,<eE,, m3,n,=number of photons of energy >, <eK. 

P(e, N, x) is the probability that a cascade of N members be found at x. 
Considering the ‘approximation A’ for the moment (e.g. full screening and 
negligible ionization loss), the function ¢ does not depend explicitly on the primary 
energy Ey. For the cross sections we use the well-known Bethe—Heitler expressions. 
and write w)(e) de dx for the probability of an electron E to emit along dx a photon 
of energy Ee in the interval Ede, and w®(e) dedx the corresponding expression 
for the emission of a pair one partner having the energy <k. We write 


a] - 

erode Soma, ie 0 ae (11) 

J0 
Because of the divergence of the Bremsstrahlung cross section the integral 
_ diverges for7=1; we avoid this purely mathematical difficulty by postulating 
W(e\=On Ome 3 ee ae (12) 


where k is a small quantity. In all physically interesting applications it will be 
possible to use the limit kR+0 for the solutions. 
The diffusion equation can thus be written 


fe, N, x) = ie exp{—a(x—&)} 2% * iy ow (5 ; Me é) 


N+ N= 


x A®) (; = NE é) w0(e') de’ dé. 
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Introducing the generating function : 
Ge, U, x) 2 Lde,N,x)U*%, sss ee (13) 
N 
where U =1y, Ug,llz, Ug, etc., we have 


‘Ge, U, x) 
“Xe rl ; 
- | exp{ —a,(x—é)} dé | 60(5, UG, é) Goo; Sy; G; é) wi(e') de’, 
0 / 0 a 


and differentiating in the same manner as before, 


Ge, U, x) 
= [,{or( U, *) oo, U, *) — Gee, U, »} w(e') de’, 
0 é Gs: 
where again-GO =0GO/0x.,- ~ 1 (15) 


We have replaced «; by Jw(e’) de’ and transferred the term «,G® from the 
left to the right-hand side and placed it under the integral. 

2. The equation (15) can be integrated for a given initial condition. We note, 
however, that owing to our approximations the number of electrons or photons 
below a given energy is infinite in the limit k>0. Physically only those states are 
of interest which contain a given number of electrons or photons above an 
energy «£, together with an arbitrary number of electrons and photons below 
Ey. We introduce therefore the probabilities 


ioe) 
Qe, Ns x)= eee, My ttle) 
V2, na=0 


“The corresponding generating functions are 
g(6u,,2)=  E $en',n", xu’o”” = Gc, u, 1,0, 1,4); 
W,n’=0 
the generating functions gare special cases of G, and thus also obey equation (15). 


To obtain the initial condition for the g we note that for x =0, 7 =1 there is one 
electron present, while for x =0, 7 =2 there is one photon present. Thus 


$Me, 10,4 =0) = dX 01 xe =O) ak exal, 
EM(e, 0,0; 4 =0) =h2X(e70;0;v=O0)= 1 e>1; 
all other 4s are 0 for x =0, and, therefore, 
EU Cy ene AUN eh esle') toa siete ree (16) 
and EVM, Uy, Ua, X= 0) =P, u,; 2, % =0) = Uy eslie eee (17) 


As at any depth there are no photons or electrons with energies exceeding the 
primary energy, we can generalize (16) to 


2%, uy, Us, x)= 1 ESSil. on Seba ene (18) 


Examining (15) we find, with help of (18), that the integrand remains finite when 
approaching «’->0, (k=0) if the specialized gs are used. Thus, considering 
electrons and photons above «Ko, the right-hand side of expression (15) has 
a finite limit for k>0. g® is the generating function of the simultaneous 
distributions of electrons and photons above «Ey. We may average over the 
photons or over the electrons and thus obtain the distributions of electrons 
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accompanied by an undetermined number of photons or the number of photons. 
accompanied by an undetermined number of electrons. We write 


$Me, n, x)= EU Pen, v, x) 
v=0 


CO 
and pis Ne, n, x) = PE, U;N, x) ; 
v=0 


the four functions $“” 7,k=1,2 represent thus the probabilities of finding 
given numbers of electrons or photons in cascades initiated by electrons or 
photons. The first index refers to the nature of the primary, the second to the 
nature of secondaries under consideration. The corresponding generating 
functions are 


ge, u, x) = > DEN ESIC re ne ie ne eee (20) 
n=0 
R= 12. 


The g“ are again special values of the G®, namely 
PON, ue) = Ge, u, 1,1,1,.), 
Pee Ux) = Ge, lolyu, 1x). 
$4. DISCUSSION OF THE GENERATING FUNCTION 
1. From the definitions (15) and (21) of the generating function, it follows that 
eee x=, AU ES ald RE ee (22) 


This condition is compatible with the diffusion equation (14) or(15), Differentiating 
equation (20) with respect to u, we find for u=1 


(i, k) oO ; 
(ee) r= DS up Be, Vv, x) 
Ou u=1 v=0 


ATG mec) ek Eee! chaste (23) 
The latter expression signifies the average number of electrons or photons with 


energy >cE at a depth » given rise to by an electron or photon, Similarly we 
can get higher momenta of the distribution 


OP ge, u, x) 
bi is SEN Pea =P Me,v, x), «ww vee 23a 
( d”(log u) Gath ego) ( ) 
or mixed momenta 
OP +I MVe, u, V, =) —Sy'Py"4 ) v0". x 23b 
Geeaieo ee ae (236) 


We note that differentiating (14) and introducing (23) we have, with the help of 
(22), replacing G by g, 


S uli, e,x)) =— (n(t,k; €,x)) ie w(e’) de’ + |, (n(1,k; =; , x))z(e’) de’ 


ile" 


i | GOs SAI eee (24) 
0 


The above equation is equivalent to the well-known diffusion equation of the 
theory neglecting the ionization term. 
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Similarly, the second derivative with respect tow of (15) gives the diffusion 
equation for the mean square spread; such equations were investigated by 
Nordsieck, Lamb and Uhlenbeck (1940). 

2. The fact that the equations (15) can be reduced in a special case to the 
ordinary diffusion equations for average numbers sheds some light on the amount 
of work necessary to integrate the system (15). To build up the solution of (15) for 
a fixed value of u is about equivalent to integrating numerically equation (24). 
Such an integration is somewhat cumbersome as it is necessary to compute 
simultaneously solutions for all values of ¢ above the required value of «. However, 
for the purpose of a complete tabulation of the functions (n(7,k; «,x)) etc. this 
procedure is not essentially more elaborate than it would be to integrate some 
ordinary differential equation numerically and to get solutions for all values of 
e one after the other. 

We are of the opinion therefore that for the purpose of tabulating the cascade 
functions comprehensively the direct numerical integration of the diffuse equations 
is not essentially more elaborate than computation by the well-known analytical 
methods. 

To tabulate the generating functions a procedure similar to direct numerical 
integration of, say, (24) has to be repeated for several values of the argument U. 

3. The actual probabilities de, 2’, n”, x) or 6“ (e,n, x) can be obtained for 
small values of n by differentiating the generating function. ‘Thus 


1 gu tn” o® U.U.X ; 
“| ( ic ea ; = AME, ite (Dee x) 
i u=v=0 


n\n eu’ avn" 
1 (# Meu 
n! ou” en fee) 


The derivatives can be obtained either by numerical differentiation or alternatively 
by successive integrations with respect to x of the derivatives into u of the 
equation (15). Dzufferentiating (15) times into u, we get a recursion formula 
which permits determination of the mth derivative of g with help of the lower 
derivatives. As an example we consider a case, of particular importance in 
practice, namely, > Ye, 1, x), 1.e. the probability of an electron passing through 
an absorber and emerging from it accompanied by photons and low energy 
electrons only. We have for this probability 


099 D(e, u, x) 
Wie oe) [poe 
p& Me, 1, x) ( By i : 


The diffusion equation for this quantity is obtained from (15) and its derivative 
into wu with the following initial conditions (see (16), (17) and (21)): 
0 eal 
ot Me = OO) = as 
| tel 


ge Ne t=O (i= 1 any € 


(eet) 1 we<<l 
du u=0 Ome 1 


(ee ma, 
eo We = any «. 


Note on the Fluctuation Problem of Cascades 247 


4. For large values of n it is impracticable to evaluate the high-order derivatives 


of the generating function step by step. Such cases can be dealt with using 
complex integrals. We have 


rn 


Pp Me, nm, x) = = > u-"lolts Me, u, x) du, 
and similar expressions for other probabilities. ‘The complex integrals can as 
a rule be deformed to a path parallel to the imaginary wu axis and the integrals 
can be evaluated by the saddle-point method. This procedure requires the 
values of ge, u, x) and both their first and second derivatives into uw. 

Similarly probabilities of the simultaneous presence of specified numbers 
of electrons and photons can be obtained by double complex integrals, which 
can be reduced to double saddle-point integrals. We hope to return to the 
details of this procedure in a later publication. 


§5. EQUATIONS FOR ARBITRARY CROSS SECTIONS AND 
IONIZATION LOSS 


1. In the previous sections we have assumed that the differential cross section 
is a homogeneous function of the energies, and we have also neglected ionization 
loss. Neither of these assumptions is essential to the formalism ; these assumptions 
have, however, the advantage of reducing the number of independent variables. 
Instead of the two variables of the exact procedure, namely the primary energy Ey 
and the final energy EL, we have in the approximate procedure only their ratio « 
as an independent variable. 

If the simplifying assumptions are dropped we have thus to deal with 
functions and generating functions * 


PO Eo, E, Nx), GO(Eo, E, U, x). 


The diffusion equations for a general cross section w(K, E’, E”) in the case of 
electron cascades without ionization loss is thus 


GO(E,, E, U,x) = { exp {—0,(E,)(x —)} dé i | GE", E, U, é) 


0 
x GEE", EU, £) (Ey; B'E")dE'dE" as (25) 
with a,(£) = | : | WE ; es a) dE! Ai Be 
0-0 


We integrate separately over E’ and E" as for actual collisions £’ + E” is not exactly 
equal to E on account of recoil of the nucleus. 
From (25) we get the differential equation 


G(Ey, E,U,x) =| | (G2(E’, E, U, x) G?-(B", B, U, x) — GEy, E, U,x)} 
0 Jo 
Kiba: he jd dk, pagesia.20) 
The equation (26) can be integrated in the same way as (15). For a fixed value 


of E one has to compute simultaneously all the values of Ey between E and the 
desired value. 


* For the electron. cascade i=1, 2 most of the following formulae apply to the case of nucleon 
cascades when the index 7 is omitted, and accordingly the number of equations reduced. 
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2. Ionization loss can be included into the formalism by remembering that 
the primary of energy Ey has at the time of its first collision at x—¢ an energy 
E,—B,(x—€), where f, is rate of ionization loss of an electron, B,=0. ‘Thus we 
have in the case of homogeneous w, i.e. for «(E)) =a; (independent of Ep), 


GO(By, E, U, x) = | exp {—a(v—2)} dé | 5 GE’, E, U, )G°-(E", E, U, €) 


x w(E, — B(x —€), EB’, E") dE’ dE’. 
Multiplying by exp («;x) and differentiating with respect to x, we find in the usual 
way 
GE HU = fe | { GOES, ES U,x)GO OE EO, 4) Ge hank, Ona) 
0:0 


@G(Ep, E, U, x) 
; dE, 


3. As the operation 0/0E, for large values of Ey reduces the order of w, the 
second term on the right-hand side is small for E,>,; thus it is small provided 
the primary energy much exceeds the initial energy. 

An approximate solution of (27) is the solution without ionization loss. 
Denoting this by Gp, we put 


GE), E, U, x) 7 G}(Eo, E, by, x) + By gi(Eo; E, U; x) + Byg}(Ep, 1 OG x) 3h oa 
In the first order of 8 we have from (27) 


ORT ond Ua) aCe alee [| topes e U, *) 
OFZ0 


x w(By, E’, BE”) dE’ dE" —B 


Ox Slee OE, 
xP OE’ E,.U, x) Po (EB, Ue )nG ie keia) 
—s) (Ey, E, U, x)\u( Ey EE) dE de" | ee (27a) 


The equation (27a) can be integrated if Gf? is known. The integration requires 
the same amount of work as the integration of (27) itself, which gives the exact 
solution directly. The practical advantage of (27a) over (27) is, however, this: 
If w(K, FE’, E") dE’ dE” is a function of E’/E, and E”/E, only, then G, is also a 


homogeneous function; further, writing 


Go, £, U8) =Gileteay e=L/E,, 
we have 
ae — m -F x function of e, 
Thus 
function of (£/£,) 
i : 
and g‘ needs to be evaluated for various values of the ratio of E/E, only. 
4. Finally, for the sake of completeness, we give the diffusion equation, 
including ionization loss for arbitrary w(Ep, E’, E’). 
The probability for a primary to lose energy by ionization only down to a 
depth x —€ is 


gE, E; U, x) a 


PO(Eo, x — €) = exp { a ie a,(Ey — Bin) an} , 
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Thus the integral equation for the generating function may be written 
GO(E,, E, U,x) = | Sag { = pee qb) an} dé | ‘ | ‘4 G(E’, E, U,£) 
x GOH, EU, &)o(E,—B(x—€)E’, EB") dE’ dE’. 22... (28) 
Multiply (28) by exp { | aE Ba) an} 
differentiate, and divide by the same factor. This procedure gives 
See) - (e f GY(E', E, U,«)G°-°(E", E, U, «)w(Ey, E', E’) 
— GO(E,, E, U, x)w!(Ey —B,x, E', E") dE’ dE" 


2 OG Ey, E, U, x) 


aE Sard Shee See (29) 


Equation (29) can be integrated numerically for a given initial condition, as all 
previous equations. Again, asin the preceding section, it is possible to approximate 
the exact solution by a series of homogeneous functions. 
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ABSTRACT. An account is given of a study of 2,150 stars in Ilford C2 plates and 312 stars 
in Kodak NT4 plates exposed under different conditions. Possible values for the cross 
section for star-producing collisions are discussed. It is concluded that this collision cross 
section is equal to the geometrical cross section of the nuclei, and that the stars producing 
particles are non-catastrophically absorbed, leading to the observed absorption cross sections. 
It is shown that nuclear transparency could not account for the observed results. It is 
estimated that 17% of all stars and 40% of stars from which two or more lightly ionizing 
particles are emitted are caused by ionizing primary particles. ‘The size distributions for 
stars caused by ionizing and non-ionizing particles are given, and possible explanations for 
the difference between them are considered. The experimental observations are in agree- 
ment with the theory that the star-producing particles are mainly neutrons and protons 
secondary to the primary cosmic rays. 

The variation with star size of the proportion of doubly charged to singly charged particles 
is found, and for stars of more than six prongs is in agreement with Le Couteur’s extension 
of the evaporation theory. A value of 0°54 is deduced for the proportion of doubly charged 
particles from disintegrations in the light elements of the emulsion. 


An increase in the proportion of lightly ionizing tracks is observed for stars of more than 
eight prongs, and is ascribed to the creation of mesons. 


SoS IN PROD CC REON 


HE fundamental problems requiring elucidation in the study of nuclear 

disintegrations caused by cosmic rays are the nature of the star-producing 

radiation and the mechanism of disintegration of the nucleus. It is now 
thought that the radiation responsible for most low energy stars (i.e.< 200 Mev.) 
consists of neutrons, and for high energy stars fast neutrons and protons, with 
a cross section for star production proportional to the geometrical cross section. 
The disintegration of the nucleus is thought to involve the emission of a number 
of relatively slow particles with energies and an angular distribution typical of 
evaporation, and a smaller number of particles of high energy (i.e.> 50 Mey.) 
which form a ‘tail’ to the evaporation spectrum. These particles probably 
consist of the incident nucleon, knock-on nucleons and mesons created in 
nucleon—nucleon collisions. 

A study of 2,450 stars formed in Ilford C2 and Kodak NT4 nuclear research 
emissions has recently been completed and the main results are summarized in 
this paper. ‘The conclusions regarding the nature of the star-producing radiation 
confirm those of other workers. The results on the mechanism of disintegration 
of the nucleus afford a test of Le Couteur’s (1950) extension of the 
evaporation theory of nuclear disintegrations, and where it can be applied this 
theory is found to be in substantial agreement with the experimental results. 
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$2. EXPERIMENTAL DETAILS 


Details of plates, exposure times etc. are given in Table 1. 


Table 1 
Type of Thickness of Area Altitude Duration of | Number of 
plate emulsion (j) (cm?) (m.) exposure stars 
Ilford C2 100 200 3450 75 days 2000 
Ilford C2 100 160 12000 2 hours 150 
Kodak NT4 200 51-4 3450 14 days 192 
Kodak NT4 400 36 2800 6 days 120 


The large group of 2,000 stars was obtained from a group of plates of which 
one-third was placed horizontally, one-third vertically, and one-third vertically 
under 10 cm. of lead. All the other plates were exposed vertically. Moreover, 
all the plates were kept at 0° c. in thin air-tight containers, and they were developed 
within a few days of the end of the exposure. Except for the Ilford plates exposed 
at 3,450 m. the fading was negligible, but even in this case it was estimated that 
less than 3% of the stars were lost. These plates were fully developed to 
give uniform development to the bottom of the emulsion, and ionizations down 
to 3-5 times minimum (corresponding to protons of energy 100 Mev.) could be 
detected. The plates exposed at 12,000 m. were under-developed to facilitate 
particle discrimination and energy measurements; for these plates ionizations 
down to 7 times minimum, or protons of 40 Mev., were recorded. ‘The electron- 
sensitive plates were developed by the freezing method of Dilworth, Occhialini 
and Payne (1948) and minimum ionization tracks were easily detected, their 
average grain density being about 400 grains per mm. 

All the plates were scanned with great care and checks indicated that less 
than 1°% of the stars were missed. ‘The electron-sensitive plates were scanned 
with oil immersion objectives. Radioactive stars and possible cases of scattering 
were excluded. 


§3. THE STAR-PRODUCING RADIATION 
(i) The Total Flux at Different Altitudes 


The Kodak NT‘4 plates exposed at an altitude of 3,450 m. under no absorber 
recorded 13-4+1 stars/cm3/day. Of these approximately 1-0/cm/day, or 8%, 
could not have been observed on the C2 plates because they would have appeared 
as single tracks or large angle scatterings. Since the fading on N’I4 plates was 
negligible and the scanning very carefully carried out, it is justifiable to assume 
that this figure represents the true rate of production of stars. The rate of 
production of stars in the Ilford C2 plates would be expected to be about 8% less, 
but in fact the measured value was almost the same, i.e. 13-5 +0-5/cm?/day for 
plates exposed under the same conditions. It is assumed that this indicates that 
fading was small and the scanning adequate. The figure of 13-5 + 0-5 stars/cm?/day 
may therefore be taken as the rate of star production at 3,450 m. Moreover, the 
number of stars with three or more prongs was found to be 12:2 + 0-6 which is 
not appreciably different from the figure of 14-2+0-5 given by Cortini and 
Manfredini (1949). 

From the observed rate of star production at 3,450 m., it is possible to calculate 
the flux of the star-producing particles at this altitude provided that the cross 
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section for star production in the emulsion is known. Consider one component r 
of the emulsion, with density p, and atomic weight A, of which there are NV, atoms. 
per cm’, Then the number of stars dn, formed in the atoms of the element r, 
per day, in 1 cm? of emulsion dx cm. in thickness is 

dn 

— =N,o,dx, 

No 
where n, is the total flux of star-producing particles/em?/day, and a, is the cross 
section for star production in the element 7, Summing over all the components. 


of the emulsion we have, since NV, = a , where N, is Avogadro’s number, 
dn oP» 
a =o A, dx,  ) ghetto ee (1) 


where dn is the total number of stars per day in 1 cm? of emulsion dx cm. thick. 
In this summation the hydrogen in the emulsion is excluded, since collisions. 
with hydrogen atoms will not in general give recognizable stars. 

Before calculating the flux of star-producing particles it is therefore necessary 
to determine the values of o,. Now the cross section for absorption of the 
star-producing radiation has been found for air, carbon and lead by George and 
Jason (1949) and for lead and aluminium by Bernardini et al. (1949). In all 
cases it is equal to half the geometrical cross section of the absorbing nuclei 
(taking the nuclear radius to be 1-42 x 10-134"? cm., where A is the atomic weight). 

An increase in absorption cross section with atomic weight may be obtained 
in two ways. In the first place the cross section for star-producing collisions. 
may be equal to the geometrical cross section of the absorbing nuclei. This 
implies that the mean free path of the star-producing particles in nuclear matter 
is small compared with the nuclear radius. ‘The fact that the absorption cross 
section is only half the geometrical cross section would then be accounted for if 
the collisions were not catastrophic, each star-producing particle causing on the 
average two disintegrations. 

The alternative explanation is that the nucleus is not completely opaque to: 
the star-producing radiation, that is, the mean free path in nuclear matter is. 
of the same order of magnitude as the nuclear radius. In order to explain the 
quantitative results it must then be assumed that the primary particle passes. 
on the average through two nuclei before producing a star, giving a cross section 
for star-producing collisions equal to half the geometrical cross section. The: 
absorption and star-producing collision cross sections would then be equal. 
This means that the mean probability for the primary particle to pass through 
a nucleus without being absorbed is equal to one half. It follows from an analysis. 
first given by Cocconi (1949, equation (7)) that if this probability is a constant 
for all elements, the mean free path of the star-producing radiation in nuclear 
matter must be directly proportional to the nuclear radius of the element. This. 
is most probable and it is therefore not possible to explain the experimental 
results satisfactorily by nuclear transparency. We are then left with the first 
explanation, namely that the cross section for star-producing collisions is equal 
to the geometrical cross section. This will give a value of the flux which differs. 


by a factor of two from that calculated using the absorption cross section (see 
George and Jason 1949). 
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Equation (1) then becomes 


dn —13)2 Pr 

ae o7( 1:42 x 10-38) An 4 ae 
Using the observed values for dn the total flux at 3,450 m. is found to be 
366 particles/cm?/day, and at 12,000 m., 11,000 particles/cm2/day. If exponential 
absorption is assumed the total flux at sea level is 28 particles/em2/day. 


(11) The Nature of the Star-producing Radiation 


It is generally agreed that the majority of star-producing particles at altitudes 
below approximately 4,000 m. are energetic neutrons, whilst the charged particles 
observed to cause stars may be protons or o mesons (George and Jason 1949), 
Using the NT4 plates exposed at 2,800 m. and 3,450 m. an attempt was made to 
find the relative numbers of stars formed by ionizing and non-ionizing particles, 
and to discover any characteristic features of the two types of star. There is 
inevitably some uncertainty in identifying a particle as the primary, but although 
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there are doubtless errors in individual cases, the general conclusions are believed 
to be valid. The criteria used to distinguish an ionizing primary were somewhat 
similar to those outlined recently by Brown et al. (1949 a), i.e. tracks less than 90° 
to the vertical, on the axis of symmetry of emitted fast particles and with ionization 
less than 1-5 times minimum. It was found that the angular distribution of 
ionizing primary particles selected by these criteria obeyed a cosine squared law, 
within the limits of statistical errors. 

In agreement with other workers it was found that ionizing primary particles 
caused 17+3% of all the stars, and 40% of the stars from which two or more 
lightly ionizing particles were emitted. If it is assumed that the cross section 
for star production is the same for charged and uncharged primary particles 
this means that at 3,450 m. the flux of charged star-producing particles is 
62/cm?/day, and uncharged is 304/cm?/day, or 1:-4°% and 6-6% respectively of the 
total charged corpuscular radiation at this height. 

The differential size distributions for stars caused by non-ionizing and 1onizing 
particles are shown in Figure 1, and are in good agreement with the results of 
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Brown et al. (1949b). This observed difference in-size distribution may be due 
to the two types of star-producing particle having different energy spectra. 
Alternatively the results might be accounted for if the cross section for star- 
producing collisions were different in the two cases or varied with the energy of 
the incident particle in different ways for charged and uncharged particles. 
Thirdly, the probability of a certain energy dE being transferred to the nucleus 
by a particle of energy E, whilst traversing a distance d/ through the nucleus, 
might not be the same for charged and uncharged primary particles. 

If the non-ionizing primary particles are neutrons and the ionizing are, for 
example, o mesons, it is probable that the observed size distribution can only be 
accounted for by considering all three factors. However, if the star-producing 
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particles are assumed to be mainly neutrons and protons the last two explanations 
are not so plausible as the first one, i.e. a difference in the energy spectra. The 
proton spectrum would be expected to contain a greater proportion of high energy 
particles than the neutron spectrum since the protons with energies less than 
10° ev. lose energy very rapidly by ionization. 

The size distribution of stars on Ilford C2 plates at 3,450 m., is shown in 
Figure 2. N (>n) is the number of stars with more than n prongs, and A is a 
constant. It appeared from the size distribution at 3,450 m. and 12,000 m. that 
the proportion of large stars increased at the higher altitude. Thus at 3,450 m. 
the proportions of stars with more than 9 prongs is 0-045 + 0-005 whilst at 12 000 
it is 0-100 + 0-028. ‘This increase would be expected if the average energy of the 
star-producing nucleons increased with height. 
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§4. THE MECHANISM OF DISINTEGRATION OF THE NUCLEUS 
(i) The Evaporation Theory 


As stated in §1 the disintegration of the nucleus involves the emission of a 
number of slow particles (i.e. <50 Mev.), and a smaller number of fast particles, 
which may consist of the primary nucleon, knock-on nucleons and created mesons. 
The slow particles arise as a result of the excitation or heating up of the nucleus 
by the fast particles and, as is well known, their emission is analogous to 
evaporation. ‘This evaporation theory for the production of the secondary 
particles was due to Bohr and Kalckar (1937) and extended by Bethe (1937) and 
Weisskopf (1937). It was applied by Bagge (1938, 1946) and in greater detail 
by Harding et al. (1949) to the low energy particles observed in stars. They 
found, using the Fermi gas model of the nucleus and making due allowance for 
cooling and barrier penetration, that the theory gave a good account of the energy 
spectrum. ‘The evaporation theory has recently been extended further by 
Le Couteur (1950) by taking into account the thermal expansion of the nucleus 
and the effect of the neutron excess. Le Couteur has obtained the energy 
distribution of evaporated particles and has calculated the probability for the 
emission of protons, neutrons, deuterons, tritons, He, and «-particles as a function 
of the excitation energy. Since the excited nucleus is assumed to be always in 
approximate thermodynamic equilibrium the calculations have not been extended 
beyond energies of the order of 600 Mev. ‘Thus the maximum star size to which 
the theory can be applied is about 14 prongs. Moreover, since this model can 
only be applied to heavy nuclei, stars which originate in the light nuclei, carbon, 
nitrogen and oxygen, should be excluded. This means the exclusion of all stars 
of less than about 6 prongs. ‘Thus, strictly, the theory may be applied only to 
stars of between 6 and 14 prongs, but in actual fact a short extension of the limits 
is permissible. 

The proof that the light nuclei do not contribute appreciably to stars of more 
than 6 prongs is as follows. (1) The plot of the logarithm of the integral size 
distribution against star size (i.e. number of prongs n) consists of two straight 
lines with a clear break at approximately n=6. (2) The size distribution of 
50 stars produced in a cloud chamber filled with air, water vapour and alcohol, 
by a beam of 90 Mey. neutrons by Horning and Baumhoff (1949) in Berkeley 
shows a negligible proportion of stars of more than 5 prongs. It is reasonable 
to assume that the size distribution of stars formed in carbon, nitrogen and oxygen 
does not vary considerably with the energy of the incident particle if this is in 
excess of 90 Mev., since these nuclei are relatively easily broken up into «-particles. 

The two main features of the theory which may most easily be checked are 
(i) the energy spectrum of evaporated particles, (ii) the ratio of doubly to singly 
charged particles. An energy analysis has been made of a group of 17 stars 
from 7 to 17 prongs in the C2 plates. The proton energy spectrum as deduced 
from this analysis combined with the results of Perkins have been shown by 
Le Couteur (1950) to be in good agreement with his theory, which also accounts 
for the a-particle spectrum. The experimental determination of the ratio of 
doubly to singly charged particles will now be described. 


(ii) The Ratio of Doubly to Singly Charged Particles 
The proportion of doubly charged particles emitted was found for all stars 
in the Ilford C2 plates and the results are shown in Figure 3. These results 
appear to be confirmed by those of Addario and Tamburino (1949). They differ 
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somewhat from the latest (unpublished) results of Perkins although the differences 
are not as much as three times the standard deviations. 

The proportion of doubly charged particles decreases with increasing star 
size, remaining more or less constant for stars of above 6 prongs. ‘There are 
392 stars with more than 6 prongs with an average size of 10 prongs and the 
value of the ratio is 0-271 + 0-008. This is in excellent agreement with the values 
obtained by Le Couteur, who found an average value of 0-27. Le Couteur 
predicts a variation of from 0-25 for 2-pronged stars in heavy elements, to 0-28 
for large stars. 

These figures may be used, with the experimental results, to give a value for K, 
the proportion of doubly charged particles from the disintegration of carbon, 
nitrogen or oxygen, if K is assumed constant for stars of from 2 to 6 prongs. 


Mean Proportion of « Particles per Star 


0 2 4 6 8 10 12 14 
Total Number of Prongs per Star 


Figure 3. 


This value may then be compared with that obtained by Perkins (1949) from a 
study of stars in a gelatine sandwich plate. 

Using Le Couteur’s figure for stars in silver and bromine the relative number 
of stars formed in the light elements of the emulsion may be calculated in terms 
of K. If the relative number of stars formed in the different elements of the 
emulsion is proportional to the geometrical cross section of the nuclei of the 
elements, approximately 27°%, of the total number of stars in the emulsion will 
be formed inlight nuclei. Harding (1949) has recently made a direct determination 
of this quantity by comparing the stars formed in gelatine and in ordinary emulsion. 
He finds approximately 36 +4-5°% of all stars are formed in the light nuclei, a 
result in reasonable agreement with the above value. 

A figure of 27% for the proportion of stars formed in carbon, nitrogen and 
oxygen is found, if K is assumed to be 0-54. Perkins (1949) from a study of stars 
in gelatine found for K a value of 0-5. 

On the Le Couteur theory, approximately 28% of all singly charged particles 
emitted from stars should be deuterons or tritons. From the relative numbers 
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of protons, deuterons and tritons from stars which ended in the emulsion it was 
estimated that approximately 25°, of all emitted singly charged particles were 
heavier than protons, and Harding et al. (1949) give a value of 30%. This may 
be taken as a further confirmation of the theory. 


§5. THE HIGH ENERGY PARTICLES 


The 312 stars in the electron-sensitive emulsion were analysed for high energy 
particles in the following way. The particles were divided according to their 
energies, those with energy corresponding to an ionization greater than 7 
times minimum, those between 7 and 3-5 times minimum, and those of less 
than 3:5 times minimum. For protons these ionization values correspond to 
energies of 40 Mey. and 100 Mev. The numbers of particles in the three groups 
are shown in Table 2. 


Table 2. Relative Number of Particles with given Value of Ionization 


Ionization (xX minimum) S7/ 5 =7/ <a3i5 

Stars with 8 or less evaporated particles 1 0-103 +0:010 0-156+0:013 

Stars with more than 8 evaporated particles 1 0-079 +0-026 0:380-40-063 
1 


All stars 0100 -+.0:010 0-:180-+0-013 


This grouping was chosen to give statistically significant results, and because 
the proportion of high energy particles was observed to increase for stars with 
more than eight evaporated particles. Charged primary particles have not been 
included. 

It was found that 56 + 10% of the particles with ionization less than 3-5 times 
minimum emitted from stars with /ess than nine prongs had ionization below 
1-5 times minimum. ‘This is in agreement with the value of 56 + 8°% calculated 
from the results of Brown et al. (1949 b) for stars with more than seven evaporated 
particles. 

It can be seen from Table 2, that whilst the proportion of particles with 
jonization between 3-5 and 7 times minimum does not vary within the limits of 
error, there is an increase in the proportion of those with ionization less than 
3-5 times minimum for large stars. Now the ratio of the number of directly 
ejected protons to the number of charged evaporated particles should not vary 
appreciably with the energy involved in the disintegration (Rosenfeld 1948). 
This experimental increase may therefore be due to the creation of mesons, 
since the probability of this occurring will be greater in high energy disintegrations. 

In the absence of scattering measurements it is dificult to estimate the 
proportion of tracks with ionization less than 7 times minimum which may 
be ascribed to created mesons and the proportion which may be ascribed to recoil 
protons. A small number of mesons with energies of a few Mev. are certainly 
emitted, and it is known from penetrating showers that very high energy mesons 
are also created. 

It is of interest to note that the proportion of particles with ionization lower 
than 3-5 times minimum which have an ionization less than 1:5 times minimum 
is the same (56%) for small and large stars. If, as suggested above, a large pro- 
portion of particles in the group may be ascribed to mesons in the case of large 
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stars, this indicates that some 60°% of such mesons have ionizations less than 
1-5 times minimum (or energies in excess of 60 Mev. if 7 mesons) and about 40% 
are slow mesons (i.e. if 7 mesons, of energies less than 60 MevV.). 

Further, the fact that the proportion of particles with ionization between 
3-5 and 7 times minimum does not appear to vary much with star size may indicate 
that only a small proportion of these are due to mesons, and the majority are 
recoil protons of energy 40-100 Mev. 


ACKNOWLEDGMENTS 


I wish to thank Professor P. M. S. Blackett for the facilities to do this work 
and for his interest and advice. I also wish to express my thanks to Dr. G. D. 
Rochester for valuable discussion, help and advice, and to Dr. Le Couteur for 
discussion on the theoretical aspects involved. Miss M. Blundell has developed 
and scanned many of the plates and assisted at all stages of the work, and I wish 
to express my thanks to her and to all those who assisted in scanning the plates. 
I am indebted to Dr. G. P. S. QOcchialini of Brussels for exposing and 
developing the plate at 2,800 m., to the Staff of the Scientific Station, Jungfraujoch, 
for their co-operation in exposing the plates at 3,450 m., and to the staff of the 
R.A.F. at Benson, Oxon., for making the exposures at 12,000 m. I should also like 
to thank the staffs of Messrs. Ilford and Messrs. Kodak, who have supplied us 
with the plates and advised us on their processing. 


REFERENCES 


Avpario, M. M., and TamBurRINoO, S., 1949, Reports of the Como Conference (in the press). 

BaccgE, E., 1938, Ann. Phys., Lpz., 33, 389; 1946, Cosmic Radiation (edited by Heisenberg) 
(New York: Dover Publications), p. 128. 

BERNARDINI, G., CorRTINI, G., and MANFREDINI, A., 1949, Cosmic Radiation, Colston Papers, 
Duzer 

BETHE, H. A., 1937, Rev. Mod. Phys., 9, 69, 245. 

Bour, N., and Katckar, F., 1937, Proc. Cop., 14, No. 10. 

Brown, R. H., Camerini, U., Fow er, P. H., Murruead, H., Powe t, C. F., and Ritson, 
D. H., 1949 a, Nature, Lond., 163, 47. 

Brown, R. H., Camerint, U., Fowier, P. H., HEITLER, H., K1nc, D. T., and Powe t, C. F.,. 
1949 b, Phil. Mag., 40, 862. 

Cocconl, G., 1949, Phys. Rev., 75, 1074. 

CorTINI, G., and MANFREDINI, A., 1949, Nature, Lond., 163, 991. 

DitwortH, C. C., OccHIALINI, G. P. S., and Payne, R., 1948, Nature, Lond., 162, 103. 

GeorcE, E. P., and Jason, A. C., 1949, Cosmic Radiation, Colston Papers, p. 31. 

HarDInNG, J. B., 1949, Nature, Lond., 163, 440. 

Harpinc, J. B., Lattimore, S., and Perkins, D. H., 1949, Proc. Roy. Soc. A, 196, 325. 

Horninc, W., and BAuMHOoFF, L., 1949, Phys. Rev., 75, 370. 

Le CouTevr, K. J., 1950, Proc. Phys. Soc. A, 63, 259. - 

Perkins, D. H., 1949, Phil. Mag., 40, 601. 

ROSENFELD, L., 1948, Nuclear Forces (Amsterdam : North Holland Publishing Co.), p. 266. 

WEIssKoPF, J., 1937, Phys. Rew., 52, 295. 


259) 


The Evaporation Theory of Nuclear Disintegrations 
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ABSTRACT. The relationship between nuclear temperature and excitation energy is 
discussed taking into account the thermal expansion of the nucleus and the effect of the 
neutron excess. Using these results, the evaporation theory is put into a form suitable for 
the treatment of highly excited nuclei. The competition between the emission of proton, 
neutron, deuteron, triton and a-particle is considered in detail, using Weisskopf’s method 
of detailed balance. 

Detailed numerical values are calculated for the disintegrations of silver bromine 
observed in photographic plates. The main points of comparison with experiment are 
(1) the energy distribution of the evaporated particles, and (ii) the ratio of emission of singly 
and doubly charged particles. The comparison yields close information about the relation- 
ship of nuclear energy and temperature; the agreement with theory is as good as can be 
expected. It is necessary to assume that the Coulomb potential barrier is reduced at high 
excitation energies. 


Sil, ONIN ROBO AD KON) 
UCLEAR disintegrations of high energy, initiated by cosmic rays, have 
been studied extensively by the photographic emulsion technique. ‘The 
recent papers of Harding, Lattimore and Perkins (1949), referred to as 
H.L.P., and Page (1950) contain sufficient statistical evidence to provide 
an effective test of some points of nuclear theory. 

The early work on disintegration stars was reviewed by Bagge (1943). From 
this and subsequent researches (H.L.P.) it has been recognized that the emitted 
particles can be divided roughly into two groups, namely (i) a few particles of 
high energy, and (ii) a greater number of slow particles with an energy spectrum 
typical of nuclear evaporation. Of course there is no clear-cut distinction 
between the groups; if the energy spectrum is plotted (i) appears as a high energy 
tail added to (ii). Qualitatively it is clear that the fast particles must consist of 
nucleons which have come into direct collision with the primary particle, often 
a fast neutron, initiating the star, and of mesons created as the Bremsstrahlung of 
these collisions. The theory of such processes is not yet in a satisfactory state, 
for example Zwanikken (1948) shows that the static nuclear forces cannot account 
for the large energy transfer from the primary particle to the nucleus: presumably 
the difference must be attributed to velocity dependent forces of which little is 
known. Before finally escaping, these fast particles must lose energy to the 
nucleus and so leave it in a highly excited state. Further particles are emitted 
until the final nucleus is stable; it is these secondary particles which one attempts 
to treat by nuclear evaporation theory. 

The statistical treatment of nuclear disintegrations, initiated by Bohr and 
Kalckar (1937) and amplified by Bethe (1937) and Weisskopf (1937), can be 
applied to the secondary process (Bagge 1938, H.L.P.) and leads to definite 
predictions. But the excitation energies involved in large stars are very much 
higher than those envisaged by the theory, which assumes that the excited nucleus 


* Now at University of Liverpool. 


260 K. F. Le Couteur 


is always in approximate thermodynamic equilibrium. ‘This requires that the 
interval between the emission of successive particles is longer than the time 
required for the convection of energy across the nucleus. This condition is 
probably satisfied since, for all but the largest stars, the mean kinetic energy of 
the evaporated particles (perhaps 13 Mev.) is much less than the kinetic energy 
(22 Mev.) of the nucleons at the top of the zero-point Fermi distribution, 
which is mainly responsible for the convection of energy. ‘There is also the 
requirement that the excitation energy should be much less than the total binding 
-energy of the nucleus. For this reason, dealing with the stars created in silver 
and bromine with an average mass number of 100, it has not been thought worth 
while to extend the calculations to excitation energies above 600 Mey. ‘Though 
this figure is a high one, it must be remembered that the nucleus cools down as 
particles are emitted so that most of the particles are emitted when the excitation 
-energy is much lower and under conditions to which evaporation theory may be 
applied. Because of these limitations, one should not apply evaporation theory 
to stars with more than about 14 prongs. On the other hand, in photographic 
emulsions, stars with less than 7 prongs may originate in the light elements and so 
-cannot be compared with the theory. However there remains a useful range of 
star sizes and the greater part of the experimental evidence actually falls within 
this range. H.L.P. have verified that, for such stars, the particles are emitted 
with random direction as would be expected for an evaporation process. 

Much of the experimental evidence to which reference has been made was 
derived from photographic plates which, for various reasons, were insensitive to 
fast particles. Energies above 40 Mev. were not measured by Page. ‘Thus the 
available data chiefly concern the secondary particles which are best treated by 
evaporation theory. ‘The main features which can be compared with experiment 
are (1) the energy distribution of the evaporated particles, (ii) the relative numbers 
of particles of charge e and charge 2e. Much more experimental evidence exists 
for (11) than for (1) because different charges can be distinguished at sight whereas 
energy measurements are difficult. By combining both types of evidence one 
obtains fairly precise information about some of the basic parameters of nuclear 
disintegration. 

The most important relationship of nuclear thermodynamics is that connecting 
the temperature 7 and excitation energy U. The ‘Fermi gas’ model, corrected 
by Bardeen’s (1937) method, has been used and is discussed in §2 where it is 
extended to take account of the thermal expansion of the nucleus and of the neutron 
excess. It is reasonable to use something resembling a free particle model at 
these high excitation energies. The alternative ‘liquid drop’ model (Bethe 
1937) is not very suitable for high excitation as it gives U =ar73 + br4 and, since 
the number of degrees of freedom is only 3A this formula must be restricted to 
low values of 7. It could be amended by Debye’s method and the ‘critical 
temperature’ is found to be about 15 Mev. ; as appreciable corrections are required 
at less than half the critical energy the method is not very attractive. The ‘ liquid 
‘drop’ formula, without emendation, was tested-by H.L.P. and found to be 
unsatisfactory. 

‘The theory has been used to calculate the expected numbers of evaporated 
neutrons, protons, deuterons and 3H, 3He and 4He nuclei. Heavier particles 
could be considered but the emission probabilities are quite small. Weisskopf’s 
(1937) method is to derive the differential emission probabilities by consideration 
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of detailed balance in a state of thermodynamic equilibrium. With the assump-- 
tion that, apart from the well-known dependence upon energy and upon the 
Coulomb barrier, all these particles have the same cross section for capture by 
the residual nucleus, it is not necessary to introduce explicit factors to represent 
the ‘probability of existence’ of compound particles within the nucleus. It 
would not be possible to obtain agreement with the observed ratio of singly and 
doubly charged particles if further small factors for the existence of «-particle 
and triton were introduced into the formulae. 


2 ee ACS hie Aly NOD ORTEXCULEDENUGLE I 
2.1 The Hartree Approximation 


The thermodynamics of an excited nucleus is completely determined if the 
total energy is known as a function of the temperature and volume. From this 
information, together with the boundary condition of zero external pressure, one 
can for example determine the equilibrium volume or radius at any temperature. 
In principle it should be possible to derive these quantities directly from the 
nuclear forces, but in practice, even for the ground state of the nucleus, the 
calculations are not very reliable. 

While the statistical mechanics of a number of independent particles in a 
fixed field of force is well known, the mechanics of particles in interaction hardly 
exists. ‘Therefore to proceed with the problem it is necessary to adopt Hartree’s 
approximation and assume that the force between nucleons can be represented 
by an average potential field. ‘Then the different nucleons are supposed to occupy 
individual quantum states in the field and, at any nuclear temperature 7 (measured 
in energy units), the probable number n, of particles in the sth state is given by 
Fermi’s formula 

1 
Ms~ exp {(e,—)/t} 41° 
Pfs J 

where ¢, is the energy of the state and <’ is the Fermi energy. Use of this result 
involves the assumption that the Hartree field is unchanged by small changes in 
the occupation numbers 7,._ In the nuclear problem there is some justification 
for this step since only very low temperatures are considered and so only a few 
particles at the top of the zero-point Fermi distribution can be excited ; this should 
not greatly change the mean potential field which is chiefly due to the particles in 
the lower part of the distribution. This use of a Hartree approximation leads to 
a close resemblance between the calculations for excited nuclei and the well-known 
calculation for the ground state of the nucleus (Bethe and Bacher 1936, 
Rosenfeld 1948). 

It is necessary to sum various functions over the distribution (1) ; by introducing 
the density in energy g(c) of the individual quantum states their sums may be 
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where «_ is the energy of the lowest individual state. The integral can be 
evaluated by Sommerfeld’s method by splitting the range of integration at «’ 
and the result is 
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’ is itself dependent on 7; if €9 is the value of ¢’ at zero temperature one has, 
neglecting terms proportional to + and higher powers of 7, 


Eo 2a d d 
OM= | gg de +(c' —€)(8)e, + (1 = +2 2) ee (4) 


The value of ¢’ — ¢, is fixed by the condition that, since the total number of particles 
is fixed, the special choice g=1 must lead to Q(r)=Q(0). Therefore (4) reduces 
to 


OG) = Oo + a Lo (4) aes (5) 


with Q,=Q(0) etc. It also follows that <’—€, is of order 7*. In particular by 
putting g=e one obtains the total energy as a function of temperature, 


qq? 


E(n)=Ey+ gy ee (6) 


In this approximation the dependence of E(7) on the volume of the nucleus appears 
as a dependence on the extensive parameter, x say, which fixes the linear scale of 
the Hartree field; it is sometimes convenient to denote this explicitly by writing 
E as E(r,x). Then from (6), one obtains the free energy F(z, x) as 
pie 
F(z, x) = E(x) = 6m: 2x): Santac (7) 

In the absence of external pressure, at temperature 7 the equilibrium value x, 
of x must satisfy the condition 

dF (1, x) dE (t,x) 7? dg(x) _ 0 9 

dn = de. ae ae Oo 


and so, if x9 is the value of x, for the ground state, 


ee (cors 
7 TO Nat eae) 
This is the law of linear expansion for the nucleus. 

It remains to calculate g(x) and E,(x) which requires a more specific approxi- 
mation. With plane wave functions and a Gaussian potential, g,(x) has been 
calculated by Bardeen (1937) and (d?K,/dx?) by Bethe (1937, equation 318), 
(d°E,/dx") determines the compressibility of nuclear matter and, according to 
Feenberg (1947), can to some extent be checked by an examination of nuclear 
regularities. It is necessary to extend Bardeen’s work a little to take account of 
the volume expansion and to apply it to nuclei with unequal numbers N and Z of 
neutrons and protons. It is found that the expression for go(x) can be put into 
a more simple form by a method equivalent to use of the virial theorem. 


2) 


2.2 Plane Wave Approximation 


The individual wave functions are approximated by plane waves periodic in a 
box of volume v=47R?/3 equal to that of the nucleus. In terms of the wave 
number k the density of individual states of proton or neutron is 

Cawak 
ge) == R? Gece 0S ee (10) 
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With N and Z unequal it is necessary to replace ¢’ by two parameters e, and «, 
given by 


2 2 
“slat-ann 7"alapeeayer oo 
and so at zero temperature all states are occupied which have wave numbers below 
Rao = Rott and | ah ae (12) 
where n=2N/A, p=2Z/A and N+Z=A=2vk,3/3n%. ...... (13) 


In the ground state the total potential energy of the nuclear forces (strictly the 
volume contribution only), 


Vena ot ee | (14) 
is determined from certain ‘mixed densities’ of the form 
1 sin kr —kr cos kr 
pti.ts) = or exp {¢k.(r,—r,)} = (“oa Be (15) 


with r= |r,—r,|, by equations («, 8=n or p) 


Vig= | pelt Po)Jae(Ti ¥2)Pp(Fir Te) 47 d°*7o, wees (16) 


In this equation J,,, represents the mean proton—neutron interaction after the 
spin summation has been carried out, J,, and J,, have similar interpretations. 
‘Thus the expectation value of the total energy ¢,(A) of a proton of wave number kis 
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For a Gaussian potential 
J n= —Bexp(—7/a*), ee CX F510") a uae (20) 
the last equation, after multiplication by v/37?, takes the simple form 
de de h? 2 ( OV ) 
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Bethe and Bacher express V in terms of a parameter «=ak« a/R and then 
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Since the range of the force is fixed, « is an extensive parameter for the nucleus. 
Therefore, if T' is the total kinetic energy, the value of 0V/dx? for the ground 
state follows at once from 


PS V0 tor (23) 


3 he 3n2A 
Now T = = 5g (Nia + Zh) = Fo ag (P+ PH) Roe, wane (24) 


which gives for x=%p, 
Date OV eer L eC mane ae ay 


34°” O28 3AOk? 5 M\ 2 
; dé\ O01 ee eT 
Thus (aE). ag lt ete te (25) 


With n=p=1, this reduces to 
yee Lot axeey 2h) ae ee eee (26) 


which, because of (23), is equivalent to Bardeen’s more complicated expression. 
By inserting (25) into (10), the basic result (6) for E(7) can be written explicitly as 
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where €=(h?k,)/2M) is the kinetic energy corresponding to ko, empirically 


€=22mev. With y,=y,=1, equation (27) reduces to Bethe’s formula for the 
free particle model. 


2.3 The Coefficient of Expansion 


To determine the expansion from (9) one needs dg,/dx which must be 
calculated from (22) as it is not permissible to differentiate (25). It is sufficient 
to consider only the simple casen=p=1. Then 
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dg,|dx 2 6 
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There (9) gives the law of éxpansion as 
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For convenience this result is stated as an expansion of the mean nuclear 
radius but, since it is derived from the simple statistical model, it must be undez- 
stood as giving the change of mean density of the nucleus with increase of 
temperature. 


2.4 The Excitation Energy 


The excitation energy U(r) required to raise the nucleus to temperature = 
is defined as 


Gna (Oey), oa Mea (30) 
Now, according to (8), 
Prim l(a Ola hae ~ soc) ina. eee (31) 
qq? 
Thus E(z, x,) = EO, Xo) + = (280(%,) — So(%o)) + O(7*) 
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since £o(x,) —£o(Xp) is of the order 7?. Therefore, from (30) and (27), 
24 nil pe 
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Since n—1=(N—Z/A)=0 is small, it is convenient to expand the right-hand 
side in powers of 6. From (25) and (26), 
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Also € and y depend on x,” which is itself dependent on 6. According to Bethe 
and Bacher (1936) *, the nuclear radius is given by 


(DEERE) PE en ae ett rs (35) 


so that, if X, is the value of x) corresponding to @=0 and (,y the corresponding 


Cy; 
y= X9(1 — 0-5 62) 
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* The numerical coefficient 0-5 in (35) is less than the coefficient 0-66 given by Bethe and 
Bacher; the smaller value is obtained by taking account of the forces between like particles. 
According to Bethe and Bacher, the interaction parameters must satisfy B+ }C=70 mev., also 
Rosenfeld (1948, § 11.11, Table 2) states that in a charge-independent theory B and C must be equal. 
Thus we take B=C=47 mev. These parameters lead to the value B’=22 Mev. for the coefficient 
(equation (43) below) which expresses the dependence of binding energy on isotopic number. 
This agrees better with the empirical value, 20-5 Mev., than does that derived by Bethe and Bacher 
with weaker forces between like particles. Since the parameter y, introduced in (38), is very 
closely related to f’, it can be assumed that y also is given fairly accurately by the theory. : 
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Therefore the final result is, from (33), (34), (36), 


(7,0) (sO) = ee (37) 

pea io eae 38 
with y= le ey 5 (14 Se ar =2 ee (38) 
and Ment Aly. > a= = a ee (39) 


The entropy S, of the nucleus is given by 


S(z, 0) = = \\) Ue (40) 


and, according to Bethe (1937, equation 281), the density p of the nuclear levels 
and the average spacing D between levels are given by 


1/2 
D(U)= aes =1(2n So) eS 
dr 


A is the most important parameter of nuclear thermodynamics. 


2.5 Numerical Values and Discussion 


The nuclear radius R is taken as R=7,A1* with 7,=1-47 x 10-4 cm.; this 
leads to €=22 Mev. 

Rosenfeld (1948, p.153) quotes a=1:9x10-%cm. for the range of the 
Gaussian potential as derived by analysis of two nucleon systems. This gives 
x= 1-5a/ro—1-93-andiy = 75. 

The present work has been carried out with the Hartree model for large 
nuclei and it may be more consistent to use force parameters chosen so that this 
model yields the correct binding energy. Bethe and Bacher, using the first 
approximation only, arrive at x=2-7 which leads to y=2:46. ‘This is probably 
an over-estimate since it is known that the method requires very strong force 
parameters; thus all that can be said is that y should lie in the range 1-75 to 2-46. 

The results of the theory have been expressed in terms of the ratio of energies 
(7/Cy) and one must regard Cy as the characteristic energy for thermal effects. 
Now the calculations deal only with the volume energy and could not easily be 
extended to include the surface energy in a rigorous way. But the latter is not 
negligible, for A=80 the volume energy is —14-46 and the surface energy 
+3-58 Mev. per nucleon. As a rough correction it seems reasonable to suppose 
that the surface contribution varies with energy in the same way as the volume 
contribution; from this point of view the characteristic energy Cy should be 
reduced in the ratio (total energy)/(volume energy), which for A =80 gives a 
factor 0:75. Leaving ¢ unchanged this is equivalent to a reduction of y from 
2°46 to 1:35 or trom: 1°75; torl-3 1° 

We shall see that a value of y of about 1:5, giving ¢y =33, is quite consistent 
with the data on disintegration stars and it alee agrees reasonably well with the 
older data on nuclear reactions at low energy. 
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‘The degeneracy parameter for the Fermi statistics is, in effect, 7/Ly. Now, 
even for an excitation energy as high as 400 mev., with A = 100 and y = 1-5, equation 
(37) gives 77 Mev. so that 7/fy = 0-22. Stoner (1936) has shown that for such 
small values of the parameter the error in the energy due to the neglect of terms 
in (7/¢y)* etc. is negligible. 

With Bethe’s values, x=2-7 and D=9-6 equation (28) gives ju/xp=4-15 and 
the law of expansion becomes 


Beets 025 Gye. 2. rs | eee. (42) 


‘Thus in the above example with 7/fy 0-2, the expansion of the mean radius is 
about 10% corresponding to a reduction in mean density of 30°%. Equation (38) 
now gives y=1-3; with this value of y the dependence of S(z7,6) on @ is almost 
negligible in all applications. 

So far nothing has been said about the Coulomb forces. For a given nuclear 
volume the total Coulomb energy is obviously independent of temperature, 
therefore equation (27) and the subsequent equations for the excitation energy 
remain true when the Coulomb interaction is taken into account. Strictly 
speaking, the Coulomb term should be included in the expression for d?E/dx? 
used to calculate the expansion, but the effect is negligible. 


{37 NUCLEAR EVAPORATION 
3.1 The Energy Balance 


It is necesssary to use the semi-empirical representation of the binding 
energy of a nucleus in its ground state 


— 2 
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For values of A in the range 70-100 with which we are chiefly concerned, it is 
best to take »=8-6Mev. The empirical value of 6’ is 20-5 Mev. which gives 
B=23 Mev. Also the value of 6 for the most stable nucleus N,, Z, of mass A is 
well represented by the empirical formula 
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Suppose that a nucleus a, with excitation energy U, emits a fragment (7, 2) 
leaving a residual nucleus } with excitation energy U,. If the kinetic energy of 
the fragment is J and its binding energy is w, the energy relation is 


U, + E,(a) =U, + 4d) -w+T 
or U,=U,— OG b)—- Tee eee (49) 
where O(a, 6) =a(n +2) —2(6—0,)(mBy—2Bp)—W soe (50) 


is the threshold energy for emission from the nucleus of a fragment (7, 2) with 
zero kinetic energy. 


3.2 The Emission Probabilities 


Weisskopf (1937) has derived a formula for the probability per unit time of 

emission of a particle with kinetic energy T from a nucleus a leaving a residual 
nucleus J, 
o Tp(>) 
% p(a) 
where y=opgm/7?h® and p is the level density. Here g is the number of spin 
states and m the mass of the particle. a, is the geometrical cross section of the 
residual nucleus and a is the total cross section for capture of the particle by b. 
For neutral particles is it usual to assume o = oy but for positively charged particles 
the capture process is opposed by the Coulomb interaction; if V is the height of 
the potential barrier the classical theory gives 


lalog=I—V, 


which in conjunction with (51) implies that particles of kinetic energy less than V 
cannot be emitted from the excited nucleus. Quantum mechanics leads to a 


larger cross section because charged particles can to some extent penetrate the 
potential barrier, thus one gets 


P(T) dT =y dla = =e (51) 
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The theory has been worked out and applied to nuclear reactions by Bethe (1937) 
and Bethe and Konopinski (1938) and has been invoked by H.L.P. to account 
for the emission of low-energy charged particles from cosmic-ray stars. With. 
this expression for o, the emission probability becomes, after use of (41), 


P(T) dT =y(T—V’)(U;{U,) 4 exp(S,—S,)dT. ...... (54) 


To determine the total emission probability and mean energy of emission it 
is necessary to evaluate 


V’ 
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with X=Yo—J. The integration may be carried out after expansion of the 
binomial. Since A =5 terms in exp(—Ayy) may be neglected for R>1 Mev. 
Thus, the emission probability per unit time is 


P(T)dT={"* P(T)aT =4yU" a : 
Me ( ) PACs exp(— 5S, +tAVR) ay 1+ gop t+): 
ee (56) 
In this equation A must be calculated for the residual nucleus 6. Similarly 
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Thus the mean energy T of emission is given by 
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The first term in this expression corresponds to the simple formula T=V' +2r, 
the correction term 6/A? is relatively small. The terms involving negative 
powers of y) or R in (56) and (58) are negligible for R>1 Mev, ; since in any case 
the integrals behave in a perfectly regular way as y)—>0 it is best to omit these 
terms altogether. 


3.3. Competition between Emission Processes 


The equations of §3.2 can be used to determine the probability of emission of 
various types of particle from an excited nucleus. ‘There is competition between 
the different possibilities ; itis sufficient to consider only the emission of neutrons, 
protons, deuterons and 3H, #He, He nuclei which are respectively denoted by 
n, p, d, t; t, «. 

According to (51) the relative weight factors y are 


Ve = Var eH Vem= ss © Ya Pls ee (59) 


and it will be seen that the factor 3 leads to an appreciable probability for triton 
emission. 

The formulae of § (3.2) involve the parameter R, which represents the excitation 
energy of the residual nucleus left after the excited nucleus a has emitted a particle 
x(n, p,....«) with kinetic energy V,’. It is useful to collect together the values 
R, of R=U-—Q-V' as 


R,= U,,—8:6 + 1-68(8 — 0) R,=U,,—174+0:88(0-8,)-Vy, °] 
R,=U,—-86—2-48(0-0,)—V,, - Ry=U,—18-1-3-28(0-8.)— Ve, | 
Ue 15=0:3p0—0.)—V; eae 2 = lone =n =, 
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Also if the initial nucleus a is A, 6 the A values for the various residual nuclei are, 
with A=A(A —1, 9), 


A, =A(1 —y0/A) A, =A(1—1/A —y@/A), | 
A, =A(1 +y6/A) A; SAQ= AS OA eee (61) 
Ne= A=) he 32) 


Now from (56), the relative probabilities of neutron and proton emission from the 
excited nucleus are 
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The factor (R,/R,)"4 in (62) may be expressed as exp(R, —R,/4R,,) and, since 
A =5, is always negligible in comparison with the last factor. ‘The terms which 
appear in (63) may be taken from (60) and (61) and so the relative emission 
probabilities take the form 
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The interpretation of the terms in these emission ratios is of interest as some of 


them represent strange effects which are important only at high energies. Thus 
equation (65 a) shows at once that 


(i) the potential barrier V favours neutron emission, 


(ii) the term 26(0—06,)=92(N—N.)/A favours neutron or proton emission 
according as N—WN, is positive or negative. It thus acts as a ‘ governor’ 
controlling the emission in such a way that N — N, never becomes very large. 
This term arises because the energy available for emission of particles is 
greatest if the isotopic number of the residual nucleus is that of maximum 


stability. In stars with several evaporated particles this effect requires that 
more neutrons than protons are emitted, 
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(iii) the term —2y@R/A favours proton emission if 6 is positive ; thus it endeavours 
to increase the value of | N— Z|, reflecting the fact that the density of levels 
of the residual nucleus increases with | N—Z]|. 


The order of magnitude of $V,’ is 1 to 3mev., depending on assumptions 
about the temperature dependence of the potential barrier; so for a nucleus of 
A=80, the first and second terms become comparable when N — JN, lies between 
land 3. ‘The third term is relatively small; with y=1-3, 0=0-:13 and A=80 it 
becomes 0-4R/100 and so becomes comparable with the first two only when 


R,.>300 Mev. 


np 
The relative probability of « emission contains the term 


(3/24 +y6/A)R ~1-67R/A 


which is important because it increases with excitation energy. The contribution 
3R/2A arises from the dependence of the level density of the residual nucleus 
upon its mass number, which discourages the emission of heavy particles. For 
excitations above 100 Mey., this is the most important single term determining 
P,/P,. Itis obvious from the form of (65 e) that P,/P,, has a rather flat maximum, 
which occurs at R~100Mey.; this has the important consequence that over a 
considerable range of star sizes the ratio of « to proto 1 emission from silver or 
bromine remains almost constant. In the absence of t e term 3R/2A, the ratio 
P,,/P, would increase steadily with energy ultimately approaching the limit 2. 
Thus the ratio of « to proton emission would increase seadily with star size, very 
much in disagreement with experiment. 


3.4 Integration of the Equations of Emission 


The emission process is treated as a continuous one. ‘The mean reduction 
in excitation due to emission of a particle of type x(x =n, p....) is given by (49) 
and (58) as 
He Or VAAN A Ro ON. stn (66) 


: : . . : . 7y7\—1 
Thus the mean number of particles emitted per unit drop of excitation is (H) 


where 


and the mean number of particles of type x emitted per unit drop of excitation . 


pe ee rr pent ei (68) 


ee Ep eee 

Of course the right-hand side of this equation depends on what particles have 
already been emitted. We have here a set of six simultaneous integral equations 
governing the emission of the six kinds of particle considered; strictly speaking 
one should integrate these numerically, working downwards from the initial 
excitation energy, a separate integration being required for each value of the 
initial excitation. However the equations are coupled mainly through the Y 
terms in (65) and it is possible to arrive at an approximate solution for 6 which is 
not very sensitive to the initial excitation energy. Then one has only to integrate 
the six independent emission equations rather than a set of six simultaneous 
equations. 
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3.5 The Equation for the Neutron Excess 
The rate of change of 0=(N—Z)/A due to the emission of particles is 

do 1 6 dA 

Frag eae ome eye) eh or 


Let 7 be the total number of particles (n, p,d....) emitted by the nucleus up to 
the time when its excitation energy has fallen to U, then dr/dt =X P, and 


eee (70) 
dr A Pape A dr 
U igs 1 Py Py es els} Py 
mow PD z log Pi Voor 2 Poe ea tee a) 
are given by (65) and may be treated as equal. ‘Then, with the notation, 
@S(P ohh Fhe Pie =! Saha) ee (71) 
one finds 
Sree! d 
H =, (8-8) = — 7 (8-6) 
aN eee 6 6 d0,.\dA 
= SAR 28-841, =2y 5 Ry +(5 +54 dr . 
By writing 
hi 2A/ 0%) d0\tidA 
= —= —|(—+——)-—=-5), ...... 72 
x=4/R,,, ; pt+qx and 7 (a+94)- ie b, (72) 
this becomes 
d 2BA 
(p+ 9x) © 28(0-0,) — “* (28(0-0,)-G(a)}=0, (73) 
where G(x) = be + 2900") Ae Se eee (74) 


In equation (73) one can treat p,q, A/A, 0,/A as constants. Then, neglecting 
derivatives of G beyond the second, the solution is, with  =2BA/Aq, 


2B(0—6,)=J(x)+K(ptqxy, see (75) 
COE ee 
q(e—1) dx * g(u—1)(u—2) dx? 


and K is a parameter determined by the initial conditions. 

In this form the solution is sufficiently general to take account of the increase 
of 6 with excitation energy and also of the possible decrease of the potential 
barrier V’ with increasing energy. According to Bagge (1938), the potential 
barrier should fall off approximately linearly with increase of energy, that is 
VeVi exe 

With any given assumption about V and A one can estimate a, b, p,q, numer- 
ically by use of equations (65); then the solutions (76) can be worked out explicitly. 
Figure 1 shows results obtained in this way using various parameters (detailed 
in §4.1 below) and for several different values of the initial excitation energy. 
The coefficient q is rather larger than the value 4/A\ suggested by (66) because as 


where J (x) = G(x) + 
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the energy is increased the emission of heavy particles becomes more probable, 
therefore H increases faster than H,, or H, and a decreases with energy; both 
factors lead to a more rapid increase of Hla with energy. 

The variation of 28(@—0@,)~N-—VN, as the nucleus cools down is easily 
understood. To maintain the neutron—proton balance of the residual nucleus 
roughly 1-5 times as many neutrons as protons must be emitted during the 


Parameter (a) Parameter (bd) 


MeV. 


MeV. 


Figure 1. Variation of the neutron excess during evaporation. 2B(8—Oc) ~ N—Ne as a function 
of x = V/U for three different initial values of the excitation energy U. Also the auxiliary 


functions Vp’, G(x) and p+qx of §3.5. 


evaporation. But, starting with 9=6, at an initial excitation energy of say 
400 mev. (x20), neutron and proton emission are at first about equally probable, 
therefore the neutron excess of the residual nucleus increases at first rapidly and 
then more slowly as the ‘controlling term’ comes into play; later when the 
excitation has fallen below about 100 Mev. the potential barrier leads to a further 
increase of neutron emission so that eventually, at x =0, the residual nucleus has a 
deficit of between 1 and 2 neutrons. It will be noticed that this final value is almost 


independent of the initial excitation energy. 
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3.6 The Emission of Low-Energy Particles 


The results of the previous section for the term 28(9—@,) have an important 
physical consequence. Figure 1(b) gives —1-7Mev. for the value of 26(6 —8,) 
at x=(. By inserting this value into the expression (50) or (60) for Q, one finds 
that the excitation energy necessary to release the last particle with zero kinetic 


energy is Q,, given by 
OQ, = 10 mev., Q, = 6:6 Mev., OZ=49 Mean Fae (77) 


Thus when the excitation energy has fallen below 10 Mey. neutrons cannot be 
emitted but charged particles can, even though they must penetrate the potential 
barrier. There is however competition from two other processes, namely 
y-ray emission and fB-decay by K capture. Estimates of the radiation width 
for y emission are a little uncertain, Bethe (1937, equation 729 b) suggests a width 
of approximately volt which represents a lifetime of 10-!°sec., very much less 
than the time required for a 8 process. Protons and «-particles can be emitted 
in less than 10-sec. only if their ‘penetration probabilities’ for the potential 
barrier are greater than 10~!° and 10-® respectively. With A =66, Z=29, which 
are reasonable values for the residual nucleus of a star of about 12 prongs, Bethe’s 
table XXXVII shows that these probabilities are exceeded if the proton and 
a-particle are emitted with kinetic energies greater than 0-5 Mev. and 4Mev. 
respectively. Therefore the threshold energies (77) for neutron, proton and 
a-emission become 


Q,, = 10 Mev., QO, =7 Mev., Ov = 9 MOVe 7 eae (78) 


However, when the excitation energy is 9 Mev., proton emission is very much 
more likely than «-emission. 

These results show that if after emission of the penultimate nucleon the 
excitation energy lies between 10 and 7 Mev. a proton is emitted, if it lies. 
below 7 Mev. y-rays are emitted. There is then a probability (10 —7)/10=0-3 
that the last heavy particle emitted by the nucleus is a proton with kinetic energy 
between 0-5 and 3-5 Mev.* 

On the average, with the nuclear parameters 5, this process contributes 0-3 
protons per star with energies in the range 0-5 to 3-5mev. With the parameters 
d or e the contribution is 0-15 protons per star with energy in the range 0:5 to 
1-5 Mev. 

Some very slow protons are observed experimentally. Classically this cannot 
be understood since Weisskopf’s formula (52) forbids emission of particles with 
energy below the barrier height V, in this case about 6Mev. Harding, Lattimore 
and Perkins (1949) have pointed out that slower particles can be emitted by 
ordinary barrier penetration, and this has been allowed for crudely by replacing 
V by V’ in (54). The penetration formulae of Bethe and Konopinski (1938) 
lead to V,’ =0-7 V, thus this process allows protons with energies down to about 
4Mev. Protons with lower energies can be attributed to the ‘trick’ mechanism 
described above. 

The mechanism of this section has been proposed independently by Fujimoto 
and Yamaguchi (1949). It is difficult to say from the short publication available 
how far their methods of calculation are similar to those of this paper ft. 

* At these low energies one should take account of the dependence of binding energy on the 


even—odd characteristics of the nuclei. However this leads to the same result. 
} The author does not agree with their statement that y-emission is negligible. 
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The experimental evidence shows that this cannot be the full theory of the 
emission of low-energy particles, for it does not explain the emission of low-energy 
a-particles. If one takes for the potential barrier V,=12Mev., the penetration 
effect gives V,’=10mev. but many «-particles of lower energy are emitted 
(H.L.P., Page (1950)). It seems necessary to assume with Bagge (1938)that the 
potential barrier is lowered at high excitation energies. 


3.7 Final Value of the Neutron Excess 
With A=70, the value —1 of 2B(@—0,) at x=0 gives N—N,=—0-75. 


Since on the average 0-2 protons are emitted in the subsequent process described 
in §3.6, one has then N=N,—0-5, a very slight neutron deficit. The residual 
nucleus, which may have 1 Mev. or so of excitation energy then reaches stability 
by y-emission or perhaps by a f-transition with K-capture. 

Thus, according to this theory, the way to estimate the unobservable neutron 
emission from a star is to calculate the charge of the final nucleus and to assume 
that its neutron content corresponds to maximum stability. The experimental 
data have usually been analysed by assuming that equal numbers of neutrons 
and single protons are emitted. The present author’s procedure would lead to 
higher estimates of the energy associated with stars of a given number of prongs. 

This prediction can to some extent be checked experimentally. For, if equal 
numbers of neutrons and protons were emitted, the residual nucleus would have 
a large neutron excess and would reach stability by B-transition with emission of 
electrons of about 1 Mev. Two or three electrons would be expected from a star 
of 10 prongs. These would be recorded in the new photographic emulsions but 
have not been observed. 


§4. NUMERICAL RESULTS 


4.1 Nuclear Parameters and Method of Integration 


The main experimental data to be compared with the theory concern ‘stars 
with an average of 10 prongs. Allowing for the preliminary ejection of a few fast 
particles, it is assumed that for this star size the evaporation proper starts at about 
A=98 and terminates at d=68. Thus during the evaporation the average 
value of A is about 83. The previous discussion has taken explicit account of 
all factors which depend on differences between values of A; therefore for the 
calculation of A which varies only as \/A it is sometimes sufficient to use the mean 
value A =83. 

The equations of emission have been integrated for various sets of parameters, 
namely 
aA=4 V,, =V4 =V, =5 Mev. Vy =V = 10 wey, 


D 
b A=5~1/(A/3°3) 
CN =5'5 ~~ 4/(A/2-7) 
d N=(A/3-3) | | Vi=0-7V, Va =0-77V, Vi=0-8V, Vi =1-6V, 
e A=/(A/2-9) ] V,/=1-66V and V=6/(1+U/200) ev. 


Ivy Va =Vi =4V{ =4V',,=5(1 — U/1000) mev. 
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The reduction of potential barrier with excitation energy assumed in 6 and c is 
roughly that predicted by Bagge (1938). The more drastic reduction d seems to 
be required to explain the emission of slow «-particles. ‘The relationship between 
V and V’ assumed in d,e is a fair approximation to the penetration formulae of 
Bethe and Konopinski. 

The evaluation of the term 28(0 —0,) was discussed in §3.5 and the results are 
plotted in Figure 1. These values can be inserted into the emission ratios 65. 
None of the ratios b, c, d, e, with the possible exception of c, is at all sensitive to 
the term in 6(0 —6,); therefore the dependence of this term on the initial excitation 
has been ignored and all the emission ratios (65) have been calculated using values 
derived from an initial excitation of 400mev. ‘The slight error in 9H emission 
thus introduced is of no importance since 3H emission is in any case small compared 
to that of protons and the experimental check (ii) on the emission ratios involves 
only (P,+P,)/=P,. Of course, (65a) is more sensitive to the 6 terms, but as 
the total neutron emission has already been determined, P, is required only for 
the calculation of XP, and SH,,P, and the approximataion is adequate for this 
purpose. ‘To illustrate these points the evaluation of (65 e) for P,/P, with para- 
meters b is given in the Table. The second line is the value of the expression 
in square brackets in that equation. 


U 30 50 100 200 300 400 500 600 
Bracket [] 1-76 1:99 2:95 4-56 6-11 7:6 8:95 9-8 
AN 
—— |]; 1-84 1:40 1-33 1:45 1-63 1:79 1:92 2:00 
Se i 
Py !Pp 0-31 0:49 0:53 0-47 0-39 0°33 0:29 0:27 


These values illustrate the existence of a maximum, already mentioned in 
§3.3, in the table of ratios P,/P, and also show the dominance of the term 
1-67 R/A at high energies. It is clear that small errors in the 6 terms are of no 
consequence. 

Integration of the equations of emission determines the total number of 
evaporated particles of each kind as a function of the initial excitation energy. 
The results are given in Figure 2, the contribution of §3.6 te the proton emission 
being included. ‘The full curves are shown as they give some insight into the 
evaporation process but the results may not be very accurate for stars of less 
than 6 prongs since the calculation has presumed initial conditions appropriate 
to larger stars. 

The theoretical ratio of the number of doubly charged tracks to the total number 
of tracks is plotted in Figure 4 as a function of star size. Experimental points 
given by Page (1950) are shown. The observed excess of «-emission for the small 
stars of six or less prongs is attributed to the disintegration of the light elements 
carbon, nitrogen and oxygen which are present in the emulsion. For stars of 
more than six prongs the agreement is satisfactory. 

H.L.P. ($46) also mention that the proportion of protons to all singly charged 
particles is 55/79=0-7. Thus the proportion of deuterons and 3H must be 
‘0-3 +0-07. For 10-prong stars the theoretical results are 


a=0-37, B= 0-28, C= (0:25; d=()-28, e=0-26, 


which is satisfactory enough. The number of singly charged particles ending 
m the emulsion is not sufficient for a precise check on this point. 


Number of Particles (7e) 
Total Number of Tracks (WN) 


0 - 0 
0 100 200 300 400 500 600 700 
Initial Excitation Energy U (Mev) 
Figure 2. The total number N of tracks and the numbers 7 of protons, deuterons, tritons, and! 


a-particles evaporated from a nucleus with initial excitation energy U. Parameters: mean 
of d ande. 
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The energy spectrum of the evaporated particles of charge e for five different values of 


Figure 3. 
the initial excitation energy. Parameters: mean of d and e 
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The energy spectrum of the emitted particles has been derived from equation 
.(54), which for this purpose can be put into the form 


PT) dT = ta exp{—-(T—V)fr}aT with 7r=(2V/R)JA, ..-.-. (79) 


2 
T 


by averaging over the whole emission process. The results are shown in Figure 3 
for various values of the initial excitation energy. All the singly charged particles 
have been grouped together since experimentally full identification of particles 
is not usually possible; similarly for the doubly charged particles. The slow 


0-40 


Q-35 


Emission Ratio 


0-20 


5 10 15 
Star Size (number of tracks) 


Figure 4. Theoretical and experimental values (Page) of the emission ratio (number of tracks of 
charge 2e/total number of tracks) as a function of star size. The rectangle indicates the 
mean and standard error of the experimental results for stars of 7 or more tracks. 


protons of §3-6 have not been included in the energy spectrum. In working 
out these curves it is worth while to take account of the dependence of A on A 
since initially, when the excitation is large, A and so A are above average ; this 
has the effect of reducing the high energy tail of the integrated energy spectrum. 

The experimental data for stars with an average of 10 prongs, due to H.L.P. and 
Page, are shown in Figures 5 and 6 for comparison with the theoretical curves. 
These stars have on the average one fast particle which falls outside the 
evaporation distribution, and the curves therefore represent the theoretical 
evaporation spectrum for 9-prong stars. The amount of experimental data 
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available for the «-spectrum is exaggerated in Figure 6 since the observations of 
H.L.P. have been increased by their correction factor f to allow for the ‘loss’ of 
long tracks. 


4.3 Discussion of Parameters 


a. The a prior’ parameters a give a surprisingly good estimate of the pro- 
portion of doubly charged tracks. The value is a little high, and the energy 
spectrum (not plotted in Figure 5 but worse than 4) contains too many high energy 
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Theoretical and experimental (H.L.P. and Page) results for the energy spectrum of 


Figure 5. 
The slow 


particles of charge e evaporated from stars with an average of 10-2 prongs. 
‘ decay protons’ from the process of §3.6 are included. 


particles. Both these results show that the A value is too small. Attempts to 
correct the emission ratio, with this value of A, by raising the potential barrier 
would lead to a very bad energy spectrum. 


b, c. Increased values of A are combined with a reduced potential barrier to 
give a small decrease in the emission ratio and a greater improvement in the energy 


spectrum. Since for the ground state of the nucleus the height of the potential 
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barrier is approximately known, the reduction must occur at high excitation 
energies only. This provides evidence, independent of that previously derived 
from the « spectrum, in support of Bagge’s calculation. Figure 5 shows that the 
theoretical energy spectrum still contains rather too many fast particles of charge e. 
Curve c is better than 6 and larger values of A would lead to further improvement 
but Figure 4 shows that the emission ratio would then be much too small. As the 
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Figure 6. Theoretical and experimental results for the energy spectrum of evaporated particles. 
of charge 2e. Stars as in Figure 5. 


a spectrum (Figure 6) contains far too few particles of low energy, a further 
reduction of the potential barrier is suggested and this will at the same time 
improve the spectrum of the singly charged particles. 

d, e. Parameters d and e have been introduced to improve the energy spectrum. 
There is now sufficient emission of low energy «-particles. In contrast, it will be 
noticed that the lowered potential barrier does not greatly increase the emission 
of slow protons since the increased contribution from states of high excitation is 
partly compensated by a reduced contribution from the process of §3.6; the high 
energy side of the spectrum is however much improved. 

Figure 4 shows that, for stars with over 7 prongs, the calculated emission ratio 
decreases linearly with star size. ‘This simplifies the comparison with experiment. 
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Page (1950) has determined the mean value of this ratio as 0-271 + 0-008 for 392 
Stars with an average of 10-2 prongs. For this size d gives 0-278 and e gives 
0-260. ‘Thus, on this evidence, the correct value of A lies midway between those 
assumed in d and e which represent the limits of uncertainty of this parameter. 
This conclusion is in full agreement with the independent evidence from the 
energy spectrum. 

The experimental data therefore indicate a value 3:1+0-2mev. for the 
coefiicient Cy/z? which determines A. This value falls well within the limits 
arrived at in §2.5. 

The value V = 6 mev. for the residual nucleus is reasonable. But the reduction 
of potential barrier with increase of energy assumed in d and e is much larger than 
would follow from the small uniform expansion calculated in §2.5 and a little 
larger than that derived by Bagge as a consequence of the formaticn of surface 
waves of large amplitude. However the values are certainly required by the 
experimental data and are not a priori unreasonable. ‘The dependence on energy 
could be less marked if one started from a lower potential barrier in the ground 
state of the nucleus and assumed the same value for the excited nucleus. The 
author’s value V=6 Mev., is actually higher than Bethe’s value V = 5-3 Mey. derived 
from his larger estimate of the nuclear radius. The above estimate of €y would 
also be consistent with a slightly larger nuclear radius, which would have to be 
combined with a value of y nearer to the upper limit of §2.5. However this work 
provides no grounds for a choice of nuclear radius. 


§5. CONCLUSIONS 

A detailed investigation of the nuclear evaporation theory shows that it can. 
give a consistent account of the main features of disintegration stars. 

It has been necessary to make rather ad hoc assumptions about the potential 
barrier but comparison of the results for b, c and d, e shows that the experimental 
estimate of the parameter A, which (equation (37)) connects nuclear energy and 
temperature, does not much depend on these assumptions. ‘T’hus it is legitimate 
to check the theoretical estimate of A against the experimental data. The agree- 
ment is as good as can be expected. 

This investigation yields information on one point of principle. Factors to: 
represent explicitly the ‘probability of existence’ of compound particles («, 
triton) within the nucleus do not enter explicitly into Weisskopf’s evaporation 
formulae. The ‘existence’ is completely taken into account through the depen- 
dence of the emission probabilities P,, P, upon the binding energies w,, w, of these 
compound particles. In the formulae, such as (56), this effect is expressed by the 
dependence of R on w. The calculations support this point of view, for it 
would not have been possible to obtain agreement with experiment if additional 
factors to represent small ‘probabilities of existence’ had been introduced into 
the emission formulae. The Table shows that to maintain the correct proportion 
of doubly charged particles, a factor } for their existence within the nucleus would 
have to be compensated by halving the value of A. This would give an incorrect 
energy distribution, the width of the proton spectrum would approximately be 
doubled. The situation in this problem is quite different from that existing in 
the theory of natural «-radioactivity (Bethe 1937) where, because of the dominance 
of the penetration factor and its critical dependence on barrier height, large 
factors against the existence of the «-particle can be ccmpensated by small changes 
in the nuclear radius. 

PROC. PHYS. SOC. LXIII, 3—A 19 


282 C. H. Collie, P. F. D. Shaw and H. F. Gale 


ACKNOWLEDGMENTS 
The author would like to thank Dr. Page and Dr. Perkins for the generous 
supply of experimental results and Professor Rosenfeld for his kind interest in 
this work. 


REFERENCES 


Baccr, E., 1938, Ann. Phys., Ltz., 33, 389; 1943, Cosmic Radiation (Ed. Heisenberg) 
(reprinted 1946 by Dover Publications, New York), p. 128. 

BARDEEN, J., 1937, Phys. Rev., 51, 799. 

BetuHe, H. A., 1937, Rev. Mod. Phys., 9, 69. 

Berue, H. A., and Bacuer, R., 1936, Rev. Mod. Phys., 8, 82. 

Betue, H. A., and Konopinsk1, E. J., 1938, Phys. Rev., 54, 130. 

Bour, N., and Karcxar, F., 1937, K. Danske Vidensk. Selskab, 24, to. 

FEENBERG, E., 1947, Rev. Mod. Phys., 19, 239. 

FERRETTI, B., 1948, Nuovo Cim., 5, 36. 

Fujimoto, Y., and YAMaGcucHI, 1949, Phys. Rev., 75, 1276. 

Harvinoc, J. B., Larrimore, S., and PERKINS, D. D., 1949, Proc. Roy. Soc. A, 196, 325. 

Pace, N., 1950, Proc. Phys. Soc. A, 63, 250. 

RoseNFELD, L., 1948, Theory of Nuclear Forces (Amsterdam: North Holland Publishing 
(COs): 

STroner, E. C., 1936, Phil. Mag., 21, 159. 

WEIsskopr, V., 1937, Phys. Rev., 52, 295. 

ZWANIKKEN, G. C. J., 1948, Physica, 14, 530. 


Note added in proof. My attention has been drawn to a paper by Ferretti (1948) 
on meson-induced disintegrations, in which the importance of the neutron 
excess in nuclear evaporation is discussed. 


The Correction of Specific B-Radioactivity for Self-Absorption 


By C..H. COLUTIE. POF: Dy SHAW aAnpeHw la Gis 
Clarendon Laboratory, Oxford 


MS. received 4th August 1949 


ABSTRACT. The self-absorption of B-rays in radioactive films between 5 and 100 mg/cm? 
superficial density has been investigated. The marked departures from the exponential 
absorption law are shown to be due to scattering of the f-rays in the film. The use of very 
thin films to avoid the correction for self-absorption is not recommended as a routine method. 


§1. INTRODUCTION 


N using f-radicactive elements as tracers one often needs to correct observed 
] readings (particles per minute per mg.) for self-absorption in the specimen to 

a specific activity fora thin film. This is particularly important in modern 
work using end-on counters since the restricted window area leads to the use of 
thicker preparations than were commonly used with cylindrical counters. An 
empirical correction curve can be constructed, but this can be very time-consuming 
if the specific activities involved are low, and a good understanding of the processes 
contributing to the self-absorption is very desirable. Extensive experimental 
work on this subject (Broda et al. 1948) has shown that for thick films an expo- 
nential absorption law is well obeyed. In the course of some measurements using 
radioactive silver iodide (Shaw and Collie 1949) two of the present authors found 


+ 
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that no appreciable correction for self-absorption was necessary below a super- 
ficial density of 50 mg/cm’, a result which greatly simplified the experimental 
work. It was later realized that this result was not in agreement with the expo- 
nential absorption usually assumed, according to which the correction should 
have been 15°. An examination of the published data showed that few measure- 
ments had been made below 25 mg/cm?, and that most of the work had been done 


on films thicker than 50 mg/cm?. ‘The object of the present work is to investigate 
this anomaly. 


§2. SELF-ABSORPTION OF THIN FILMS 


‘To measure the self-absorption a stock of dry radioactive powder was prepared 
and suitable amounts were weighed out into the specimen holder. This 
consisted of a brass plate with a lcm. diameter hole in it closed by a piece of 
‘Cellotape’ placed across the bottom of the plate. After introduction of the 
powder the cavity was closed with a further piece of ‘Cellotape’ and the whole 
pressed with a well fitting plunger to make the preparation as homogeneous as 
possible. Specimens made in this way gave results which were well reproducible. 
The specimen holder fitted into a lead castle below a G.E.C. end-on counter 
(Type G.M. 2) for measurement. The usual corrections for decay and paralysis 
time were made, and the sensitivity of the counter was controlled by the repeated 
measurement of an uranium oxide source. ‘Iwo positions for the source were 
available, an upper and a lower one, distant 0-4 cm. and 2-4cm. respectively from 
the window of the counter. 

The results obtained are shown in Figure 1, curves ) and g, and disclose a very 
shallow maximum in the self-absorption curve. Similar curves, c, a and f, were 
obtained using films of radioactive silver bromide and barium phosphate 
respectively. It was at first thought that the phenomenon was to be explained in 
terms of some specific property of silver iodide which is reported (Broda et al. 1948) 
to show an absorption curve anomalous in other respects. However, since 
radio-phosphorus emits no y rays, and radio-bromine and radio-iodine have very 
different B-ray spectra, the phenomenon is clearly a general one, not to be explained 
in terms of some specific property of silver iodide. ‘The maximum is sufficiently 
shallow to require rather accurate counting to establish it with certainty, but an 
examination of the self-absorption curve for radio-sulphur, *°S, in barium 
sulphate obtained by Hendricks et al. (1943) shows an ill-defined maximum at 
about 0:3 mg/cm?. A similar phenomenon has been briefly reported by Solomon 
et al. (1947) for thin films of barium carbonate, Ba'*CO;. The latter at first 
attributed the phenomenon to the loss of the very small amounts of barium 
carbonate (0:65 mg/cm?) involved, but Glendenin (see Solomon et al. 1947) 
suggested that the observations could be explained in terms of scattering. 

The existence of amaximum indicates the interplay of two opposing phenomena, 
the one causing an increase in specific activity with superficial density and the 
other, which is presumably the exponential absorption dominant for thicker 
films, causing a decrease. 

It was at first thought that the phenomenon causing the increase was back 
scattering. For a thin film the electrons which are counted are those initially 
projected forward together with those scattered back from the support, which is 
either ‘Cellophane’ and air or ‘Cellophane’ and brass, according to the method 
of supporting the film. Since a large part of electron scattering is nuclear and 
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so, according to Rutherford’s law, is proportional to Z?, the back scattering in 
silver iodide is more than in brass, and it is possible that this might account for the 
phenomenon. ‘The exponential external absorption of f rays is sufficiently well 
established by experiment for it to be usual to quote the absorption coefficient 
u/p in the literature. The theoretical explanation of this is uncertain since a 
stream of homogeneous electrons is not absorbed exponentially. ‘The conditions 
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Figure 1. Self-absorption curves. 

a. **P Bs in barium phosphate (lower position, no backing). 

b. 8] 6s in silver iodide (upper position, brass backing). 

c. **Br fs in silver bromide (upper position, brass backing). 

d. Calculated self-absorption curve for 1°81 f’s in silver iodide assuming exponential absorption 
(upper position). 

e. Calculated self-absorption curve for 181 f’s in silver iodide using measured external absorption 
curve (upper position). 

f. ®*°P B’s in barium phosphate (upper position, no backing). 

g. ‘*8I B's in silver iodide (upper position, no backing). 


under which exponential absorption will be observed have been discussed inter 


alia by Feather (1944). The occurrence of exponential external absorption over 


a limited range accounts for the observed self-absorption, provided that all effects 
due to scattering are neglected, as they can be for thick films. A theoretical 
treatment of the effect of scattering is extremely complicated, but the self-absorption 
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curve can be constructed by graphical integration if the external absorption curve 
and the back-scattering curve are known. Accordingly the effect on the observed 
activity of a thin film of active silver iodide of backing it with increasing thicknesses 
of inactive silver iodide was measured. The maximum amount of back scattering 
is approximately 26°%, and the shape of the curve agrees generally with what was 
to be expected from the work of Yaffe (1949) and Zumwalt and Lloyd (1947). 
By combining this curve with an exponential external absorption curve one can 
determine what the self-absorption curve should be. Imagine a layer of active 
material divided up into m equal sections, and consider one of the inner sections 
of activity 27. From the mth layer radiation J will proceed towards the counter ; 
%,,1 will reach the counter, and 8,,J will be scattered back. ‘Thus the downward 
radiation from the layer is /(1+8,,), of which /(1+8,,,)y,, will be scattered back. 
The total radiation reaching the counter which started in the mth layer is thus 
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and the total activity of the film is given by 


m=n 1 +B ) 
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in which expression «,, is to be read from the external absorption graph and 
f,, and y,, from the graph of back scattering against superficial density at the thick- 
ness corresponding to the value of m. ‘The result of such a graphical integration 
{using 5mg. steps) assuming exponential absorption in the layer is shown in 
Figure 1, curve d. Although the effect of taking back scattering into account is 
seen to be considerable, there is no maximum in the curve, as is found experi- 
mentally. The assumption that the external absorption remains exponential for 
thin layers of absorber can be investigated by direct experiment. This was done 
using weighed layers of silver iodide and a thin active silver iodide source. In the 
upper position the use of silver iodide was impracticable and sheets of tin foil, 
which has the same scattering power, were used instead. The results are shown 
in Figure 2, curves a, b, c; it is at once apparent that the simple absorption law 
obeyed for thick absorbers breaks down for thin ones, and indeed that the effect 
of a thin absorbing foil is to increase and not decrease the measured activity. 

If the experimental absorption curve is combined with the observed back- 
scattering curve to give a calculated self-absorption curve, the result, shown in 
Figure 1, curve e, is seen to be in as good agreement with experiment as could be 
expected when the rather coarse steps in the graphical integration are taken into 
account. 

It is quite easy to account for the maximum in the external absorption curve 
in terms of electron scattering. Electrons are heavily scattered in layers of material 
in which the absorption is small. It is therefore not unreasonable to consider a 
small thin film in which absorption can be neglected. If the film is thin enough 
to allow scattering to be neglected as well, the polar diagram of the electrons 
emitted into asolid angle Q will be independent of the angle 0 between the direction 
of and the normal to the film. 

For a thicker film the angular distribution of the electrons will approximate 
to the usual cos@ distribution: since there is no absorption in the film the total 
number of electrons emitted is constant. Asa result of this, a counter subtending 
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any angle less than 27 at the film will count a greater fraction of the electrons from 
a thick film, and in the extreme case of a counter subtending a small solid angle 
about the normal to the specimen the observed activity will be increased, due to 
scattering, by a factor two. This is a large effect well able to mask the effect of 
absorption and so give a maximum in the absorption curve. ‘The general 
correctness of this explanation is confirmed by the increased sharpness of the 
maximum which is shown when thin foils of platinum (which has a larger atomic 
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Figure 2. External absorption curves. 


vase B's in silver (lower position, brass backing). 

fait B's in silver iodide (lower position, no backing). 
ae B's in tin (upper position, no backing). 

oe B's in silver (lower position, aluminium backing). 

*2P B's in platinum (lower position, aluminium backing). 
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number and, therefore, shows more large angle scattering) are used as the absorber 
(Figure 2,curvesd,e). ‘The matter may be looked uponin another way. Electrons 
which start off towards the counter have a good chance of getting there because the 
thickness of film through which they have to pass is small. On the other hand 

electrons emitted nearly parallel to the film are quite likely to be scattered ‘hrouen: 
a large angle out of the film into the counter. Thus increasing the thickness of 
the scattering layer will increase the recorded count until the film is thick enough 


; 
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for the probability of scattering to be large whatever the initial direction cf the 
B particle. To verify that this effect was taking place in films under consideration 
whose superficial density lay in the range 0-50 mg/cm? the apparatus shown in 
Figure 3 was constructed. 

It consisted of an evacuated space in which a small source could be mounted 
as freely as possible on a light wire frame, which could be rotated on a ground 
Joint. Electrons in a selected solid angle entered an external counter after 
traversing the thin mica window at the top of the apparatus. An aluminium 


Figure 3. Apparatus for measuring the angular distribution of electrons from a soutce. 


A. Geiger-Miller counter. B. Diaphragm. C. Source. D. Aluminium frame. E. Ground 
joint. F. Scale of degrees. G. Inspection window. H. Vacuum pump. 


diaphragm ensured that electrons scattered from the aluminium-lined walls 
could not enter the counter unless they had been scattered at least twice. The 
sources consisted of active barium phosphate on a lcm. diameter filter paper 
mounted on a ‘Cellophane’ strip. Coherence was given to the specimens by a 
final wash with dilute collodion, and the superficial density was determined after 
the measurement by incineration in the usual way. 

The electron distributions obtained are given in Figure 4 and show the 
expected change in the angular distribution cf the electrons from independence 
of @ to dependence according to a cos @ law in just the range of superficial densities 
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considered. There can be no doubt that this phenomenon must be taken into 
account in explaining the observed self-absorption curves, and that it will give 
rise to a maximum if the true absorption is small and the solid angle subtended 
by the counter is appreciably less than 27. 


$3, DISCUSSION 
The accuracy with which the polar diagrams represent the true angular distribution 
of the electrons from the source is rather difficult to assess. Removal of the 
diaphragm increased the readings by about 30° and began to affect the shape of 
the curves at about 60° azimuth. This indicates that the bafHe is effective in 
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Figure 4. Angular distribution of B-particles from a plane source. 


Each quadrant refers to a separate experiment, and @ increases uniformly from 0-90° clockwise. 
Readings were taken every 10° and are represented by radial lines whose length represents 
the statistical error of the readings. 


(a) *P on 12 « * Styrafoil’. (6) **P in barium phosphate 8 mg/cm?. (c) *2P in barium phosph 
Hesse : : 
22 mg/cm’. (d) **P in barium phosphate 70 mg/cm?. Full line represents a cos 0 aiscaeceoee i 


preventing scattered electrons from entering the counter. Another test was made 
using a very thin active phosphate film made by evaporating a drop of separated 
phosphate on to a 12 ‘Styrafoil’ sheet. This was supported on a frame in such 
a way that the scattering due to the frame when the film was nearly vertical was 
intercepted by the walls and the baffle. The results obtained are shown in 
Figure 4(a). ‘The bulge between 45° and 80° represents the effect of scattering in 
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the film, of electrons initially directed in its plane. Evidence in favour of this 
view can be obtained by a graphical integraticn of the polar diagram to obtain the 
total number cf electrons emitted from the source. If w is the mean solid angle 
subtended by the counter at the source, and 7 is the rate of counting of particles 
in the spherical part of the polar diagram, then the total number of particles 
emitted forward from the source is 


Aare. 2. 
he —— (1388 + 8) electrons/minute. 


This number can also be cbtained by graphical integration of the complete polar 
diagram. This gives (27/w) (1387 + 8) particles/minute. Agreement within the 
statistical error is fortuitous, but does show that general scattering from the walls 
of the containing vessel cannot be seriously distorting the true angular distribution 
of the f particles. 

It remains to see to what extent one can reconstruct the self-absorption curve 
from the now empirically known elements of absorption scattering and back 
scattering. 

The most suitable case for a quantitative comparison is that of barium phosphate 
observed with the source some distance from the counter, for which the self- 
absorption curve is given in Figure 1, curve a. ‘The maximum occurs at about 
80 mg/cm? when the apparent specific activity is about 30°% above that found for 
a very thin film. Using our empirical data one can calculate the solid angle 
which the counter must subtend at the specimen to give this result and compare 
this with the measured solid angle. 

The thin source (of activity say 100 units) will be giving an approximately 
spherical distribution while the thick sources will have passed over to the cos @ 
distribution. The 130 units of radiation observed with the thick source have 
been reduced by self-absorption, and if the film had been an ideal one showing only 
scattering the observed activity would have been 157 units. ‘This increase is, 
on our assumptions, entirely due to changes in the angular distribution of the 
electrons. 

Let @ be the effective angle subtended by the counter and let 47/ f particles 
be emitted from the ideal source. Then for a uniform angular distribution 
2nI(1—cos@) will enter the counter and after redistribution by scattering 
7I(1 —cos 26) willenter the counter. Thus in the case being considered 4 is given 
by the equation 

a1(1 —cos 28) —157. 
2mI(1 —cos 8) 
whence 0=35°. 

The measured value of 6 is approximately 26°. This agreement is as close 
as can be expected when the approximate nature of the assumptions and of the 
methods of calculation involved is considered: for example, no correction was 
made for scattering due to the brass plate surrounding the source; a reasonable 
allowance for this makes the agreement between measured and calculated angles 
much closer. 

It remains to discuss the practical application of these results to tracer work. 
It is quite clear that the specific activity obtained by correcting a measurement 
made on a thick preparation of moderate activity is different from, and cannot be 
compared at all easily with, the specific activity measured on a thin specimen. 
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It appears that there is no particular advantage in working with extremely thin 
preparations as this does not avcid the need for applying a correcticn. The mest 
satisfactory prccedure seems, in comparative work, to be to measure near the 
maximum in the apparent specific activity curve. By arranging for the ccunter 
to subtenda large sclid angle this maximum can be made so flat that except for 
work of very great accuracy using very strong preparations (to avoid the statistical 
errcr) no correction is needed. 
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The Magneto-Resistance of the Alkali Metals 


By D. K. C. MacDONALD 
Clarendon Laboratory, University of Oxford 


Communicated by F. E. Simon; MS. received 4th January 1950 
ABSTRACT. The magneto-resistive effect has been examined experimentally in the: 


alkali metals at low temperatures. It is concluded that deviations from the free electrom 
gas model become considerable in the case of the heavier metals, caesium and rubidium.. 


electrons in a magnetic field (e.g. Titeica 1935), then the change of resistance: 


I one neglects the possibility of quantization of the orbits of the conduction. 


of a metal in a magnetic field arises from the departure of the electronic 
behaviour from that of the ideal free electron model. Sommerfeld and Frank 
(1931) analysed the very small effect that does in fact arise for free electrons, 
but Peierls (1930, 1931) first showed that the observable change springs from the 
non-ideal behaviour which may also be expressed as a departure of the Fermi 
surface from the spherical surface in momentum space characteristic of a free 
electron gas. ‘The detailed analysis of the magneto-resistive effect is still today 


by no means complete, but measurements of its magnitude should offer a means. 


of comparing how far different metals of similar type depart from the free electron 
model. 


I have examined experimentally the relative magneto-resistive effect in all 


the alkali metals; these would be expected a priori to present the closest. 
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approximation to the ideal case. The magnitude of the change anticipated 
depends essentially on the ratio of the field strength, H, to the electrical resistance, 
Ra, at the temperature 7. This dependence, often known as Kohler’s Rule, 
was indicated originally by Peierls (1930, 1931) and by Bethe (1931). Low 
temperatures are then particularly valuable in yielding a satisfactory magneto- 
resistive change even for moderate fields (103-104 gauss) because of the low value 
of resistivity attainable in pure samples of metal. The experiments were carried 
out at temperatures between 4-2 and 20°K., and the results are shown in the 
accompanying figure. The ordinate exhibits the percentage change in 
resistance observed, 100(AR,/R,), while as abscissa we have H/r, where rp is. 
the ratio R7/Ro.,. To provide a perhaps more valid comparative abscissa for- 
different metals one might use instead of 7, the parameter t,=R/,R, where 
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Magneto-resistance of the alkaii metals. 


@ is the Debye temperature for the metal concerned; difficulty arises, however,. 
due to the uncertainty in the value of ©. Justi (1948) has used this parameter’ 
in a recent publication in this field, but unfortunately so many slips occur in 
the relevant part of his paper that useful comparison of the results has not been 
found possible. Another choice of abscissa which suggests itself (cf. MacDonald 
and Mendelssohn 1950) is obtained through multiplication of rp by the factor 
Q/Iog., Where Q is the atomic volume of the metal concerned and / the electronic 
mean free path. This factor yields the values below for the alkali metals, when 
normalized to unity for sodium: 


Metal Iki Na K Rb &s 
Oars 1-6 (EO ITA 3:6 6:2 


The modification of the abscissa, however, in either of these ways only serves 
to emphasize more strongly the difference between the alkali metals demonstrated 


in the Figure. 
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The results on all specimens except Nal0 may be regarded as essentially 
appropriate to polycrystalline metal since no appreciable variation in AR,/Rp 
was observed with angular rotation of the specimens in the field. 

The increasing departure from the free electron model as one goes from 
sodium towards the heavier members of the alkali group is in agreement with 
the conclusions drawn from detailed experimental comparisons of the low 
temperature variation of the conductivity with the theoretical Bloch—Grtineisen 
law by MacDonald and Mendelssohn (1948, 1950). It was found there in 
particular that the theory was obeyed very closely by sodium, but deviations— 
becoming particularly strong in the case of rubidium and caesium—were 
observed in all the other metals. 

Further work in this field, particularly on single crystals, is at present in 
progress, having as its ultimate goal a more precise knowledge of the characteristics 
of the Fermi surface structure. 
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The Radiations of Praseodymium-142 


The 19-3 hour activity of 142Pr has been investigated by a number of authors (De Wire et 
-al. 1942, Cork et al. 1948, Jurney 1949, Mandeville 1949) who have given rather inconsistent 
results for the energies of the beta and gamma radiations. This activity, in a sample of 
spectroscopically pure Pr.O, irradiated in the Harwell Pile, has been examined here in a 
semicircular magnetic spectrometer of 12 cm. radius and geometrical resclution of 2°% in 
the course of other work, and the results of this investigation are given below. 

Figure 1 shows the beta-ray spectrum which, when examined by means of a Fermi plot, 
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Figure 1. The beta-ray spectrum of !42Pr. 


gives end points for the soft and hard components of 0:66-+0-02 and 2:23 +0-04 Mev. 
respectively. Since a very thin source was used, and since the counter window was of mica 
of the order of 1 mg/cm? thick, the soft end of the spectrum should not be appreciably 
distorted and a comparison of the areas under the two components indicates that the hard 
beta-rays outnumber the soft by a factor of about four to one. It is also interesting to note 
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that the end point of the soft component (0°66 MeV.) is not inconsistent with the experimental 
points of Jurney from which he derived an end point of 0-35 Mev., especially as it can be 
shown that, with a coincidence circuit of any reasonable resolving time, it would be extremely 
difficult, using an absorption method, to determine the end point of a spectrum of this 
nature because of the very large number of accidental coincidences due to the background of 
hard beta-rays. In the case of the hard component of the beta-ray spectrum, the agreement 
between the experimental values for the end point given by the various authors is much more 
satisfactory and the value obtained here is in reasonable agreement with the only other 
spectrometer determination, that of De Wire et al., of 2:14-0:02 mev. 

The spectrum of the secondary electrons ejected from a thick copper cylinder covered 
with a thin (0-0015 in.) lead foil by the gamma-rays of “2Pr is shown in Figure 2 and the 
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Figure 2. Momentum spectrum of electrons ejected Figure 3. Energy level diagram for the 
from a thick-walled copper cylinder covered with beta-decay of 14?Pr. 
a thin lead foil by the gamma-rays of !4?Pr. 


estimated energy of the hard gamma-ray is 1:°59-+-0:04 mev. There is a suggestion of the 
presence of a weak soft radiation with an energy of about 0-14 Mev., and this is confirmed by 
the presence of a very weak conversion peak in the beta-ray spectrum which gives the energy 
of this soft line as 0:135-+-0-005 Mev. and shows that even if this is of electric dipole origin, 
its intensity cannot be greater than about 20°% of that of the hard gamma-ray. Since this 
peak appears to decay with a similar period to that of the other radiations, it seems likely 
that it is, in fact, the same line as that reported by Cork etal. as having an energy of 0:134Mev., 
the conversion lines of the harder radiations at 0:329, 0:490 and 0-624 Mev. being invisible 
in this instance on account of their much smaller conversion coefficients. 

Since the four soft gamma-rays reported by Cork et al. add up in energy to 1:58 Mevy., it 
seems likely that the transition from the level at 1:59 Mev. to the ground state of the *Nd 
nucleus may either be made directly or else by way of three intermediate levels, thus giving 
rise to four soft gamma-rays in cascade. ‘These cascade gamma-rays would account for the 
small number of gamma-—gamma coincidences reported by Mandeville and Jurney. It is 
also possible that there are subsidiary beta-ray transitions directly to these intermediate 
levels which would be too weak to be observed against the strong 2:23 Mev. component. 
It would therefore appear that the complete decay scheme must be similar to that shown in 
Figure 3 in which one possible arrangement of the intermediate levels is given. ‘The tran- 
sitions actually observed are indicated by full lines and the possible subsidiary beta-ray 
transitions by dotted lines. 


Department of Natural Philosophy, E. R. Rar. 
University of Glasgow. 
29th November 1949. 
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The Angular Distribution of Gamma-Rays from ’Li(p, p’)’Li* 


The nuclear level of 7Li at 0-48 Mev. is of considerable theoretical interest; owing to its 
relative simplicity, this nucleus should form a useful test case for any proposed model of 
nuclear structure. ‘The most convenient interpretation of the low-lying level is that it forms 
a 2P spin doublet (Inglis 1936) with the ground state, which is known to have spin See 
On this basis the spin of the excited state should be 1/2. There are, however, some experi- 
mental facts and theoretical considerations difficult to reconcile with this view (Devons 1949, 
Hanna and Inglis 1949). 

To study the problem further, the angular distribution of gamma radiation from the 
reaction 7Li(p, p’)’Li* has been measured. If the spin of “Li* is indeed 1/2, this distribution 
must be isotropic in the centre of mass system, whereas in general some correlation between 
the directions of the bombarding proton and the gamma-ray would be expected for all spin 
assignments other than 1/2. 

The gamma detectors used were scintillation counters, and the bombarding energies on 
the thin (~80 kev.) Li,CO; target ranged from the effective threshold value (800 kev.) to 
1,080 kev. Under these conditions the contribution to the counting rates from the 17 Mev. 
gamma-ray due to ’Li(p, y)®Be was negligible. 
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The accompanying diagram shows the experimental results at 980 kev.; essentially the 
same distribution was obtained at the other bombarding energies used. The points have 
been superimposed on the curve [(@)=I,(1+ A cos @), with A=0-016; this would be the 
expected distribution in laboratory space if the gamma-ray were emitted isotropically from a 
nucleus still moving with the full ‘ recoil’ velocity. The results of Elliott and Bell (1948) 
indicate that the nuclei are not slowed down much in the target material during the life-time 
of the excited state. Approximate isotropy is assumed for the inelastically scattered 
protons, so that the momentum carried away by them is, on the average, zero. 

The experimental points fit the curve well, so that the upper limit of any undetected 
anisotropy can be placed at A=0-005. Although a full theoretical analysis has not yet been 
made, it would appear difficult to reconcile this result with any spin assignment to 7Li* 
other than 1/2. 

Full details of this work will be published in due course. JI am indebted to Dr. S. Devons 
for suggesting the problem and for valuable discussion. 


Cavendish Laboratory, ROM. cee 
Cambridge. ; : 
16th December 1949. 
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On the Production of Penetrating Showers 
in Nuclear Collisions 


In a recent paper Heitler and Janossy (1949) have given a rigorous theory of meson 
production, and applied it to the absorption of the particles that produce local penetrating 
showers. Their calculations are based on the assumption that the probability w of an 
energy loss € by a nucleon of energy £ in a collision follows a law of the form 


we, EB) de=w(e/E) de/E, 


and it is shown that in order to comply with experimental results z must be a fairly broad 
distiibution, favouring rather low values of «€/E. This assumption for w was suggested by 
analogy only; other laws are being considered as well. It is the purpose of this communi- 
cation to point out additional evidence which could be interpreted as supporting that 
particular choice. 

In the analysis of experimental results on penetrating showers, in general only two 
factors are taken into account: the fact that while the collision cross section for fast nucleons 
is equal, or very nearly equal, to the geometrical nuclear cross section, the absorption 
cross section is only about one-half of that value; and the absence of an appreciable geo- 
magnetic effect (Appapillai and Mailgavanam 1948), at least in sea-level observations. One 
would believe that the first is explained by any theory which permits the average primary 
to undergo about two collisions before it becomes incapable of further shower production, 
and the second by any theory in which the bulk of the shower-producing particles at sea 
level are descendants of high-energy primaries capable of surviving repeated collisions with- 
out too much energy loss. A third criterion appears to be more sensitive, namely, the total 
intensity of the shower-producing radiation. 

Recent experiments carried out at an altitude of 3,260 m. gave a result of 


(1:-6+0-3) x 10-5 cm~? sec sterad—! 


for the vertical intensity of the charged particles which produce, in lead, showers containing 
at least three particles capable of penetrating more than 200 gm/cm? Pb. ‘This value, while 
slightly larger than the estimate of Hazen, Randall and Tiffany (1949), seems to agree well 
with experiments both at sea level and mountain altitudes of other authors (Broadbent and 
Janossy 1947, Cocconi 1949), if in comparing them inlot’s (1948, 1949) value for the 
absorption mean free path in air is used. However, extrapolation with the same absorption 
coefficient to the top of the atmosphere gives a’primary flux of shower-producing particles 
of only 0:55 10% cm~? sec™ sterad~t, or only about one-tenth of the total primary 
intensity. Thus it would appear that even if no multiplication process of nucleons takes 
place, only a fraction of all primaries give rise to penetrating showers. In fact, however, we 
know from the existence of such showers produced by neutral particles that fast nucleons are 
ejected in nuclear collisions : a process disregarded in the computations of Heitler and 
Janossy, but probably of considerable importance. 

This surprisingly small intensity could be explained in three ways : Firstly, low-energy 
primaries, that is, primaries of less than about 1:5 x 101° ev. if the integral energy spectrum 
follows a power law with an exponent of the order of 1:8, may be more rapidly absorbed 
than faster particles, and may thus not contribute to showers observed even at mountain 
or aeroplane altitudes. Secondly, the total energy which a shower-producing particle 
dissipates in the production of a penetrating shower in lead may be much higher than the 
2-3 -< 10° ev. observed in the form of charged szcondaries. ‘Thirdly, particles in the lower- 
energy region may undergo predominantly collisions in which no hard secondaries are 
emitted, so that again their contribution to the production of penetrating showers is 
negligible. 

The first assumption of more efficient elimination of slower than of fast particles could 
be understood if in each collision on the average there was dissipated a constant amount of 
energy rather than a more or less constant fraction of the incident energy. However, this 
would make the absorption by nuclear collisions of the primaries resemble ionization 
losses, and would not lead to an exponential absorption curve. Only with very artificial 
assumptions could one combine ‘ selective ’ absorption of lower-energy primaries with an 


exponential absorption law. 
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The second assumption would obviously require a much larger fraction of the total 
incident energy to go into neutral instead of charged particles. Even if all observed charged. 
shower particles were protons, the inclusion of neutrons could only about double the 
amount of 2-3 x 10° ev. carried by the charged particles; a large fraction of neutral mesons 
would be required to account for the rest, and owing to the decay of these particles all 
penetrating showers would be accompanied by an energetic cascade shower. While such 
events are known to occur, they are by no means common and typical. 

There remains, finally, the third assumption which demands that the slower primaries 
in the majority of their collisions will not transfer sufficient energy to generate a penetrating 
shower. ‘The same fractional energy transfer ¢/E would, however, suffice for the production 
of an observable penetrating shower if the primary energy was high. In other words : 
while every collision of fast particles initiates a local hard shower, only a small fraction of 
the collisions of slower particles does so, while the rest of these collisions will be of a more 
star-like character. This assumption corresponds exactly to the distribution function of 
Heitler and Janossy, with large probability for small e/E. Hence the experimental evidence 
that only a fraction of all incident primaries appears to produce penetrating showers lends 
support to that rather than to other transfer laws. 


Physics Department, Syracuse University, Ke Sires 
Syracuse, New York. 
28th December 1949. 
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Some Decay Effects of Mesons in Penetrating Showers 
observed in a High-Pressure Wilson Cloud Chamber 


Penetrating showers produced by cosmic-ray particles are now known to contain a small 
number (i.e. one to thirty) of fast particles that are created in a single nuclear explosion or 
star. In these stars, slow 7z-mesons are sometimes observed (Brown et al. 1949), and there is 
further evidence of the presence of neutral particles (Rochester and Butler 1947, Kaplon, 
Peters and Bradt 1949). The Plate shows a photograph of a penetrating shower of six 
particles which was obtained with a counter-controlled Wilson cloud chamber containing 
nitrogen at a pressure of 60 atmospheres. ‘The shower, which has originated in a lead plate 
placed above the cloud chamber, shows two interesting events. Firstly, one of the particles 
(track ABC) can be identified as a slow 7-meson decaying in flight at B, and secondly, a 
particle (track DE) starts in the gas as if it came directly (or nearly so) from the star. 

The velocity v of the particle along AB can be estimated in two ways : (a) from the number 
of short electron branches on AB, (b) from the branch electron near A. The first method 
gives the limiting values of v as 0-46¢ and 0-62c, whilst the second yields 0:5¢ and 0-7e. 
The second value is probably less reliable than the first, and it is not unreasonable to take 
0:-55e¢ as the mean. ‘To estimate the mass WM of the particle, the deflection due to multiple 
scattering along the length AB has been measured (Williams 1940), and the limiting values 
of 0-5e<v<0-7e give 60m,<M<600m,. The particle can therefore be safely assumed 
to be a meson. 

At B there is a single deflection of 23°, although the probability of the scattering of a 
meson through an angle greater than 20° is only 10-4 per cm. of path. Thereafter along 
BC the track is perceptibly thinner and straighter than before and, furthermore, only one 
electron branch is observed in the 3-3 cm. length BC, whereas one per cm. is observed in AB. 
This is consistent with the assumption that the particle velocity of the meson, just before B 
is less than its value just af er, and thus it is reasonable to assume that it has suffered en 
at B. Since the time of flight of the meson in the chamber is only of the order of 10~° sec. 
the event is unlikely to represent the decay, in flight, of a u-meson whose half-life is of ine 
order of 2 10~* sec.; on the other hand, it can well be an example of the decay, in flight, of 
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t t 
E C 
A high-pressure Wilson cloud-chamber photograph of a penetrating shower, showing (qa) track 
ABC, representing the decay of a 7-meson in flight, (b) track DE, a particle starting in the 
*gas as if it came directly from the origin of the shower. 
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Figure 1. Oxidized zinc-blende cleavage face; beam parallel to 
cube-face diagonal. 


Figure 2. Oxidized zinc-blende cleavage face of Figure 1; beam 
along the cube edge. 
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am-intoap-meson. It can be estimated that the velocity of the latter along BC is within the 
range 0-6c to 0-8c, and a simple calculation gives the velocity, Vy, Of the z-meson in the 
frame of reference in which the 7m-meson is at rest to be 0:29e<v,<0:50c; this velocity has 
been found by other methods to be 0:28e (Lattes, Occhialini and Powell 1948). Itcan be 
shown that with these values of uv, the 4-meson would be projected at right angles to the 
original direction of the 7-meson. 'To sum up this much can be said: it is exceedingly 
probable that the track ABC represents the decay, in flight, of a meson whose half-life cannot 
greatly exceed 10~* sec., and that most likely it represents the decay of a 7-meson. 

Turning now to consider the track DE, the particle starting in the gas leaves a track 
similar to the other shower particles, and this means that it is either an electron whose energy 
is of the order of 507mc? or a heavier particle whose velocity is at least 0-8c. It cannot 
represent the track of a proton, since the deflection in this case due to multiple scattering 
would be negligible. The particle is therefore assumed to be either an clectron or a meson. 
If it is an electron, then one explanation of this event would be that it is a Compton recoil 
electron. It can be shown that in this case the energy of the incident photon must be less 
than 100mc? and greater than 50mc?. If its energy is greater than 100me?, it is more likely 
that an electron pair would be produced, whilst if it is less than 50mc? it would only rarely 
happen that the recoil electron would be projected at such a small angle (i.e. 5° or less) to the 
original direction of the photon (Heitler 1935). Now if this event can be interpreted as a 
Compton effect, it is interesting to consider the source of the y-radiation. If it comes from 
the star, it is difficult to understand why there are so few electrons coming from the walls of 
the chamber, which are of steel of thickness 6 cm. It is only in large electron showers that 
recoil electrons of such energy are observed starting in the gas, and then only in small 
numbers representing less than 2° of the total number present. It therefore seems that the 
photon is created either in the steel or in the chamber itself, or, in other words, it must have 
been created in some secondary process in which the other fragments of the star are involved. 
Now no experimental evidence has yet been found for the existence of y-radiation arising 
from the decay of charged 7- or sx-mesons; on the other hand, it has been suggested that 
neutral mesons, if they exist, would, in decaying, give rise to two photons, and on this view 
the event under consideration might be interpreted as due to the scattering of one of these 
photons. While this is a possible interpretation, there are a number of objections to it. 
Firstly, the number of photons produced in this way should be about equal to the number of 
charged mesons in the shower; secondly, the scattering cross section in nitrogen for photons 
of energy 50mc? is small. Finally, the striking feature of the photograph is the fact that the 
particle is projected forwards as if coming directly from the star. On the whole it appears 
that this hypothesis is unsatisfactory, and it may well be that some hitherto unsuspected 
nuclear reaction is taking place, whereby a neutral meson, coming directly from the star, 
is able to acquire a charge without at the same time exchange of energy taking place. If 
there were an exchange of energy, other effects must result, and it would certainly not be 
expected that the recoil particle would be projected forwards. 

The experiments are being continued in an attempt to ascert2in the identity of the recoil 
particle. 

The work was carried out with the aid of a research grant from the Department of 
Scientific and Industrial Research; one of us (T.C.G.) is also indebted to this Department 
for a maintenance allowance. We wish to thank Professor P. M. S. Blackett of the 
University of Manchester and Professor R. M. Davies of the Department of Physics at 
Aberystwyth for their help and encouragement in this work. 


Physics Department, G. R. Evans. 
University College of Wales, T. C. GRIFFITH. 
Aberystwyth. 
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Epitaxial Strain and Disorientation in Crystals Growing 
on Single-Crystal Substrates 


The electron-diffraction patterns, Figures 1 and 2, were obtained by one of us 
(D.M.E.) from a zinc-blende {110} cleavage face after heating in air. The crystal was 
a transparent pale yellow and was raised from room temperature to approximately 520° c. in 
air during 45 minutes, kept at 520° c. for a further 30 minutes, and then allowed to cool in 
the oven for about 10 minutes before removal and transference to the electron-diffraction 
camera. Figure 1 was obtained with the electron beam grazing the surface along the 
cube-face diagonal of the ZnS substrate, and Figure 2 was from the same surface with the 
beam along the ZnS cube edge (for Figures 1 and 2 see Plate opposite p. 297). 

Apart from the pronounced arcing (6°) shown in Figure 1, apparently not observed 
previously, the diffraction patterns confirm the {103} ZnO orientation indicated by the 
experiments of Yamaguti (1935), Aminoff and Broomé (1936, 1938), and Uyeda, Takagi 
and Hagihara (1941) (see also Seifert 1940). This difference is probably due to the higher 
temperature (red heat) of oxidation which they used. 

The initial pattern from the crystal which yielded Figures 1 and 2 showed the high 
perfection and atomically smooth cleavage surface of the crystal; and after etching away 
the ZnO—ZnS surface layer by nitric acid a sharp spot pattern was obtained, showing the 
zinc-blende substrate to be still undistorted at a short distance below the interface. 

Recent work submitted for publication in July 1949 (also Wilman 1950) by one 
of the authors (H.W.) demonstrated the common occurrence of rotational slip as a new kind 
of deformation process in crystals. The above results (Figures 1 and 2) are now explained 
satisfactorily as being due to this kind of deformation, caused by the epitaxial stress between 
the substrate and overgrowth crystals, due to their different lattice constants. Thus the 
strong arcing in Figure 1, but the inappreciable arcing in Figure 2, shows clearly that the 
lattice disorientation is, in this case, simply a rotation about an axis paralle! to the cube-face 
diagonal of the {110} ZnS surface. ‘This could therefore correspond to rotational slip of the 
{103 }-orientated ZnO crystals, or of parts of the ZnS substrate, on the planes which are 
normal to this rotation axis, i.e. {1120} in the ZnO (since (100) ZnO is parallel to (110) 
ZnS) or {110} in the ZnS. In the ZnS, at least, this plane is a direction of especially 
perfect cleavage, and thus according to the previous results it should indeed be a direction 
of easy rotational slip; and it should also be similarly easy on the corresponding {1120} 
planes of the ZnO, which has a structure differing from that of zinc blende only in the 
hexagonal type of superposition of its (0001) layers. 

This constitutes the first clear example of rotational slip caused by epitaxial stresses, 
since we have also verified experimentally that zinc blende does undergo easily a lattice 
rotation on {110} planes. For example, when a cleavage face was abraded along the 
cube-edge direction, a strong (~30°) lattice rotation occurred in the surface region, about 
an axis in the surface and normal to the abrasion direction, corresponding to rotational slip 
on the {110} planes which were parallel to the abrasion direction and normal to the surface. 

_We have similarly demonstrated extensive rotational slip in a copper crystal on {001}, 
{110} and {111} planes, by abrasion of an anodically polished {110} face in directions 
parallel to these planes respectively, which were normal to the surface. In this case there 
is no possibility of interpretation in terms of flexural translational slip, since such slip in 
copper is well authenticated as being along (110) in {111}. 

Rotations of this kind, due to abrasion of {001} cleavage faces along (100), were noted in 
PbS by Germer (1936) and Beukers (1939) and also in NaCl by Raether (1947). Raether 
followed Germer’s interpretation, which entirely neglected the possibility of lateral cohesion 
between the rotated blocks (in these cas<s cleavage blocks) on their common faces which are 
normal to the rotation axis. 

In order to present a possible picture of the way in which the rotational slip arises in 
orientated overgrowth crystals, let us consider Figure 3 (a): ABCD and A’B’C’D’ represent 
the corners of neighbouring rotational-slip lamellae as part of an overgrowth crystal before 
the stress is relieved by its subdivision. The part near CD adjoining the substrate is held 
under tension (or compression) paralle! to CD owing to the atomic forces of the substrate 
on the ov ergrowth crystal, while further from the substrate, i.e. along AB, the lattice 
dimensions tend to their normal values. It thus seems possible that relief of the tension 
(or compression) near CD, by a slip at the overgrowth-substrate interface, may occur 
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unevenly over the contact surface, so that, for example, the lamella ABCD may remain held 
near C and contract (or expand) from D towards C, while A’B’C’D’ may be held near D’ 
and contract (or expand) in the direction C’D’. This motion, oppositely directed in the 
two parts of the crystal, may thus result in a relative rotational slip which in general will not 
be fully reversed as the interface slip proceeds; thus there will be left both a residual 
rotational-slip deformation and a residual but reduced strain near the interface (Figure 3 (d)). 


A B 


Overgrowth 
Crysta 


ee 
— Substrate Crystal 
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Figure 3. Schematic representation of a possible way in which an epitaxially strained overgrowth 
crystal (a) may have its strain partially relieved by a lamellar subdivision by rotational 
slip (6) on the plane ABCD. 


The nearly uniform intensity along the arcs in Figure 1 suggests that, in actual fact, a 
fine-lamellar subdivision of the ZnO crystals occurred, initiated by a sudden release along 
one lamellar region and thereby passing on rotational motion to neighbouring regions as in 
normal rotational slip, but also triggering their partial release at the interface region. 

It is hoped to publish a fuller account of the above experiments in the near future. 


Applied Physical Chemistry Laboratories, D. M. Evans. 
Imperial College, London S.W.7. H. WILMaAN. 
28th December 1949. 
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Corrigendum 


“Time Dependence of Electronic Processes in Dielectrics”, by H. FROnLICcH 
and J. O’Dwyer (Proc. Phys. Soc. A, 1950, 63, 81). 


Equation (8), page 83, should read 
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REVIEWS OF BOOKS 


The Structure of Matter, by F. O. Rice and E. Terier. Pp. xii+361. Ist 
Edition. (New York: John Wiley and Sons; British Agents: Chapman 
and Hall, 1949.) 30s. 


The great successes achieved by quantum theoretical methods in the rationalization and 
understanding of physical and chemical phenomena have raised inevitable difficulties in the 
teaching of physics and, particularly, of chemistry. The phenomena of molecular structure 
and of chemical reactions involving simple systems are now amenable to quantitative 
interpretation in terms of concepts such as energy levels, wave functions, degenerate states, 
bonding orbitals, spin functions, reaction complexes, and so on, but the understanding of 
these concepts requires some knowledge of quantum mechanics, and this in turn requires a 
mathematical equipment well beyond the elementary calculus with which most students of 
physics and chemistry can be assumed to be equipped. 

The work under review aims at providing an introduction to the study of theoretical 
chemistry and related topics while keeping the mathematical formalism to a minimum. As 
an indication of the success achieved, the authors manage to express in a lucid and imagin- 
ative way some of the most fruitful (and difficult) concepts of quantum theory without once 
stating the Schrédinger Equation. 

The book however deals, not with the concepts themselves, but with their application 
to atomic and molecular structure, the chemical bond, atomic and molecular spectra, 
molecular vibrations and related topics. Chapters are also included on nuclear chemistry 
and the state of matter in stars. ‘The presentation is stimulating throughout and the book 
represents a notable contribution to the teaching of quantum chemistry. Students of 
physics and chemistry in the third year of their course, and graduate students alike, will 
find it of great value. 

This book is the first in the Wiley “ Structure of Matter’ Series and sets a very high 
standard both in subject matter and in technical excellence of production. E. H. S. B. 


The Adsorption of Gases on Solids, by A. R. MILLer. Pp. xi+133. Cambridge 


Monographs on Physics, Ist Edition. (Cambridge: University Press). 
[2s. 6d. 


Ten years ago J. K. Roberts wrote a monograph called Some Problems in Adsorption 
giving a live account of his own recent researches, experimental and theoretical, on the 
adsorption of gases (hydrogen and oxygen) as monolayers at the surface of a clean metal 
(tungsten). This monograph was then, is still, and will remain a valuable, instructive and 
faithful first-hand record of how an outstanding physicist plans and directs the exploration of 
a new field. Unfortunately the monograph is out of print. 

The editors instead of reprinting Roberts’ monograph asked his pupil A. R. Miller to 
revise’ it. The result is the book now under review. The ‘ revision’ consisted of 
(a) omission of some bits of theory and substitution of others (the essentially new matter is 
equivalent to rather more than a chapter), (4) some changes in the order of presentation, 
(c) replacement of an accurately descriptive title by a thoroughly misleading one, (d) sub- 
stitution of the reviser’s name for the author’s. 

The reviewer is not at all clear what is gained by thus revising this kind of monograph. 
The product is no longer an authoritative description of how a great researcher works. 
It is certainly not the sort of book that Roberts himself might have written now had he been 
alive, for his outstanding charactetistic was progress, not stagnation. He would have 
discarded the clumsy statistical methods used ten years ago in favour of the simpler and 
more direct methods developed meanwhile and incidentally used elsewhere by Miller him- 
self. Still less can the book be regarded as a comprehensive textbook on the adsorption of 


gases on solids. Perhaps it is merely the book which the editors or the reviser would have 
wished Roberts to write. E. A. G 


‘ 
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Les Spectres Moléculaires, Colloques Internationaux, Paris, May 1947. Pp. 176. 


(Paris: Centre Nationale de la Recherche Scientifique; British Agents: 
H. K. Lewis, 1948). 16s. 


Collected in this book are the 46 papers which were presented at an international 
conference held in Paris in May 1947, and which have since been reprinted either in the 
Journal de Physique or in the Fournal de Chimie Physique. A few papers discuss topics in 
the spectroscopy of diatomic molecules—for example, the demonstration of septet—septet 
transitions in the manganese halides, and the question of interchange of rotational and 
translational energy between AIH and He—but most of the work can be described as the 
application of infra-red and Raman spectroscopy to structural problems of polyatomic 
molecules. Notable here are the studies of crystals which seek information about questions 
as diverse as those of symmetry, dynamics and hydrogen-bonding. Earlier established 
subjects are not neglected, and there are articles on the improvement of vibrational potential 
functions, measurements of absolute intensities and the use of isotopic substitution. The 
progress here is illustrated by the account of the first singlet excited state of benzene, 
which fits so happily with what is known theoretically of the electronic structure of this 
species. 

Although the work of some of the most active laboratories in this field is not represented, 
so that there is, for example, no mention of microwave spectroscopy, the book does provide 
both an interesting picture of the directions in which the subject is moving and a useful key 
to some of the literature—there are about 300 references, many to work done or published 
during the war years. R. F. B. 


Electron Tubes, by the LABORATORIES OF THE RADIO CORPORATION OF AMERICA. 
Volume I (1935-41). Pp. x+475. Volume II (1942-48). Pp. x+454. 
Ist Edition. (Princeton, New Jersey: R.C.A. Review, 1949.) $2.50 
each volume. 


‘These two volumes form the ninth and tenth in the R.C.A. Technical Book Series and 
are the first to be devoted entirely to electron tubes. They are made up of reprints of some 
of the most important papers published by workers at the Radio Corporation of America, 
together with abstracts of some other papers and a complete bibliography of all the publi- 
cations of the Laboratory on this subject. 

The material presented in Volume I, which covers the period 1935-41, is by now fairly 
well known, and the collection here presented will probably be of most interest to workers 
in the subject, partly as an historical account and partly as containing some of the key papers 
on important aspects of tube design. The series of five papers by Thomson, North and 
Harris on ‘“ Fluctuations in Space-charge Currents at moderately high Frequencies ”’ are 
by now classicaland provide a good introduction to the problems of ‘ noise’ in electron tubes. 
The paper by Linder on “‘ Excess-energy Electrons and Electron Motion in high-vacuum 
Tubes ”’ deals with excess energy electrons in magnetrons and the excessive ‘ noise’ pro- 
duced by them. The subject is of present-day interest, not only because magnetrons are 
still ‘ noisy ’ but also because mechanisms of the type described here are being suggested as 
possible causes of ‘ noise ’ radiation from sunspots. The papers on “‘ Beam Power Tubes ”’ 
by Schade and “‘ The Secondary Emission Multiplier”? by Zworykin, Morton and Malter 
are well known. 

Volume II (1942-48) is likely to be of more general interest, not only because the subject- 
matter is newer but also because many radio workers, not immediately concerned with 
problems of vacuum tubes, had not time during the war years to follow the literature on the 
subject. For them it is valuable to have a collected series of reprints of the kind here 
presented. The papers on ‘“‘ Grounded Grid Amplifiers’? (M. C. Jones), “ Excess Noise 
in Cavity Magnetrons”’ (R. L. Sproull), “Maximum Efficiency of Reflex Klystrons # 
(Linder and Sproule), ‘“‘ Beam Deflection Control for Amplifier Tubes ” (Kilgore), “* Small 
Signal Analysis of Travelling-wave Tubes” (Shulman and Heagy) and “ Analysis of a 
Simple Model of Two-beam Growing Wave Tube” (Nergaard) all deal with important 
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fundamental principles. It is of interest to note that work on the : growing wave tube’ 
appears to have started simultaneously and independently at R.C.A., by Haeff at the Naval 
Research Laboratory, and at the Bell Telephone Laboratories. 

Naturally the volumes present a one-sided view of the development of electron tubes, but 
since the R.C.A. Laboratories have played such a great part in this development the papers 
cover most aspects of the subject, with only a few important gaps, amongst which the problem 
of the cavity magnetron is most noticeable. These volumes will find a place in all labora- 
tories dealing with the design of electron tubes, and Volume II can well be read by the 
radio engineer with wide interests, at least until the subject is really well written up in 
textbook form. J. A. RATCLIFFE. 


Reprints of Papers in the ‘Proceedings of 
the Physical Society’ 


In general, reprints of papers published in the Proceedings of the Physical 
Society are not available at the Offices of the Society, but may occasionally be 
obtained from the authors of the papers concerned. ‘The address of the author 
is given at the head of each paper. 

Exceptions to this are the special Lectures delivered before the Society, such 
as the Guthrie Lecture, the Presidential Address, the Charles Chree Lecture and 
the Rutherford Lecture. 
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ABSTRACTS FOR SECTION B 


The Deposition of Ice Crystals on Cooling Surfaces in Low Temperature Plant, by 
G. G. HASELDEN. 


ABSTRACT. A study is made of the conditions governing the deposition of ice crystals 
from a moist air stream on to the inside surface of a cooled tube. The promoting effect of 
impurities such as oil vapour in the air stream is demonstrated. Series of experiments were 
made to establish the effect of air velocity, and of the temperature difference between the 
tube wall and air stream, on the amount of deposition. In the flow range tested (Reynolds’ 
number 11,000 to 27,000) it was observed that increasing the air velocity increased the per- 
centage deposition; at the lower velocities there was no evidence of progressive deposition. 
The effect of temperature difference on the percentage adhesion in the range 7° to 14° c. was 
small, but an experiment at a value of 4:6° c. showed a decrease in the amount of condensed 
snow deposited. In most experiments a glass tube was used but in others use of a brass 
tube gave similar results. 


A Relationship between the Refractive Index and the Infra-Red Threshold of 
Sensitivity for Photoconductors, by 'T. 5. Moss. 


ABSTRACT. It is shown that for photoconductive compounds the long wavelength 
‘ threshold ’ of sensitivity A should be related to the refractive index n of the photoconductor, 
and analysis of the available data shows that for the more highly refractive compounds 
n*/X always has a value close to 77. As a consequence of this relation, an explanation can 
be given for the change of the threshold wavelength with temperature which is observed for 
lead sulphide and similar materials. Suggestions are made as to what materials ought to be 
suitable for photoconductive detectors at longer wavelengths in the infra-red, and some 
confirmation of the ideas is shown by results obtained with cadmium arsenide, the properties 


of which are briefly described. 


Capacitance Measurements on Selenium Rectifiers: Evidence of Anomalous Dis- 
persion, by R. Cooper. 


ABSTRACT. ‘The paper presents evidence showing that, contrary to the generally 
accepted view, the capacity of a selenium rectifier at audio frequencies may vary with the 
frequency of measurement. The magnitude of the variation depends on the metal used as 
counter-electrode. Thus in the frequency range 300—2,000 c/s. little variation occurs with 
antimony or bismuth, but with cadmium there is a considerable decrease in capacitance. 
The effect is considered due to polarization arising from the relative motion in the selenium 
lattice of the negatively charged impurity centres and positive ions which diffuse into the 


selenium from the counter-electrode. 


On Dissociation Processes in certain Gases of High Dielectric Strength, by J. W. 
Warren, W. Hopwoop and J. D. Craccs. 
ABSTRACT. The dissociation processes, due to electron collisions, in certain polyatomic 


gases such as CCl,, CCI,F, are discussed with reference to the unusually high dielectric 
strength of certain of these gases. Possible explanations of the latter are then suggested in 


terms of these dissociation processes, 
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The Static Coefficient of Friction and the Area of Contact, by R. C. PARKER and 
D. Harcu. 


ABSTRACT. An apparatus for measuring the static coefficient of friction has been 
designed in which the forces tending to disturb the contact area are minimized and which 
further enables the contact area to be viewed throughout the experiment. Application of a 
tangential load, less than that required to give limiting friction, has been shown to give rise 
to a growth of both the real and apparent contact areas and to produce a relative displacement 
of the surfaces to an extent governed by the tangential load. An examination of the surface 
topography, by means of an electron microscope, has shown that the earlier type of slip 
ceases because surface scratches extend the boundaries of real contact. Both the force of 
adhesion and the tangential force vary linearly with the real and apparent areas of contact. 
These results are in accord with a cohesive theory of friction, and offer no support for the 
theory on which friction is attributed to the formation of welded junctions. 


The Thermal Etching of Single Crystals of Cadmium, by E. N. pa C. ANDRADE 
and R. F. Y. RANDALL. 


ABSTRACT. The heating of single crystal wires of cadmium in circumstances that permit 
free evaporation produces a development of crystalline planes. When the basal plane makes 
a small angle with the axis of the wire there is a characteristic formation of hexagonal pits, 
the bottoms of which are mirror-like basal planes: when the basal planes make larger 
angles, characteristic elliptical traces appear, reminiscent of those formed by mechanical 
extension. The thermal etching reveals, by the formation of steps, thin laminae existing 
in the unstrained wire : the thickness of these, of the order of 1 jx, is the same as the spacing 
of the glide planes revealed by mechanical extension. The etched figures, in general, indicate 
that the substructure revealed by straining is already present in the unstrained crystal. The 
secondary pyramidal faces revealed agree with the calculations of Stranski, Kaischew and 
Krastanov. 

The mirror-like basal planes are sufficiently large to permit the application of an optical 
method to reveal a scatter in their inclination to the axis. "This method shows that there is a 
variation, of round about 20 minutes, of orientation among the basal planes about the direction 
of growth of the crystals, which were grown by the method of Andrade and Roscoe and 
showed good single crystal properties by the usual criteria. 


The Lattice-Parameters of Clear Crystalline Quartz, by H. D. Ketru. 


ABSTRACT. Some experimental results are presented showing that the calibration of 
Debye-Scherrer x-ray diffraction cameras by means of clear crystalline quartz may lead to 
appreciable errors unless the parameters of the quartz are first determined by a direct 
method. ‘The parameters of a sample of clear Brazilian quartz have been found to be in 
serious disagreement with the accepted values, and variations in the parameters of synthetic 
quartz, with an apparent dependence on growth temperature, have also been observed. 
Taken in conjunction with previous results of other workers, the experiments show that 
quartz parameters vary from sample to sample by amounts as large as 0:01% according 
to the impurity content. The previously proposed quartz calibration standard is rejected. 
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